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ABSTRACT 

The  objective  of  the  STAR  project  is  to  test  and  space 
qualify  a  new  continuous  cycle  cryogenic  refrigeration 
system  for  cooling  of  sensors  and  electronics  which  is 
based  upon  the  newly  discovered  thermoacoust ic  heat  pumping 
effect.  The  new  refrigerator  has  no  sliding  seals,  a  cycle 
frequency  of  about  300  Hz,  and  uses  acoustic  resonance  to 
enhance  the  overall  power  density  and  efficiency.  This 
thesis  is  concerned  specifically  with  the  design  and 
testing  of  the  thermodynamic  element  (or  stack),  which  is 
responsible  for  the  thermo-acoust ic  power  conversion,  and 
the  testing  of  binary  inert  gas  mixtures  as  working  fluids. 
Using  the  refrigerator's  coefficient  of  performance 
relative  to  the  ideal  Carnot  coefficient  of  performance  as 
a  measure  of  efficiency,  we  have  achieved  a  93%  improvement 
over  previous  designs. 
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I.   INTRODUCTION 

A.   BACKGROUND 

This  being  the  second  in  a  series  of  theses,  covering 
different  aspects  of  the  Space  Thermoacoust ic  Refrigerator 
(STAR)  project,  a  great  deal  of  the  introductory  material 
has  already  been  written.  In  order  for  this  thesis  to  be 
complete,  the  introductory  material  that  is  felt  to  be 
necessary  will  be  repeated  here.  With  the  exception  of 
minor  corrections.  Chapter  I  and  sections  A. 1  and  A. 2  of 
Chapter  II  have  been  reproduced  from  the  work  of 
Fitzpatrick  (1988). 

1.   History 

Thermoacoust ics  can  generally  be  described  as  the 
study  of  the  interaction  between  heat  and  sound.  Scientific 
interest  in  this  interaction  is  not  new.  Lord  Rayleigh 
(1878  and  1945)  discussed  various  qualitative  examples  of 
the  production  of  sound  by  heat.  In  one  of  these  examples 
he  describes  investigations  by  Sondhauss  (1850)  of  an 
effect  observed  by  glassblowers :  heating  a  bulb  of  glass  at 
the  end  of  a  hollow  rod  occasionally  led  to  the  generation 
of  sound.  Faraday  (1818)  demonstrated  that  this  effect 
occurred  with  gases  other  than  air.  Sondhauss  proved  that 
the  Vibration  of  the  glass  itself  did  not  generate  the 
sound,   but  he  offered  no  explanation  as  to  what  did.  In  his 
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description  of  Sondhauss'  work,  Rayleigh  stressed  the 
importance  of  the  phase  difference  between  temperature  and 
particle  motion.  Rijke  (1859)  built  an  open  tube  with  a 
wire  mesh  inside.  When  the  mesh  was  heated,  the  tube 
produced  a  sound  of  impressive  intensity.  The  functioning 
of  this  oscillator  is  slightly  different  since  it  requires 
steady  gas  flow  for  its  operation.  A  similar  sound 
production  effect  was  observed  by  Taconis,  et  al.  (1949)  in 
hollow  tubes  immersed  in  liquid  helium  (4.2'^K).  The  Taconis 
oscillations  were  addressed  quantitatively  by  Yazaki,  et 
al.  (1980).  The  work  described  above  deals  primarily  with 
prime  movers.  These  are  devices  that  convert  a  temperature 
gradient  into  sound  energy.  We  call  this  the  classical  era 
of  thermoacoust ics. 

Interest  in  thermoacoust ics  was  renewed  when  the 
idea  of  the  use  of  acoustical  devices  as  refrigerators 
appeared.  Gifford  and  Longsworth  (1966)  described  the 
pumping  of  heat  along  a  surface  caused  by  a  periodic  change 
in  the  pressure  of  the  adjoining  gas.  Such  a  change  can  be 
produced  by  the  oscillations  of  a  sound  wave.  In  their 
experiments  Merkli  and  Thomann  (1974)  explored  the  heating 
and  cooling  effects  on  the  wall  of  a  gas-filled  resonant 
tube.  They  found  that  heat  was  transported  from  a  region 
near  the  velocity  antinode  (or  maximum  value)  of  the 
sound's  standing  wave  to  regions  near  the  adjacent  pressure 
antinodes.   Thermoacoust ic   prime   movers  have  also  received 


recent  attention,  Kempton  (1976)  investigated  the  excess 
noise  of  aeroengines  above  that  predicted  by  theory.  He 
determined  that  the  sound  was  produced  by  unsteady  heat 
transfer.  Each  of  these  latter  three  groups  of 
experimenters  used  some  theory  for  comparison,  but  it  was 
mostly  qualitative.  None  provided  the  complete  quantitative 
theory  that  would  explain  their  experimental  observations. 

The  theoretical  breakthrough  in  the  thermoacoust ics 
field  was  made  by  Nicklaus  Rott  (1969,  1974,  1975  and 
1980).  His  theory  combined  the  basic  principles  from 
physics,  thermodynamics  and  acoustics  to  quantitatively 
describe  the  effects  found  in  both  types  of  thermoacoust ic 
devices:  prime  movers  and  heat  pumps.  Prime  movers  use  a 
temperature  gradient  to  create  sound,  as  in  the  effects  dis- 
cussed by  Sondhauss,  Rayleigh,  Taconis  and  Kempton.  Heat 
pumps,  or  refrigerators,  use  the  oscillating  pressure  of  a 
sound  wave  to  produce  a  temperature  difference,  as 
described  by  Gifford  and  Longsworth  and  Merkli  and  Thomann. 
Rott  described  the  effect  found  by  Merkli  and  Thomann  as 
thermoacoust ic  streaming.  It  is  this  effect  that  makes 
thermoacoust ic  refrigeration  possible. 

Inspired  by  Ceperley's  (1979)  traveling  wave  Ster- 
ling Cycle  heat  engine  and  Rett's  quantitative  theory,  the 
team  of  Wheatley,  Hofler,  Swift,  Migliori  and  Garrett 
(1982,  1983a,  1983b,  1985  and  1986)  developed  a  series  of 
thermoacoust ic     experiments    at    Los    Alamos    National 


Laboratory  in  New  Mexico.  They  investigated  the  basic 
thermoacoust ic  effects  in  both  prime  movers  and 
refrigerators  and  compared  their  experimental  results  to 
Rett's  theory. 

2.   Thermodynamics 

We'll  digress  from  history  here  to  explain  the 
thermodynamic  distinction  between  prime  movers  and 
refrigerators.  This  discussion  follows  Sears  and  Salinger 
(1975).  A  prime  mover  receives  heat  at  a  high  temperature, 
T„  (the  hot  reservoir),  does  work  on  its  surroundings,  and 
rejects  heat  at  a  lower  temperature,  T-,  (the  cold 
reservoir),  as  shown  in  Figure  I-l.a.  The  first  law  of  ther- 
modynamics  tells   us   that   the  heat  flow  (Q„)  from  the  hot 
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reservoir  must  equal  the  work  done  on  the  surroundings  (W) 
plus  the  heat  flow  ( Qp )  to  the  cold  reservoir: 


Qj^-W+Q^ 


(  I-l) 

The  second  law  of  thermodynamics  states  that  in  every 
process  the  entropy  of  the  universe  must  either  remain 
constant  or  increase.  The  entropy  change  for  the  isothermal 
heat  transport  to/from  the  reservoirs  is  defined  as  the 
heat  flow  divided  by  the  reservoir  temperature.  For  a  prime 
mover  the  second  law  can  be  written: 


.S^iAS^   or   Q^/T^iQ^/T^. 


(1-2) 

The   thermal  efficiency  of  the  prime  mover  is  defined  as  the 
ratio  of  the  work  output  to  the  heat  input: 

^=W/Q^=(Q^-Q^)/Qj^=l-(Q^/Q^).  (1-3) 


FIRST   LAV:    Qh  =    W    -    Qr 


SECOND   LAV;      9li  ^    -9^ 


Tu  T. 


EFFICIENCY: 


V   - 


V 


T. 


(a) 


V 


FIRST  LAV:  Qh  =  V  +  Q, 


SECOND  LAV 


Qm     Or 


Th  "  T, 


COEFFICIENT  OF  PERFORAANCE : 


COP  =  -i  ^  ^ 

V 


Th-  Te 


(b) 


Figure  I-l    The  Thermodynamics  of  (a)  Prime  Movers 
and  (b)  Ref r 1 oerat ors 


Comparing  equations  1-2  and  1-3  gives  us  a  limit  on  the 
efficiency  based  on  the  temperatures  of  the  reservoirs: 

7il-(T^/T^).  (1-4) 

This  limit  is  known  as  the  Carnot  efficiency. 

In  a  refrigerator  the  process  is  basically  reversed 
(see  figure  I-l,b.).  Heat  (Q^)  is  removed  from  the  low 
temperature  reservoir,  work  (W)  is  received  from  the 
surroundings,  and  heat  ^Qu^  is  delivered  to  the  higher 
temperature  reservoir.  The  expression  for  the  first  law  is 
the  same  as  for  the  prime  mover.  Equation  I-l.  For  the 
refrigerator  it  means  the  heat  flow  into  the  hot  reservoir 
must  equal  the  heat  flow  from  the  cold  reservoir  plus  the 
work  done  on  the  system.  The  second  law  for  this  case  is: 

(  1-5) 

The  efficiency  parameter  for  the  refrigerator  is  called  the 
coefficient  of  performance  (COP).  It  is  defined  as  the 
ratio   of   the  heat  flow  from  the  cold  reservoir  to  the  work 


'^H-^^C      °^   V^H-Qc^^C 


input  to  the  refrigerator,  or 


COP=Q^/W=Q^/  (  Qj^-Q^  )  iT^/  (  Tj^-T^  )  . 


The    limit    is    known 


as 


(  1-6) 
the   Carnot   coefficient   of 


performance. 

The  work  in  this  thesis  deals  solely  with  refriger- 
ators. Even  though  prime  movers  are  mentioned  from  time  to 
time,  the  primary  focus  of  the  rest  of  our  discussion  will 
be  refrigerators. 


3.   Overall  Efficiency 

According  to  Hofler  (1986),  the  early  thermoacoust ic 
work  at  Los  Alamos  focused  on  experimental  refrigerating 
engines.  The  performance  of  these  engines  fell  short  of 
expectations,  leading  to  simple  experiments  on  basic 
thermoacoust ic  effects  and  a  proof-of-pr inc iple  experiment 
on  thermoacoust ic  refrigeration.  Hofler  then  applied  the 
Rott  theory  to  the  experimental  systems  and  solved  the 
resulting  equations  numerically.  For  his  doctoral 
dissertation  from  the  University  of  California,  San  Diego, 
Hofler  designed  and  constructed  a  completely  functional 
thermoacoust ic  refrigerator.  He  also  made  accurate 
measurements  of  its  thermodynamic  efficiency,  and  used  this 
efficiency  to  make  comparisons  to  the  Rott  theory.  After 
receiving  his  doctorate,  Hofler  came  to  the  Naval 
Postgraduate  School  ( NPS )  as  a  post-doctoral  fellow,  and 
brought  his  prototype  refrigerator  with  him. 

The  purpose  of  this  thesis,  in  conjunction  with  the 
work  of  several  other  students,  is  to  modify  Hofler's 
refrigerator  design  in  order  to  improve  its  overall 
efficiency  and  make  it  suitable  for  space  cryocooler 
applications.  The  basic  design  of  the  thermoacoust ic 
refrigerator  is  shown  schematically  in  Figure  1-2.  The 
driver  (A),  which  produces  the  sound  waves,  is  coupled  to 
the  resonator  (D)  using  bolts,  via  a  reducer  cone  (B)  and 
bellows   (C).   Inside   the   resonator   is  a  stack  of  plastic 
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Figure  1-2.  Basic  Thermoacoust ic  Refrigerator 
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plates  (E)  and  their  heat  exchangers  (F).  The  heat 
exchangers  allow  heat  to  be  removed  from  the  hot  end  of  the 
stack  and  applied  to  the  cold  end.  It  is  the  interaction 
between  the  sound  waves  and  the  plastic  plates  that 
produces  a  temperature  difference  across  the  plates  and 
allows  heat  to  be  pumped  from  the  cold  heat  exchanger  to 
the  hot  heat  exchanger.  A  brief  description  of  efficiency 
would  be  useful  here.  Figure  1-3  shows  a  diagram  of  the 
chain  of  energy  transformations  which  occur  in  a 
thermoacoust ic  refrigerator.  There  are  three  conjugate 
variable  pairs  (V  and  I,  P  and  U,  and  LT  and  Qp ) ,  and  two 
transformations  (electrical-to-acoustical  and  acoustical-to- 
thermal  ) , 

We  start  with  electrical  input  power  to  the  driver, 
which  can  be  calculated  (using  r.ra.s.  values)  with: 

(  1-7) 

where  I  is  the  current  into  the  e lectroacoust ic  driver,  V 
is  the  voltage  across  the  driver  terminals,  ^  is  the  phase 
angle  between  I  and  V,  and  cosw  is  known  as  the  "power 
factor."  An  e lectrodynamic  driver  converts  the  electrical 
power  to  acoustic  power  (the  first  transformation). 
Analogous  to  the  previous  definition  of  the  electrical 
input  power,  the  acoustic  power  is  given  by: 

W=PUcos*,  (1-8) 

where  P  is  the  acoustic  pressure,  U  is  the  volumetric 
velocity,   which  equals  the  product  of  particle  velocity  and 


P  ,   =IVcose, 
e  lee 


DRIVER 

REFRIGERATOR 

ELECTRIC 

V.     I 

p 

elec 

VIcosG 

ELECTRIC  \ 

ACOUSTIC 

P.    U 
W  = 

PUCQS0 

ACOUSTIC\ 

THERAAL 
AT.    Q, 

TO         / 
ACOUSTIC 

TO        / 

THERAAL 

ELECTROACOUSTIC 
EFFICIENCY 


REFRIGERATOR 

COEFFICIENT  OF 

PERFORAANCE 


A 


W 


EA   P, 

si  ec 


COP 


a, 


w 


OVERALL  THERAOACOUSTIC  EFFICIENCY 


n 


Q, 


TA     p 

e  1  ec 


n     •  COP 

EA 


Figure  1-3.  Chain  of  Energy  Transf oraiat ions 
in  a  Ther»oacou3tic  Refrigerator 


0 


resonator  cross-sectional  area,  and  i>  is  the  phase  angle 
between  them.  The  efficiency  of  this  e lectroacoust ic  trans- 
formation is  given  by: 

7^  =W/P  ,   .  (1-9) 

' EA    elec 

In  the  second  transformation,  the  refrigerator  con- 
verts the  acoustic  power  to  a  temperature  gradient  (  T)  and 
heat  flow  (Q).  As  discussed  previously,  the  efficiency  for 
this  transition  is  given  by  the  coefficient  of  performance: 

COP=Q^/W.  (1-1(2)) 

The  overall  thermoacoust ic  efficiency  is  therefore  the 
product  of  the  efficiencies  for  the  two  transformations: 

Improvements  in  the  overall  refrigerator  efficiency 
can  thus  be  divided  into  two  distinct,  though  interrelated, 
efforts.  The  first  is  the  subject  of  this  thesis: 
modifications  to  improve  the  coefficient  of  performance. 
The  second  effort  involves  modifications  to  the  driver  to 
optimize  the  electroacoust ic  efficiency  and  is  addressed  by 
LT  Michele  Fitzpatrick,  U.S.  Coast  Guard,  in  her  master's 
thes  is. 

B.   SPACE  THERMOACOUSTIC  REFRIGERATOR  (STAR) 
1 .   Hot  i vat  ion 

In  addition  to  improving  overall  efficiency,  our 
experimental  refrigerator  will  be  designed  and  built  with 
the   intention   of  flying  it  on  the  space  shuttle  as  part  of 


the  National  Aeronautics  and  Space  Administration's  (NASA) 
Get  Away  Special  (GAS)  program.  As  of  3  February  1988, 
there  is  a  Memorandum  of  Agreement  between  the  Naval  Post- 
graduate School  (NFS)  and  the  Air  Force  which  provides 
funding  for  the  space  flight  and  which  assigns  the  NASA 
payload  number  G-337  to  this  project.  The  reason  for 
testing  the  thermoacoust ic  refrigerator  in  space  is  that 
the  immediate  applications  of  the  STAE  are  expected  to  be 
the  cooling  of  electronics,  high  T^  superconductors  and 
infrared  detectors  in  space. 

There  are  currently  two  primary  methods  for  cooling 
satellite  systems  in  space:  (1)  evaporation  of  expendable 
cryogens  (liquid  helium,  nitrogen,  ammonia-methane  and 
solid  hydrogen,  etc.),  and  (2)  closed  cycle  refrigerators 
(Stirling  cycle,  Vuilleumier  cycle,  etc.),  which  involve 
large  reciprocating  piston  masses  operating  at  low 
frequency.  [Walker  (1983)  and  Smith,  et  al.  (1984)3  The 
disadvantages  of  these  current  cooling  methods  are  their 
short  lifetimes  (expendable  cryogens)  and  high  vibration 
levels  and  low  reliability  (closed  cycle  refrigerators). 
The  advantages  of  the  thermoacoust ic  refrigerator  over 
these  and  other  systems  lie  in  its  simplicity  and 
reliability.  The  STAR  has  no  sliding  seals,  relatively 
efficient  performance,  low  vibration  levels  and  (hopefully) 
a  longer  life  span. 


The  thermoacoust ic  refrigerator  needs  to  be 
surrounded  by  a  vacuum  and  insulating  material  to  minimize 
heat  loss  due  to  thermal  conduction  and  radiation.  The 
vacuum  of  space  will  provide  the  best  insulation  to 
determine  the  refrigerator's  ultimate  efficiency.  The 
absence  of  gravity  removes  the  possibility  of  thermal 
convection  as  an  additional  nuisance  heat  transport  effect. 

2.  Get  Away  Special  (GAS)  Program 

NASA's  GAS  program  allows  small,  self-contained  pay- 
loads  to  fly  on  the  Space  Shuttle  in  GAS  canisters  for 
relatively  low  cost  ($10,000)  (Get  Away  Special  Team, 
1984).  The  GAS  can  is  five  cubic  feet  in  volume  and  can 
house  a  payload  of  up  to  200  pounds.  Each  payload  must 
contain  its  own  electrical  power,  control,  data  acquisition 
and  storage  facilities,  etc.  The  Shuttle's  astronauts  will 
turn  a  switch  on  or  off  at  designated  times  during  the 
flight,  but  are  otherwise  not  involved  with  the  experiment. 

3.  Shared  Subsystems  from  NASA  G-313 

Another  group  of  NFS  students  and  faculty  (Boyd,  et 
al.,  1987)  have  taken  advantage  of  the  GAS  program  to 
measure  the  resonant  acoustic  modes  of  the  shuttle  payload 
bay  and  the  ambient  acoustic  environment  produced  as  a 
result  of  main  engine  and  booster  operation  during  launch. 
Their  experiment  is  titled  "The  Space  Shuttle  Cargo  Bay 
Vibroacoust ics  Experiment"  and  is  designated  by  NASA  as 
payload   G-313.   Several   subsystems  that  were  developed  for 


NASA  G-313  Will  be  used  by  NASA  G-337  (STAR).  A  schematic 
of  the  STAR  in  its  GAS  can  is  shown  in  Figure  1-4.  One  of 
the  systems  borrowed  from  NASA  G-313  is  the  controller 
system  used  to  run  the  experiment  and  record  the  data.  The 
recorder  system  consists  of  the  INTEL  model  BPK  5V75 
magnetic  bubble  memory  module  and  an  NSC  800  micro- 
processor-based controller.  Two  other  NPS  students,  LT 
Charles  B.  Cameron,  USN  and  CRT  Ronald  Byrnes,  USA,  will  be 
designing  the  analog  electronics  and  software  to  integrate 
these  systems  into  the  STAR  experiment  for  their  master's 
theses.  A  block  diagram  of  the  electronics  is  shown  in 
Figure  1-5. 

Another  borrowed  system  is  the  power  supply,  which 
consists  of  Gates  brand  lead  acid  battery  cells  (five 
ampere-hour,  two  volts  each).  These  gelled  electrolyte 
batteries  are  ideal  for  the  STAR  due  to  their  high  power 
density,  low  cost  and  the  low  level  of  outgassing  during 
discharge  cycles.  NASA  G-313  used  a  one  layer  battery  of 
68  cells  providing  880  watt-hours  of  energy  and  weighing 
about  80  pounds  (including  the  cells'  support  structure). 
We  will  be  using  two  battery  layers  with  as  many  cells  as 
we  can,  up  to  136.  This  will  give  us  a  minimum  of  680  watt- 
hours  and  a  maximum  of  1,360  watt-hours  of  available 
electrical  energy,  and  a  total  battery  weight  between  80 
and  160  pounds.  The  number  of  cells  we  can  use  will  depend 
on   the   total   weight   of  the  driver-resonator  assembly  and 
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its  auxilliary  equipment  (vacuum  can,  gas  reservoir,  etc.) 
and  electronics,  as  we  have  to  stay  within  the  200  pound 
limit. 

The  use  of  the  GAS  can  imposes  certain  restric- 
tions. Since  we  are  using  batteries  to  supply  the  power, 
the  refrigerator  has  to  be  energy  efficient.  Also,  the  GAS 
can  setup  requires  the  STAK  to  be  compact  and  lightweight. 
Figure  1-6  shows  a  photograph  of  Hofler's  prototype 
refrigerator.  This  setup,  which  includes  parts  that  will  be 
unnecessary  in  space,  is  approximately  six  feet  in  height. 
In  comparison,  the  maximum  payload  height  for  the  GAS  can 
is  28.25  inches,  or  less  than  2.5  feet.  These 
considerations  played  a  major  role  on  the  choice  of 
equipment  for  and  the  design  of  the  STAE. 

C.   SCOPE 

This  thesis  is  primarily  experimental  in  nature.  The 
purpose  being  to  improve  the  thermodynamic  efficiency  of 
Hofler's  prototype  refrigerator.  Theory  is  presented  only 
as  necessary  to  provide  the  reader  with  an  understanding  of 
the  thermoacoust ic  refrigeration  process  and  to  explain  why 
certain  modifications  were  made.  For  this  reason,  no  direct 
comparison  of  overall  refrigerator  performance  to  theory  is 
presented.  Hofler  (1986)  presents  a  quantitative  comparison 
of  the  refrigerator's  performance  to  the  detailed  Rott 
theory. 
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Figure  1-6.  Photograph  of  Hofler*s  (1986) 
Prototype  Refrigerator 
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Chapter  II.  begins  with  a  fairly  in  depth  discussion  of 
thermoacoust ic  theory,  to  provide  the  above  mentioned 
understanding  and  to  illustrate  possible  areas  for  the 
improvement  of  efficiency.  This  is  followed  by  the  details 
of  how  and  why  specific  modifications,  such  as  gas  mixtures 
and  stack  designs,  were  made. 

Chapter  III.  provides  all  of  the  information  necessary 
for  the  construction  of  a  thermoacoust ic  stack. 

Chapter  IV.  provides  a  discussion  of  the  overall 
experimental  apparatus  and  procedures.  We  have  basically 
used  the  prototype  system  of  Hofler  (1986)  for  our 
experiments.  The  material  in  sections  A  and  B  of  this 
chapter  are  derived  primarily  from  Hofler's  doctoral 
dissertation  (1986),  with  emphasis  on  the  few  areas  that  we 
have  changed. 

Chapter  V.  presents  the  methods  used  for  analysis  of  our 
data  and  the  overall  results  of  our  experiments. 

Chapter  VI.  provides  conclusions  and  recommendations  for 
further  development. 


II.   THEORY 

A.   THERMOACOUSTIC  THEORY 

As   stated   in   Chapter  I.,  the  first  two  sections  (A.l. 
and   A. 2)   of   this   Chapter  are  essentially  reproduced  from 
the  master's  thesis  of  LT  Michele  Fitzpatrick  (1988). 

Thermoacoust ic  theory  has  been  developed  in  detail  by 
Rott  (1969,  1974,  1975  and  1980)  and  adapted  to  the  thermo- 
acoustic  refrigerator  by  Wheatley  et  al.  (1982,  1983a, 
1983b,  1985  and  1986),  Wheatley  and  Cox  (1985),  Hofler 
( 1986)  and  Swift  ( 1989) . 

The  space  thermoacoust ic  refrigerator  (STAR)  basically 
consists  of  an  acoustic  driver  producing  sound  waves  in  a 
resonant  tube  (see  Figure  1-2).  This  resonator  is  filled 
with  a  mixture  of  helium-xenon  gas  (~12.5  atomic  %  xenon) 
pressurized  to  ten  atmospheres,  and  contains  a  stack  of 
plastic  plates.  This  section  of  the  Chapter  will  present  a 
qualitative  model  for  the  thermoacoust ic  heat  pumping 
process  followed  by  a  quantitative  development  for  heat  and 
work  flow  at  a  single  plate.  This  is  followed  by  a  more 
rigorous  development  for  the  heat  flow,  work  flow  and 
coefficient  of  performance  for  a  stack  of  plates  that  are 
much  shorter  than  one  quarter  wavelength.  Section  A. 
concludes   with   a   discussion   of   some  of  the  changes  that 
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could    improve    the   coefficient   of   performance   of   the 
refrigerator . 

1.   Qualitative  Picture 

Consider  a  parcel  of  gas  that  moves  back  and  forth 
along  one  of  the  plates  at  the  acoustic  frequency  (see 
Figure  II-l),  As  it  moves,  the  gas  parcel  will  experience 
changes  in  temperature  and  volume.  Part  of  the  temperature 
changes  come  from  the  adiabatic  compression  and  expansion 
of  the  gas  by  the  sound  pressure,  and  part  as  a  consequence 
of  the  local  temperature  of  the  plate  itself.  A  temperature 
gradient  may  develop  along  the  plate  as  a  result  of  the 
operation  of  the  refrigerator.  The  temperature  and  volume 
changes  can  be  described  by  six  separate  steps  (the  changes 
in  volume  are  indicated  by  changes  in  the  size  of  the 
square  parcel  in  the  figure). 

Assume  the  plate  is  at  a  uniform  temperature,  T.  The 
parcel  starts  at  position  1  (X=0)  with  a  temperature  of  T. 
The  sound  wave  moves  the  parcel  to  the  right  to  position  2. 
The  parcel  has  undergone  adiabatic  compression  and  its 
temperature  is  now  T  .  Since  the  temperature  of  the  gas  is 
higher  than  that  of  the  plate,  heat  (Q)  will  flow  from  the 
gas   parcel   to   the  plate.  The  parcel's  volume  is  decreased 


and   its   temperature   lowered   to   T 


This  heat  flow  also 


causes  the  plate's  temperature  to  increase  at  the  position 
X  .  Parcel  position  3  is  actually  in  the  same  location  as 
position   2,   but   is  displaced  vertically  in  the  figure  for 
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clarity.  The  parcel  now  moves  to  position  4  and  expands 
adiabat ical ly  to  a  new  temperature  of  T  ,  A  repeat  of  this 
process  puts  the  parcel  in  position  5  with  a  temperature  of 
T  .  Now  the  temperature  of  the  gas  is  lower  than  that  of 
the  plate  and  heat  (QM  flows  from  the  plate  to  the  gas 
parcel,  expanding  the  parcel  and  raising  its  temperature  to 
T  .  This  heat  flow  causes  the  plate's  temperature  to 
decrease  at  position  X  .  The  parcel  moves  to  the  right, 
again  under  adiabatic  compression,  and  we  are  back  where  we 
started,  in  position  1  with  temperature  T  (adapted  from 
Wheatley,  et'al.,  1985). 

This  gas  parcel  cycle  is  repeated  at  the  resonator's 
operating  frequency.  For  the  STAR  this  frequency  is  between 
approximately  250  and  600  cycles  per  second,  depending  on 
the  gas  types  and  mixture  (specifically  the  speed  of  sound 
in  the  gas)  and  the  length  of  the  resonator. 

Notice  that  the  temperature  of  the  gas  parcel  at 
position  X=0  is  different  depending  on  which  direction  the 
gas  parcel  is  moving.  It  is  this  phase  shift  in  temperature 
relative  to  motion  that  produces  the  thermoacoust ic  effect, 
as  we  will  show  in  the  calculations  that  follow. 

It  is  also  important  to  note  that  in  order  for  heat 
to  flow  between  the  gas  parcel  and  the  plate,  the  parcel 
must  be  vertically  located  within  about  one  thermal 
penetration  depth  of  the  plate.  The  thermal  penetration 
depth   (&«)   is   the   distance  that  heat  can  diffuse  through 
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the   fluid   during   a   time  1/tj  ,       where   u^   is  the  acoustic 

angular  frequency. 

Wheatley,   et   al,   (1986)   describe   thermoacoust ic 

engines   as   consisting  of  long  trains  of  these  gas  parcels, 

all   about   a   thermal  penetration  depth  from  the  plate.  The 

parcels   draw   heat   from   the  plate  at  one  extreme  of  their 

oscillatory   motion   and   deposit  heat  at  the  other  extreme. 

Adjacent   heat   flows   cancel,   except   at   the   ends  of  the 

plate.   The   net   result   is   that   an  amount  of  heat,  Q,  is 

passed   from  one  end  of  the  plate  to  the  other.  It  should  be 

noted  that  the  above  discussion  assumes  zero  viscosity. 

2.   Single  Plate 

Swift  (1989)  illustrated  the  basic  principles  of  the 

thermoacoust ic   engine   by   examining   a  simple  example  of  a 

single   plate   in   a   gas-filled   tube   such  as  we  described 

above.   In   this   example   the   acoustic   and   thermodynamic 

effects   are   nearly  distinct.  In  the  absence  of  the  plastic 

plate,   the   plane   standing   wave  sustained  in  the  tube  has 

well    known   acoustic   properties.   It's   pressure   and   x- 

direction   velocity  (along  the  length  of  the  tube)  are  given 

as  (using  peak  values): 

p  =P^sin (x/x)  ( I I-l) 

1   A 

and     u  =i(P./e  a)cos(x/^)  (II-2) 

1     A  ^  m 


where  P.  is  the  pressure  amplitude  at  the  pressure 
antinode,  ^  =  V2T=a/u^  is  the  reduced  wavelength,  a  is  the 
wavelength,   a   is   the   speed  of  sound  in  the  gas,  u;  is  the 
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angular   frequency   and   f    is  the  mean  density  of  the  gas. 

The   subscript  1  indicates  the  first  order  time  oscillation, 

and   the   subscript   m   indicates   the  mean  value.  The  total 

pressure  can  be  given  as: 

P=P„+Pie^"'^.  (II-3) 

A   Similar   expression   may   be   used  to  represent  the  total 

velocity.   Following   Swift,   we   will  assume  that  the  first 

order   in   the   acoustic   amplitude   is   sufficient   for  all 

acoustic   and   thermodynamic   variables  (pressure,  velocity, 

temperature,  density  and  entropy). 

The   sound  wave  is  adiabatic  and  has  an  oscillatory 

temperature  given  by: 

T  =(T  P/P  c  )p,,  ( II-4) 

1    m   '  m  p  -^1 

where   P=-(df/dT)  /p        is   the   isobaric   thermal   expansion 

coefficient   and  c   is  the  isobaric  (constant  pressure)  heat 

P  ^ 

capacity   per   unit   mass.  Note  that  T.  and  p.  are  in  phase. 
For  ideal  gases,  thermodynamics  tells  us  that: 


T  ^/o   c    =(y-l)T  /yp  , 
m   5m  p       m   ^m 


( II-5) 


where   ^ =c  /c    is   the   ratio   of   isobaric   to   isochoric 
P   V 

specific  heats.  ^  equals  5/3  for  monatomic  gases  and  is 
smaller  but  still  greater  than  one  for  other  gases. 
Combining  Equations  II-4  and  II-5  gives  us: 

(1 1-6) 

The  introduction  of  the  plastic  plate  into  the 
standing  wave  modifies  the  original,  unperturbed  tempera- 
ture  oscillations.   This   modification   is   due  to  the  heat 


T,/T  =C (y-l)/y]p,/p 
1   m  ^1  *^m 
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flow  between  the  gas  and  the  plate,  as  described 
previously.  The  temperature  is  modified  in  both  magnitude 
and  phase  for  gas  about  a  thermal  penetration  depth  away 
from  the  plate.  According  to  Swift,  this  results  in  two 
important  effects:  l)  a  time  averaged  heat  flow  near  the 
surface  of  the  plate,  along  the  direction  of  acoustic 
vibration,  and  2)  the  generation  or  absorption  of  real 
acoustic  power  near  the  surface  of  the  plate.  In 
refrigerators  the  effect  is  an  absorption  of  acoustic 
power.  In  prime  movers  the  acoustic  power  is  generated  by 
the  temperature  gradient  present  in  the  plate. 

There  are  several  assumptions  we  will  make  to 
simplify  calculations,  again  following  Swift.  We  assume 
that  the  plate  is  short  enough  compared  to  a  reduced 
wavelength  (hx«X)  and  far  enough  from  both  velocity  and 
pressure  nodes  that  p.  and  u,  can  be  considered  uniform 
over  the  entire  plate.  We  assume  the  gas  has  zero 
viscosity,  so  that  u.  does  not  depend  on  y  (we  had  already 
made  this  assumption  by  looking  at  u.  in  the  x-direction 
only).  We  assume  that  the  plate  has  a  large  enough  heat 
capacity  per  unit  area  that  its  temperature  does  not  change 
appreciably  at  the  acoustic  frequency.  We  assume  that  the 
plate  has  a  mean  temperature  gradient  in  the  x-direction 
^T  .  Finally  we  neglect  the  plate's  and  gas'  thermal 
conductivity  in  the  x-direction. 
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Applying   these   assumptions,   we   see  that  the  mean 

fluid   temperature  (T  (x))  is  the  same  as  that  of  the  plate. 

m 

Swift  calculated  the  oscillating  fluid  temperature  using 
the  general  equation  of  heat  transfer.  He  kept  only  first- 
order  terms,  neglected  thermal  conduction  along  x  and 
applied  the  boundary  condition,  T.(0)=0,  imposed  by  the 
plate.  The  resulting  equation  is*. 

T,(y)=C(T  P/P  c  )p,-(VT  /w)u,]Cl-expC-( 1+i )y/&w]].  (II-7) 
1-^      m'rap^l     m     1       ^       y" 

To   interpret  this  equation,  we  look  at  it  in  the  limit  that 

the   gas   is   far   enough   from   the   plate   (y»iK)   to  make 

negligible  thermal  contact  with  the  plate.  This  gives: 

T,(y)~C(T  ^/<?  c    )p,-(VT  /w)u,].  (II-8) 

1       m^mpl    m     l 

The    first    term   in   Equation   II-8   is   the   temperature 

oscillation    due    to    the    adiabatic    expansions    and 

compressions   of   the   gas   (see   Equation  II-4).  The  second 

term   comes   from   the  mean  temperature  gradient  in  the  gas. 

As   the   gas   oscillates  along  x  with  displacement  amplitude 

u./>*^,   the   temperature  at  a  given  point  in  space  oscillates 

by   an   amount  VT  \i./^      even   if  the  temperature  of  a  given 
■^  ml  CO 

piece  of  fluid  remains  constant.  The  actual  temperature 
oscillations  are  a  linear  superposition  of  these  two 
effects. 

The  y  dependent  part  of  Equation  II-7  is  complex.  It 
approaches  1  for  y»&Ki  sis  given  previously.  It  approaches 
zero  for  y<-*K»  where  the  plate  imposes  the  condition  T  =0. 
Host   importantly,   for   y~&K,  its  magnitude  is  still  on  the 
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order  of  1,  but  it  has  a  substantial  imaginary  part.  This 
phase  shift  in  the  oscillating  temperature  of  the  standing 
wave  at  y~&K»  ^^®  "^o  '^^^  thermal  presence  of  the  plate,  is 
an  important  result  because  it  leads  directly  to  the  time 
averaged  heat  flow  in  the  x-direction.  This  is  the  same 
phase  shift  that  was  emphasized  in  the  description  of  the 
gas  parcel  motion. 

a.   Heat  Flow 

Swift  argues  that  since  we  are  neglecting  ordi- 
nary thermal  conductivity  in  the  x  direction,  the  only  way 
heat  can  be  transported  along  x  is  by  the  hydrodynamic 
transport   of   entropy,   carried  by  the  oscillatory  velocity 


u 


1 


( II-9) 


q^=T  f  s ,u,=( 1/2) P   c    ImCT,]u,. 
^2m'mll       ^mp     1   1 

The  line  above  the  quantity  s.u,  indicates  the  time-average 
of  the  first  order  entropy  and  velocity  product.  The  heat 
flow  is  a  second  order  quantity,  signified  by  the  subscript 
2, 

The  total  heat  flow  Q^  along  the  plate,  in  the  x 
direction,  is  found  by  integrating  q^  over  the  y-z  plane: 

( 11-10) 
where   TT  is   the   perimeter   of  the  plate  in  the  y-z  plane. 


0  =JT 
^2      0J 


q2dy 


Substituting  for  T.  and  performing  the  integration  gives: 

Q2=-(  1/4)176^  (T^?)p^u^(  r-1)  (11-11  ) 

where   "115^   is   the  thermodynamical  ly  active  area  in  a  plane 
perpendicular   to   the   longitudinal  acoustic  motion,  Tmp  is 
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the   heat   parameter   of   the   gas,   and  r  is  defined  as  the 

ratio    of   actual   temperature   gradient   to   the   critical 

temperature   gradient   (VT/vT   ,^).  The  critical  temperature 
•^  cr  1 1 

gradient  occurs  when  the  temperature  change  along  the  plate 
just  matches  the  temperature  change  due  to  the  adiabatic 
compression  of  the  gas,  and  no  heat  flows  between  the  gas 
and  the  plate.  This  is  the  boundary  between  the 
refrigerator  and  prime  mover  functions  of  thermoacoust ic 
engines. 

Equation  11-11  shows  that  when  r<l,  heat  flows 
up  the  temperature  gradient  from  cold  to  hot  and  work 
(acoustic  power)  is  absorbed,  as  for  a  refrigerator.  When 
r=l  there  is  no  heat  flow.  When  r>l,  heat  flows  down  the 
temperature  gradient  from  hot  to  cold  and  acoustic  power  is 
produced,  as  for  a  prime  mover. 

Note  that  the  total  heat  flow  is  proportional  to 

the   area  TT&k  »  and  to  T  P  (which  equals  1  for  ideal  gases). 

m 

It   is   also  proportional  to  the  product  PiU.,  and  so  equals 

zero   if   the  plate  is  at  either  a  pressure  or  velocity  node 

of   the   standing   wave.   The   maximum   value  of  the  product 

occurs   halfway  between  the  nodes.  Finally,  the  heat  flow  is 

proportional   to   the   temperature   gradient  factor  r-l.  For 

vT    >\7T   .^,   r-l   is   greater  than  zero  and  the  heat  flow  is 
m    cr it  ^ 

toward   the  pressure  node,  while  for  VT  <VT   ,^,  r-l  is  less 

m    cr  i  t 

than   zero  and  the  heat  flow  is  away  from  the  pressure  node. 

If   suitable  heat  exchangers,  at  T„  and  T^,  are  installed  at 

n       L 
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the   ends   of   the   plate  (with  T„-T^=  T  Ax),  this  heat  flow 
carries  heat  from  one  exchanger  to  the  other. 

The  heat  flow  is  small  under  ordinary 
Circumstances.  However,  in  closed  resonators  it  is  possible 
to  achieve  sound  amplitudes  many  orders  of  magnitude  higher 
than  those  of  ordinary  conversation.  Since  Q^  is 
proportional  to  Pa^'  ^^^  since  in  practical  acoustic 
engines  the  entire  cross  section  of  the  standing  wave  is 
filled  with  plates  (spaced  roughly  46k  apart),  very  high 
heat  flows  and/or  a  large  T„-T-  may  be  achieved, 
b.   Work  Flow 

The  work  flow  (i.e.  acoustic  power)  is  given  by 
the  work  per  cycle  times  the  rate  at  which  that  work  occurs 
(the  acoustic  frequency  f).  From  thermodynamics,  the 
average  acoustic  power  produced  per  unit  volume  is: 


w  =-(u//^  )  ip^  =-(  l/2)-'Pp  ImCT^].  (11-12) 
The  gas  about  a  thermal  penetration  depth  away 
from  the  plate  "breathes",  because  of  thermal  expansion  and 
contraction,  with  the  right  time  phase  with  respect  to 
oscillating  pressure  to  do  (or  absorb)  net  work.  This  is 
exactly  the  same  gas  that  we  have  seen  is  responsible  for 
the  heat  flow.  Gas  elsewhere  is  ineffective  in  doing  (or 
absorbing)  work.  The  density  oscillations  for  y^&K  ^^^  for 
y»6K  are  in  phase  with  the  pressure  oscillations,  and  hence 
do  (or  absorb)  no  net  work. 


30 


The   total  acoustic  power  produced,  W^,  is  found 

by  integrating  u^   over  all  space,  as  with  the  heat  flow: 

W^=(  l/4)Tr&H  (^-1)  (p,  V?  a)(r-l)iix/3(,       (11-13) 
z  I'm 

where  all  the  terms  have  been  defined  previously. 

The  acoustic  power  is  proportional  to  the  volume 

TTShAx   of   a  fluid  that  is  about  a  thermal  penetration  depth 

from   the   plate.   It   is   proportional   to   Pi^i  ^'^^   so  is 

quadratic   in   the  acoustic  amplitude  (as  was  the  heat  flow) 

and   vanishes   at  pressure  nodes.  Finally  U^    is  proportional 

to   r-l,  the  same  temperature  gradient  factor  as  appeared  in 

Q^      (Equation   11-11).  When  VT  =7T   ,^,  r-l=0,  and  there  are 
z  m    cr  i  t 

no   temperature   oscillations  in  the  fluid,  other  than  those 

due   to   adiabatic   processes,   and   no   acoustic   power   is 

absorbed   or   generated.   For  7T  >vT   ,^,  r-l>0  and  acoustic 

m    cr  it 

power   is   produced   near   the   plate.   Whether   this   power 

increases   the   amplitude   of  the  standing  wave,  is  radiated 

away   to   infinity,   is  simply  absorbed,  or  flows  through  an 

acoustic-to-electric   transducer   to  generate  electric  power 

depends   on   details   of   the   resonator,   not   on  the  plate 

itself   or  the  standing  wave  near  the  plate.  For  vT  <vT      .^, 

^  ^  m    cr  i t 

r-l<0  and  acoustic  power  is  absorbed  near  the  plate.  For  a 
tube  without  plates,  which  has  a  diameter  less  than  the 
wavelength,  at  constant  temperature  (vT=0),  this  work  flow 
is  responsible  for  the  ordinary  thermal  attenuation  of  a 
sound  wave  (Kinsler,  et  al.,  1982). 


31 


c.   Efficiency 

We   combine  Equations  11-11  and  11-13  to  get  the 

efficiency   of    the    plate   with   no   viscous   losses   or 

longitudinal  conduction  losses: 

'^=W-/Q^=[(^-1)/(T  P)]('^Axp,  )/(  f  a^u,).     (11-14) 

Since  u  =u^sin  (x/x) ,  and  p-=   au^cos  (x/ 5f)  ,  we  get: 

^=C(!^-l)/(Tmp)](  Ax/(Xtan(x/5t)  )  )  .        (11-15) 

For  x«x,  this  reduces  to: 

'7=C(^-l)/(TmP)](Ax/x),  (11-16) 

and    7  =  7     ^/r.  (11-17) 

^  (carnot 

This  efficiency  approaches  the  Carnot  efficiency 
as  the  power  output  and  heat  transfer  rates  approach  zero. 

We   can   make   a   similar   calculation   for   the 

refrigerator    mode    of    operation,    where   the   relevant 

efficiency   is   the   coefficient  of  performance  (COP).  Given 

that  C0P=Q2/W2,  we  find  that: 

COP=r(COP     ^).  (11-18) 

carnot 

3.   Short  Engine 

Having  developed  the  simplified  case  of  a  single 
plate  we  now  turn  our  attention  to  a  more  realistic  model, 
with  an  entire  stack  of  plates.  Since  Swift  (1989)  has 
incorporated  the  thermoacoust ic  theory  of  Rott  (1989,  1973, 
1974,  1975  and  1980),  we  will  once  again  follow  his 
development  closely.  The  assumptions  that  the  length  of  the 
stack  is  much  shorter  than  a  wavelength  and  that  the 
temperature    spanned    is    much   less   than   the   absolute 
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temperature   will   be  retained.  However,  this  more  realistic 

model   includes   viscosity,  longitudinal  thermal  conduction, 

finite   (instead   of   infinite)   plate   heat   capacity  and  a 

stack   of   plates   (instead   of   just   one).  The  goal  in  our 

development   is  to  include  these  effects,  which  produce  many 

complications,   in  such  a  manner  as  to  maintain  a  connection 

to    the    more   simplified,   but   more   easily   understood, 

development  of  the  previous  section. 

Swift   defines  a  co-ordinate  system  where  the  x-axis 

is   along   the  direction  of  acoustic  vibration  (longitudinal 

in   the  stack)  and  the  y-axis  is  perpendicular  to  the  planes 

of   the  plates  (using  a  parallel  plate  geometry)  with  the  y- 

origin   centered   between   two   plates.   This   gives  a  plate 

spacing   of   2y^  and  a  plate  thickness  of  21,  He  begins  with 

an   equation   of   motion   and   a  continuity  equation  for  the 

fluid: 

^'C  (^  v/Jt )  +  (  V  -v)  v]=-vp+vv^v+(^  +  V/3)7(V-  V)  ,     (  1 1-19) 

and     (nJr/Jt )  + v-(  ?v)=0,  (11-20) 

where  f      is   density,   v   is   velocity,  p  is  pressure,  v  is 

kinematic   viscosity  (Swift  uses  F  for  this  parameter)  and  ^ 

is   second  viscosity.  He  imposes  the  boundary  condition  that 

v=0   at  the  fluid-plate  interface.  He  then  uses  the  equation 

for  heat  flow  in  a  fluid  and  solid: 

fTC<is/=)t  +  v-^s]=  v-(;^VT)  +  (quadrat  ic    velocity    terms),        (11-21) 

and  p  c     ((^T    /c)t)^A   v^T    ,  (11-22) 

)s    s         s  s         s 

where  ;\      is   thermal   conductivity,   s   is  entropy  per  unit 
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mass,  T  is  absolute  temperature  and  the  subscript  s  refers 
to  the  solid.  The  applicable  boundary  conditions  are  that 
T=T  and  A^T/Jy=A  JT  /^y  at  the  fluid-plate  interface. 
Equations  for  each  of  the  needed  parameters  are  then 
linearized,   assuming  that  they  all  oscillate  at  the  angular 


frequency  u^; 


P=P  +p. (x)e 
^   m  ^1 


i  u>t 


^=fm-'?l^^'y^^ 


icut 


v=xu^(x,y)e    +yv^(x,y)e 


T=T  (x)+T,(x,y)e 
m      1    -^ 


iu)t 


T  =T  (x)+T,   (x,y)e 
s   m      l,s    "^ 


lu^t 


and    s=s  (x)+s,(x,y)e    , 
m      1    "^ 


11-23) 
11-24) 
11-25) 
11-26) 
11-27) 
11-28) 


where  the  subscript  m  denotes  a  mean  value  and  the 
subscript  1  refers  to  a  first-order  quantity.  Swift  then 
manipulates  these  equations,  by  integrating  with  respect  to 
y,  into  the  form  of  a  wave  equation  for  p.(x)  and  an 
equation  for  energy  flux  H^  (which  includes  both  heat  and 
work  since  H=Q+W)  along  x.  The  wave  equation  is:  (11-29) 


(J-Df. 


1  + 
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'  m 


^  d   1-f, 


dp 


Pi+ i 
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u^^  dx 


dx 


f  -f 


-P- 


dT   dp. 
ra   ^  1 


=<d 


L   fm    _^ 
and  the  equation  for  energy  flux  is 
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^'  (  l-j ) (  1+  ^  )   dx   dx 
s 
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dp 


2^f 


Imi 
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dx 
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1-^- 


.  T^P(f,  -f.  ) 


(  1+^   )  (  l+(7) 
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^y^^p 


dT^  dp^  dp^ 


Im 


2u.'  P  (  1-iT)   dx   dx   dx 
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fv-^ 


(f«  -f^  )(l4-t^)(f^/fj 


(  1+^   )  (  1+(T) 
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dT 


■irry^A+U^) 


m 


dx 


( 11-30) 


We  have  used  'T='7c  /?.    as  Prandtl  number, 

'    P 

f,  =tanhC(  l  +  i)y^/iK]/n  1  +  i  )y(^/-iH  :i , 
f^  =tanhC  (  1  +  i  )y^/S^]/[  (  1  +  i  )  y^^,/ iv  3, 


(11-31) 
(  11-32) 


and  t  =   c  i.  tanhC  (  1  +  i  )y,-y  S„  ]/ p  c  i  tanh  [  (  1  +  i  )  1  /  i<- ]  ,  (  I  I -33  ) 
smp*^  •^^"'sss  ^ 

where  i  is  a  specific  penetration  depth  denoted  by  the 
appropriate  subscript,  k  for  thermal  penetration  depth  in 
the  fluid,  s  for  thermal  penetration  depth  in  the  solid  and 
V  for  viscous  penetration  depth  in  the  fluid  (which  is 
roughly  the.  distance  that  shear  momentum  diffuses  in  time 
1/w;).  Swift  then  applies  some  simplifying  assumptions  to 
produce   equations   that   can   be   related  to  the  results  of 


section   2.   A   "boundary   layer"  approximation  (i.e.,  y 


0 


and   l"'^^.)   sets   the  hyperbolic  tangents  to  unity.  Since  in 

practical   engines   the   plate   spacing   is   2   to  4  thermal 

penetration   depths   (or  greater)  and  the  hyperbolic  tangent 

function   is   within   a  few  percent  of  unity  at  y^=2t.K,  this 

is     a     reasonable    approximation.    A    "short    stack" 

approximation   (i.e.,   /ix  A)   yields  relations  for  p.  and  u. 

that    are   identical   to   those   of   section   2 ,   with   the 

exception   of   a   (1+1/y^)   term  in  u,  that  accounts  for  the 
^  -^0  1 

velocity  being  higher  inside  the  stack.  This  approximation 
also  gives  us  the  fact  that  p  ,  u  and  all  thermophys i ca 1 
properties  are  independent  of  x  within  the  stack.  Finally, 
the   zero   viscosity   assumption  is  made  (although  viscosity 
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will  be  added  back  in  later)  thereby  setting  <r  and  f^  equal 
to  zero.  Using  these  approximations  Equations  11-29,  11-30 
and  11-33  can  be  reduced  to: 


Pi* 


o  a^   d 
Jm 


u;2   dx 


1 


1  dp 


1 


f    dx 


L  'm 


()f-l)  &« 


(  1  +  i ) (  l+(  )y 


•(r-l)p^  , 


0 


dT 


T  P 
H.^  =  --Trs,  -Ji-p  \,    ^  (  r-  1  )  -Tr(  y  A  +  1  A^  )  -'^    , 
"4    1+^;  '"   dx 


( 11-34) 


( 11-35) 


and 


?  c  u 
'm  p 


( 11-36) 


9    C      i> 
'S     S     S 

where   all   terms   have   previously   been  defined.  It  can  be 

seen   that   the  first  term  in  Equation  11-35  is  identical  to 

the   hydrodynamic  heat  flux  Q^  derived  for  the  single  plate, 

except   for   the   1/(1+^.  )   term.  This  term  accounts  for  the 
^  s 

finite   heat  capacity  of  the  plate  and  modifies  the  boundary 

condition   for  T.  at  the  fluid-plate  interface  to  a  non-zero 

value,   thereby   reducing  the  heat  flux.  For  the  short  stack 

W^   is   small   (W^/Qo=      ^/r  =  iiT/T  r«l)   so   heat  flux  and 
2  2   2   carnot       ra 

energy  flux  are  approximately  equal.  The  second  terra  in 
Equation  11-35  accounts  for  ordinary  conduction  of  heat 
down  the  temperature  gradient  by  the  fluid  and  solid. 

The  next  step  is  to  derive  an  equation  for  the 
acoustic  power.  Since  net  acoustic  power  cannot  enter  or 
leave  the  engine  in  the  y-direction  the  acoustic  power 
generated    or   absorbed   must   involve   the   difference   in 


average   acoustic  intensity  p,u.  between  the  two  ends  of  the 
stack.   Using   this   and  the  fact  that  dp./dx^-i"^^  u.,  Swift 
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manipulates   Equation   11-34   to   arrive   at   a  relation  for 
acoust  ic  power  *• 


1        (y-l)u, 


P^l  Mr-l)  , 


(  11-37) 


f  aM  1+ t  ) 
^m      s 

which   is   the  same  as  Equation  11-13,  with  the  exception  of 

the   l/(l+t  )   term   which   was   discussed  above.  From  these 
s 

results  it  is  clear  that,  under  these  simplifying 
assumptions,  the  short  engine  performs  in  a  manner  similar 
to  that  predicted  by  the  single  plate  results. 

We  now  follow  Swift  in  repeating  this  derivation, 
but  include  non-zero  viscosity.  The  inclusion  of  viscosity 
adds  some  complexity  to  the  situation.  Under  these  circum- 
stances u.  becomes  a  function  of  y  and  the  concept  of 
critical  temperature  gradient  becomes  less  well  defined. 
Since   T.   now   depends  on  y,  in  an  irregular  fashion,  there 

is   no   natural   definition   of  ^T   .^.  However,  we  will  use 

cr  i  t 

(as   did   Swift)   r=7T  /VT       .^      and  VT       , ^=1    pp,^/P  c  <u,">, 

m    cnt         crit   m  ^1   )m  p   1 

where  the  superscript  s  defines  the  value  for  the  standing 
wave.   The   addition   of  viscosity  modifies  the  relation  for 


dp  /dx  to 


dp 


1_ 
dx 


J  in       1    ]  m   1 


( I I-3S) 


1-f, 


1-f, 


which  is  complex  even  though  <u   >  is  real. 

Substituting  into  Equations  11-29  and  11-30  and 
applying  the  boundary  layer  approximation  yields  the 
following  relations: 
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ra  ^  1    1 


(1+t^)  (i  +  o  ( i-iv/y(2,+^(  '-/y0^' 


1+  I.T  '+  ir+  iJ 


1+f 


1  +  cr  ' 
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dT 
dx 


HI 


(  11-39) 


and 


Pi" 


0     A 


1^'   dx 


5     dx 
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(  1  +  i  )(!+':  )y^ 


(  1  +  j  ')  (  l-f .  ) 


-1 


( 11-40) 


For   the  zero  viscosity  case,  cr=0,  these  relations  reduce  to 
Equations  11-34  and  11-35. 

Work   flow   is   calculated,   using   the  short  stack 
approximation,  from  the  wave  equation.  This  gives: 


(»-l)u>(p^^)  i 


—  1 


4       f  aMl+^;  )   (  1  +  cr^)  (  l-5y /y^+iv  ^/2y^' ) 

/m       sL  ^■'0"^0    J 
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Jm   1 


l-S^/y^+5yi/2y^^ 


( I  1-41) 


which  reduces  to  Equation  11-37  when  (r=0. 

From   these   relations   it   is   clear  that  viscosity 

makes   the  situation  much  more  complex.  Since  the  goal  is  an 

expression    for   the   coefficient   of   performance,   it   is 

necessary   to   determine   Q  /W  .  We  know  that  H^=Q_+W  ,  such 

that   C0P=Q„/W^=(H^-W^)/W.-  =  (H„/W^)-1.   In   order  to  maintain 
z   ^   ^     z.        z.        z.     z. 

the  essential  effects  of  viscosity,  while  reducing  the 
complexity  of  the  resulting  equations,  we  will  include 
viscosity  only  to  the  lowest  order  present,  that  is  <x''. 
Additionally  we  will  assume  that  y  -  i^ ,  which  is  the  case 
for   most   practical   engines   as  previously  discussed.  This 
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gives   us  that  iv/y^x-^  ^v  /  *k  =*^ '•  Applying  these  approximations 


to  Equations  11-40  and  11-41  yields 


1    T^pp.^<u^>  dT 

H^=--Tri, -£i — i L__(r-i)-F(y^A  +  lA^)--i!i  , 


(1+^  )(l-j') 
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0' 
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dx 


and 


1     ( y  - 1 )  u.  ( p  ^ )  i     1 
W=-Tr£^iix (  r-i)--Tr6^iixa^^  <u^>'  ,   (11-43) 

^4       e  aM  1+c  )        4        ^ 
>m  s 

Now  we  can  take  the  ratio  of  these  two  relations  and 

get   a   useable,  understandable  relation  for  the  coefficient 

of   performance,   which   is   what   we   are  trying  to  improve 

experimentally.  The  result  is: 

(1  /4  )  TTS K T„  pp  .  ^<  u  .  ^  X  r-  1)  dT„ 

ml    1         TT,    ,  ,  T  ,   ra 

Tr(y^A  +  lA^; 


COP= 


(1+t  )(l-(r^') 
s 


dx 


-1  .  (  11-44) 


(l/4)Tri^iix(^-l)u^(p   )Mr-l)  1 


--TTi^  Ax^^P  <u,  >  ' 
A  'ml 

4 


p  a^ ( 1+  ^  ) 
'  m      s 

The  above  relations  are  approximate  expressions,  which 
should  apply  to  realistic  situations,  with  the  only 
remaining  assumptions  being  1)  the  short  stack 
approximation,  2)  the  boundary  layer  approximation,  and  3) 
a  small  Prandtl  number  approximation.  We  have  previously 
explained  that  the  boundary  layer  approximation  is 
reasonable  for  the  plate  separations  of  our  stacks.  The 
short  stack  approximation  is  also  not  particularly 
restrictive  for  realistic  engines.  With  a  Prandtl  number  of 
approximately  0.027  for  helium-xenon  we  have  a  relatively 
small   Prandtl   number,   so   this  approximation  is  also  very 
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reasonable.  In  his  comparison  to  theory,  Hofler  (1986) 
shows  that  in  relatively  "long"  acoustic  engines  Equations 
11-29  and  11-30  predict  the  performance  within  20  to  30 
percent . 

4.   Design  Considerations 

The  purpose  of  this  section  is  to  discuss  the 
primary  modification  that  we  have  made  based  on  theory, 
that  is,  the  reduction  of  the  Prandtl  number  of  the  system. 
To  evaluate  how  this  affects  performance  we  will  first  look 
at  Equation  11-44.  We  know,  from  experience,  that  the 
second  term  in  the  numerator  is  small  compared  to  the  first 
term.  We  also  know  that  i  is,  in  general,  small  compared 
to  unity,  so  we  will  set  i  =0  for  evaluation  purposes.  The 
product  T  p  is  known  to  equal  unity  for  ideal  gases. 
Applying   these   conditions,   with  Swift's  iv=S^(r*'  and  y^.^^^ 


0 


we  can  reduce  Equation  11-44  to 


( 11-45) 


COP  = 


{C-HTTy^p/<u,^>(  r-1)  ]/(  l-(r^S}-C'^■(y^^+l;\  )dT  /dx] 
■^0-^11  ■'0sm 


s  , 


-1 


{CHfTy^Ax(Jf-l)p^^  (  r-l)]/(^^a^  ) -C  HTTy^j  ''ixu/p^<u  ^^ >  ^ 
It  is  clear  from  this  that  as  <s  is  reduced,  l-a'  is 
increased  resulting  in  an  overall  decrease  in  the  numerator 
of  our  relation,  and  a'  is  reduced  resulting  in  a  decrease 
in  the  denominator  of  our  relation.  These  two  effects,  when 
evaluated  fully,  show  that  COP  will  decrease  with  reduced 
Prandtl  number.  However,  we  know  based  on  our  experiments 
that  this  cannot  be  true  for  our  system.  This  effect  that 
is   apparently   the   opposite  of  what  we  expect,  and  know  to 
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be    true,   is   therefore   due   to   some   (or   one)   of   the 

approximations   that   ue   have   made   being   invalid  for  our 

system.   As   we   are   about   to  show,  the  relations  given  in 

Equations    11-39    and    11-41   yield   the   proper   result, 

indicating   that   the   inclusion   of   viscosity   only  to  the 

lowest   order  does  not  provide  adequately  for  the  effects  of 

Viscosity  in  our  system. 

Looking   at   Equations   11-39   and  11-41  we  see  that 

they   are   fairly   complex.   We  will  use  the  simplifications 

that  K    =0,   T  p=l   and  ?   a'=yp  ;   and  we  will  consider  the 
s       m  '  m     ^m 

thermal   conduction   term  in  the  heat  flow  to  be  negligible. 
We   also  know  that  6^  =  6^/ff",    and  we  will  use  p.  =p^coskx  and 

3 

<u,  >=a  (  1/ JC)  ( p^/p  )sinkx   to  further  simplify  our  relations. 
1  ^0  -^m  sr         J 

This  allows  us  to  re-write  Equations  11-39  and  11-41  as: 


4Ho^P       acoskxsinkx 
z   m 


■^SkP^'  (  l+cr)Cl-i^/y^+*i(  i^/y^)  n 


1+0-  ■'+  cr 


l+(r 


- — ( 1+^  '-■'v/y0) 


(11-46)  and 


4W2YP 


m 


P0''"'^K 


=  (  V  -  1  )u;ixcos  '^k: 


(  1+f  •^)(  i-iv/y^+^(iy/y^)^ 


-1 


0 


0 


>^  ix  ff  s  in  ^kx 


ti-iv/y0+^(*v/y^)^] 


( I  1-47) 


Now  we  can  combine  Equations  11-46  and  11-47   to  get  a 
relation  for  COP  as  follows: 


coskxs  inkx 


COP  = 


kx(  1+0  (  l-i^/y^+^(   ^i;/y,J  ') 


-  —  -1  , 


'^'  ^0 


(  1  1-48) 


where  the  brackets  { }  represent  the  following  relation 


(  JT-  1  )cos  ^kx 


Vi. 


-1 


(T '  s  i  n  '  k  X 


1-6^/y^+^A  ^v/y^)  'j 


_(l+(r  ^)  (l-i^/y^+3t.(  6,/y^)  n 
and   we  have  used  a/u>  ^x''a/wx=  1/kx.  It  is  not  clear  from  this 


relation  how  COP  will  vary  with  Prandtl  number.  However,  if 
we  define  a  set  of  variables  and  perform  a  numerical 
analysis  we  can  see  the  desired  effect.  We  will  define  a 
quantity,  b    =[t^by^)    =constant,  as  the  geometric  mean  of  the 


H  ^ 


penetration  depths,  such  that  6^  =  6  u  and  ^v/y/7)=(  ^vi/y/x ) ''' '^• 
This  roughly  corresponds  to  our  experiment.  Setting 
reasonable  values  for  the  four  parameters  in  our  relation: 
kx=0.2,  r=0.7,  S  /y  =<d,'5  and  V  =  i.667,  we  can  see  how  COP 
varies  with  Prandtl  number.  Figure  II-2  is  a  plot  showing 
how  heat  ( H^ )  ,  work  (W^)  and  COP  vary  with  Prandtl  number, 
using  Equations  11-46,  11-47  and  11-48.  It  is  clear  that  as 
Prandtl  number  is  reduced  COP  increases.  It  should  be  noted 
that  as  Prandtl  number  becomes  very  small  we  actually  see 
COP  decrease  with  reduced  Prandtl  number,  as  we  estimated 
using  the  approximation  with  only  the  lowest  order 
viscosity  terms.  Since  Swift's  engine  (1989)  with  liquid 
sodium  has  a  Prandtl  number  of  approximately  0.004,  his 
approximation  is  in  fact  valid  for  his  system,  although  it 
is  not  for  ours,  as  we  discussed  above. 
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B,   BINARY  GAS  MIXTUURE  THEORY 
1.   Sound  Speed  Theory 

One  of  the  parameters  of  importance,  when  consi- 
dering refrigerator  performance,  is  the  velocity  of  sound 
in  the  fluid  being  used.  The  purpose  of  this  section  is  to 
show  how  this  value  is  calculated  for  binary  gas  mixtures. 

In  general,  the  sound  speed  "a"  in  a  gas  is  given  by 

the   relation:    a=(J'R^T/m)  ,   where  "^(=0  /c   is  the  ratio  of 

0  P   V 

isobaric   to   isochoric   specific  heats,  R^  is  the  universal 

gas   constant,   T   is   the  absolute  temperature  and  m  is  the 

molecular   weight   of   the   gas.   This   relation  holds  for  a 

binary   mixture   of   gases   as  well.  However,  Y and  m  may  be 

different   for   a  mixture  and  must  be  calculated.  From  Huang 

(1976)   we   get  a  relation  for  the  specific  heat  capacity  of 

a    binary    mixture:   (c  )..  =  Cx.M.c  .+x^H,c  ,  ]/Cx , M . +x  .  M , ] , 

P  ij    1  1  Pi   J  J  PJ     1  1   J  J 

+  h 

where   x.  is  the  mole  fraction  of  the  i    species,  H.  is  the 
1  1 

+  Vi 

molecular   weight  of  the  i    species  and  c    is  the  isobaric 

«•  pi 

specific   heat   capacity   of   the   i     species,   A   similar 

relation   holds   for   c    as   well.  The  equivalent  molecular 

v  ^ 

weight  ra,  of  the  mixture  is  given  by:   m=x.H  +x,H,. 

1  1   J  J 

Combining   the   above  relations  we  get  a  final  equa- 
ion  for  the  speed  of  sound  of  a  binary  mixture,  as  follows: 
a={CR^T(x.H.c  .+x,M.c  .  )  D/C  (  x  .  M  . +x  .  H  .  )  (  x  .  M  .  c  .+x.M.c  ,  )  ]  >  *'\ 

0       1   1   pi     J   J   PJ  1   1     J   J       1   1   VI     J   J   VJ 

Since  ^  is  approximately  5/3  for  any  ideal  monatomic 

gas   it  is  clear  that  any  variation  due  to  c   and  c   will  be 
°  -^  p       V 

small.   Therefore,   the   primary   effect   on  sound  speed,  of 
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mixing  two  gases  is  that  it  is  inversely  proportional  to 
the  square  root  of  the  total  mass.  At  300  degrees  Kelvin 
the  sound  speeds  of  the  mixtures  used  in  our  experiment  are 
as  f o  1  lows : 


a -1008  m/s 
a-::614  m/s 
a^460   m/s 


pure  helium 

helium-argon  (19%  Ar ) 

helium-xenon  (12%  Xe ) 
2.   Prandtl  Number  Theory 
a.   General  Discussion 

Before  going  into  the  details  of  how  Prandtl 
number  varies  for  binary  gas  mixtures,  it  is  important  to 
understand  the  general  concept  of  Prandtl  number,  and  to 
discuss  the  expected  variation  based  on  this  concept. 
Prandtl  number  is  a  diraens ion  less  ratio  that  relates 
transport  properties  of  fluids.  In  Fluid  Mechanics  it  is 
given  the  symbol  "Pr"  and  in  Thermodynamics  we  use  "u". 
Prandtl  number  can  be  represented  in  many  different  ways. 
For  calculat ional  purposes  we  define  Prandtl  number  as 
(T=i^c  /A,   where  r^      is  dynamic  viscosity,  c   is  the  isobaric 

c  r 

specific  heat  capacity  and  A  is  the  thermal  conductivity  of 
the  fluid.  This  relation  is  used  because  these  parameters 
are  well  tabulated  and  more  readily  calculated  than  others 
that  are  available.  For  an  understanding  of  the  physics  of 
Prandtl  number  and  how  it  relates  to  the  refrigeration 
process   we   define   Prandtl  number  as  .t=v/k,  where  -v  =  r?/^     is 

the    kinematic    viscosity    and  K  =  Vfc    is   the   thermal 

'     P 
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diffusivity.   This   definition   is   used   because   kinematic 

viscosity   represents   momentum   diffusion  through  the  fluid 

due   to   a  gradient  in  the  velocity,  and  thermal  diffusivity 

represents    the    diffusion    of    thermal   energy   due   to 

temperature   gradients   in  the  fluid.  This  says  that  Prandtl 

number    is    a   ratio   of   momentum   diffusion   to   thermal 

diffusion.   A   third   way   of   describing   Prandtl  number  is 

(r=(  iy /ij^  )  ' ,   where   Sv   is  the  viscous  penetration  depth  and 

&K   is   the   thermal   penetration   depth.   The   capacity  for 

diffusion   of   a   given   quantity  is  directly  related  to  the 

magnitude    of    the   associated   penetration   depth.   Since 

5v  =  (2V/u»)     is   a   function   of   kinematic   viscosity,   and 

&K  =  (2K/uj)    is   a   function   of  thermal  diffusivity,  we  will 

restrict  the  discussion  to  v  and  K  effects. 

Hirschfelder   et   al.,   (1954,   p. 14)   derive   a 

relation,   using   simple   kinetic   theory,   for   the  thermal 

conductivity   of   a   monatomic   gas,  ^=15R^>7/4H,  where  R^  is 

the   universal   gas   constant   and  H  is  the  molecular  weight 

ofthe   gas.   From   this   relation  and  the  knowledge  that  c  - 

c  =R   and  ^=c    /c    ,   where   R=R^/M,   it  is  easily  shown  that 
V  p   V  0  ^ 

iT  =  '^c  /A=4  y/ (  15(  y- 1 )  )  for  a  monatomic  gas.  It  is  also  known 
that  y  =5/3  for  any  monatomic  gas,  which  results  in  a-=2/3 
for  all  monatomic  gases.  Calculations  using  rigorous 
kinetic  theory  also  show  that  Prandtl  number  is  approxi- 
mately 2/3  for  monatomic  gases.  This  theory  will  be 
discussed  in  the  next  section. 


4G 


When  two  monatomic  gases  are  combined  the 
molecular  interactions  produce  changes  in  the  physical 
properties,  resulting  in  Prandtl  number  deviating  from  the 
2/3  value.  For  a  mix  of  a  heavy  gas  and  a  light  gas,  as  we 
have  used,  the  Prandtl  number  drops  rapidly  as  the 
percentage  of  heavy  gas  increases  from  zero.  This  is  due 
primarily  to  a  very  rapid  rise  in  the  molecular  density 
(  =grams/mole)  of  the  gas  as  the  concentration  of  the 
heavier  gas  increases  from  zero  to  some  small  percentage. 
As  an  example,  a  mix  of  two  percent  xenon  with  98%  helium 
produces  a  64%  increase  in  this  molecular  density.  Since 
the  kinematic  viscosity,  V  =  '^/^,  is  inversely  proportional 
to  this  molecular  density,  while  dynamic  viscosity  is 
relatively  invariant,  there  is  a  rapid  drop  in  kinematic 
viscosity  and  therefore  momentum  diffusion.  Thus,  we  see 
that  the  momentum  diffusion  process  is  dominated  by  the 
presence  of  the  heavier  gas  molecules.  The  thermal 
diffusivity,  K=A/fc  ,  on  the  other  hand  is  inversely 
proportional  to  the  product  fc  ,  which  is  essentially 
constant  with  changes  in  concentration,  and  directly 
proportional  to  the  thermal  conductivity,  thereby  changing 
strictly  due  to  variations  in  thermal  conductivity.  Re i f 
(1965)  represents  thermal  conductivity  by  A=c(kT/m)  ',  where 
c  is  a  heat  capacity  per  molecule  (which  is  essentially  a 
constant),  k  is  the  Boltzmann  constant,  T  is  the  absolute 
temperature   and  m  is  the  molecular  weight.  It  is  clear  from 
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this  relation  that  the  thermal  conductivity,  and  therefore 
thermal  diffusivity  and  thermal  energy  diffusion,  will  be 
much  greater  for  a  lighter  gas,  and  their  values  will  be 
dominated  by  the  concentration  of  the  lighter  gas  in  the 
mixture.  Combining  these  two  effects  it  is  clear  that  both 
momentum  and  thermal  energy  diffusion  rates  decrease  with 
increasing  concentration  of  the  heavier  gas.  However,  the 
rate  of  change  in  momentum  diffusion  is  much  more  dramatic 
at  the  lower  concentrations  and  drops  off  rapidly  as  the 
concentration  of  the  heavier  gas  increases,  while  the  rate 
of  change  of  thermal  energy  diffusion  is  much  more  gradual 
and  continuous  over  the  concentration  range.  This  results 
in  the  ratio  of  the  two,  which  is  effectively  the  Prandtl 
number,  dropping  rapidly  from  the  2/3  value  for  the  pure 
light  gas  as  the  heavier  gas  is  initially  introduced, 
reaching  a  minimum  as  the  momentum  diffusion  steadies  out 
and  rising  back  to  the  2/3  value  as  energy  diffusion 
continues  to  decrease  while  the  heavier  gas  concentration 
approaches  100%. 

b.   Kinetic  Theory  and  Quantitative  Analysis 

The  purpose  of  this  section  is  to  discuss  the 
theory,  with  its  assumptions  and  limitations,  and  the 
equations  used  to  calculate  the  Prandtl  number  of  a  binary 
gas  mixture.  Much  of  the  information  presented  here  is 
reproduced  in  Appendix  A  in  order  to  allow  the  removal  and 
use   of   Appendix   A  without  reference  back  to  the  main  body 
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of  the  text.  Hirschf e Ider ,  et  al.  (1954)  introduce  a  basic 
kinetic  theory  in  Chapter  one,  develop  a  rigorous  kinetic 
theory  in  Chapter  seven  and  present  equations  for  the 
calculation  of  thermal  conductivity  and  dynamic  viscosity 
for  binary  gas  mixtures  in  Chapter  eight  of  their  text. 
Using  these  relations  with  another  for  the  isobaric 
specific  heat  capacity  of  the  mixture  we  can  determine 
Prandtl  number,  <j=7c  /^  .  Rather  than  reproduce  the  work  of 
Hirschf e Ider ,  et  al.  (1954)  we  will  present  a  discussion  of 
the  methods  they  used  in  deriving  their  equations, 
including  the  assumptions  and  limitations  that  apply.  Then 
we  will  present  all  of  the  equations  necessary  for  the 
calculation  of  Prandtl  number  for  a  binary  gas  mixture,  and 
a  sample  curve  of  Prandtl  number  as  a  function  of 
concentration  that  was  created  using  the  output  data  of  the 
program  presented  in  Appendix  A,  which  uses  these 
equat  ions. 

The  development  of  the  kinetic  thoery  of  gases 
is  based  on  knowledge  of  the  distribution  function  that 
represents  the  number  of  molecules  of  a  specific  species, 
in  a  unit  volume  element,  about  a  point  in  space,  with 
velocities  in  a  unit  range  about  a  specified  velocity,  at  a 
given  instant  in  time.  This  distribution  function  is 
represented  by  f,(r,v.,t).  H irschf e Ider ,  et  al.  derive  the 
Boltzmann  equation,  in  terms  of  this  distribution  function, 
which   is   the   basis  for  discussion  of  transport  properties 
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of  gases.  The  Boltzmann  equation  is  given  as: 


if 


af 


1  /   if 


.+  V 


iU-  X.-i  =2 


Jt 


9r      m 


^V 


fi'^'-^^  «ij J 


bdbd^dV . (11-49) 
/   J    VI   J    1  jy  ^  1 J 

i  \   '  "  i 
where   the   superscript   "1"  denotes  a  post-collision  value. 


g  is  the  pre-col 1 ision  initial  relative  velocity,  b  is 
the  impact  parameter,  m.  is  the  mass  of  the  i  species,  X 
is  an  external  force  on  the  molecule,  V  is  the  velocity  of 
the  molecule  and  i  is  a  rotational  impact  parameter.  The 
hydrodynamic  equations  of  change  (conservation  of  energy, 
conservation  of  momentum,  etc.)  are  derived  from  the 
Boltzmann  equation. 

Hirschf elder ,  et  al.,  use  a  perturbation  method 
of  successive  approximations,  which  they  attribute  to 
Enskog,  to  obtain  solutions  to  the  Boltzmann  equation. 
Using  this  method  they  derive  a  form  for  the  perturbation 
function  j^,  that  depends  on  space  and  time  through  the 
quantities,    number   density,   mass   average   velocity   and 


(11-50) 


temperature,  and  their  spatial  derivatives; 

,«>=-(  A,  ^lnT/<ir)-(B,^V^/o)r)+nS^  (c/ ^  ^d  J. 
i     i  i   0        J   i     J 

The  coefficients  A  ,  B,  and  C  "^  each  represent  rather 
involved  integral  equations  which  must  be  solved.  They  are 
solved  by  a  series  approximation  method,  using  the  infinite 
Sonine    polynomial    series.   The   Sonine   polynomials   are 


defined  as 
^  (m) 


n 


(x)=sr(-l)J (m+n) !x^ ]/C (n+j ) ! (m-j ) ! j 


( 11-51) 


Linear   combinations  of  a  finite  number  of  these  polynomials 
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are  shown  to  be  the  solutions  to  the  integral  equations  for 
A  ,  B.  and  C.  ^  .  However,  the  transport  coefficients 
(which  is  what  we  are  after)  are  expressed  in  terms  of  the 
Sonine  polynomial  expansion  coefficients.  These  Sonine 
polynomial  expansion  coefficients  involve  complicated 
combinations  of  integrals,  called  "bracket"  integrals. 
Chapman  and  Cowling  (1939)  express  these  bracket  integrals 
as  linear  combinations  of  a  set  of  integrals: 


(1  s  ) 


^ 


■"0  /-<X> 


(e 


(11-52),   where 


0  -'o 


^  )y         (  1-cos  (A))bdbdl{   , 


ij 


is   the   reduced  mass  of  the  colliding 


molecules   i   and   j,  Y  is   the  reduced  initial  relative 

speed   of  the  colliding  molecules,  b  is  the  impact  parameter 

and   X   is   the   angle   of   deflection   of   the   molecules. 

(1  s ) 
Hirschf elder ,   et  al.,  have  tabulated  values  of  these  « 

integrals   for   frequently   encountered  bracket  integrals  in 

Table   I-H   of   their   text.   Hirschf elder ,  et  al.,  go  on  to 

(1  s ) 
define   the   transport  coefficients  in  terms  of  these  ft   ' 

integrals.   It   is   these  relations  that  we  have  used  in  our 
determination  of  Prandtl  numbers  for  binary  gas  mixtures. 

Before  presenting  these  relations  it  is  impor- 
tant to  understand  the  assumptions  used  in  this  derivation 
and  the  limitations  they  impose,  if  any.  Since  the  usual 
definitions  of  transport  coefficients  only  apply  under 
conditions  only  slightly  displaced  from  equilibrium,  our 
results  are  limited  by  this  restriction. 
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Only  binary  collisions  are  considered,  which 
limits  us  to  densities  low  enough  that  three-body 
collisions  are  unimportant.  The  theory  predicts  that 
Viscosity  is  independent  of  pressure  at  constant 
temperature.  Hirschf elder ,  et  al,,  state  that  experiments 
on  nitrogen  have  shown  that  only  a  4%  deviation  in 
Viscosity  occurs  as  pressure  is  varied  from  1  to  60 
atmospheres.  Our  operating  condition  of  10  atmospheres  is 
well  Within  this  range  and  should  therefore  closely  conform 
to  the  theory  in  this  respect. 

The  theory  limits  itself  to  Classical  Mechanics 
and  thereby  precludes  low  temperature  situations  where 
quantum  effects  become  significant.  These  quantum  effects 
are  less  than  1%  for  helium  above  200  Kelvin  (and  even 
smaller  for  heavier  isotopes).  The  coldest  temperature  we 
have  achieved  to  date  has  been  193  Kelvin,  under  a  no  load 
condition,  so  our  normal  operations  are  all  above  200 
Kelvin  and  we  are  not  concerned  with  quantum  effects. 

The  Chapman-Enskog  (1939)  theory  uses  a  first- 
order  approximation  in  their  perturbation  method  solution 
to  the  Boltzmann  equation.  This  is  only  valid  when  the 
gradients  of  the  physical  quantities  (density,  velocity, 
temperature,  etc.)  are  small.  This  means  the  changes  in 
these  properties  over  one  mean  free  path  must  be  small  with 
respect  to  unity.  Based  on  the  mean  free  paths  for  gases  at 
pressures   above  one  atmosphere,  this  condition  is  satisfied 
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except  under  conditions  of  extreme  gradients,  such  as  those 
present  in  shock  waves.  Since  the  observed  pressure 
waveforms  in  our  resonators  are  sinusoidal,  we  do  not  have 
any  shock  waves  and  we  are  not  restricted  by  this 
assumpt  ion . 

The  theory  further  assumes  that  the  dimensions 
of  the  containing  vessel  and  any  obstacles  within  it  are 
large  compared  to  the  mean  free  path,  which  provides  no 
limitations  for  our  gases  and  pressures.  Although  the 
Chapman-Enskog  theory  applies  strictly  to  monatoraic  gases 
the  results  .  have  been  shown  to  be  good  even  for  polyatomic 
gases,  provided  that  the  molecules  are  not  too  non- 
spherical.   Finally,   the   development  of  the  theory  results 

in   a  set  of  equations  for  transport  coefficients  which  rely 

(1  s  ) 
upon    the   previously   discussed   ii   '     integrals.   These 

integrals   are   dependent   upon   the  potential  function  that 

represents    the    molecular   interactions.   Throughout   the 

derivation   by   Hirschf e Ider ,   et   al.,   and   in   all  of  the 

equations   used   for  our  calculations,  the  Lennard-Jones  ( 6- 

12)  potential  is  assumed. 

Now   that   we   have  presented  the  theory  and  its 

limitations   we  will  present  the  equations  necessary  for  the 

calculation   of   dynamic  viscosity,  thermal  conductivity  and 

isobaric   specific   heat   capacity   of   a   binary  mixture  of 

monatomic   gases.   From  these,  as  we  have  already  discussed, 

we  get  Prandtl  number,  where  <t  =  )^c  /7^. 
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is  given  by 


where 


The   dynamic   viscosity  for  a  binary  gas  mixture 


u     1  -11^ 

X„= — i — + — !— i-+— £- 


in,\   cjigii  t^s^i 


,>, 


and 


Zl  =  -*12*k{^^2XjX2 


^'"l^"2''\r''l2^l/1l2^   ^ 


4HjH2 


In  these  equations  x-  and  'x.^  represent  the  mole  fractions 
of   species   1   and   2,   M.   and   H^  represent  the  molecular 

.u^      r         .     1     JO     JA    ♦.(2,2)*,.,(1,1)* 

weights  of  species  1  and  2  and  A.^  =^^  /^^  ,  where 
the  superscript  ♦  refers  to  a  value  that  takes  into  account 
the  deviation  of  a  particular  molecular  model  (in  our  case 
the  Lennard-Jones  (6-12)  potential)  from  the  idealized 
rigid  sphere  model  of  molecular  interactions.  The  values 
for   ^'7i^i  ^^^      ^^^-'i   ^^®  given  by  the  equation  for  a  pure 


substance 


C^]^xl0'^=266.93(MT)  *^/C  (r  ^  ft  ^  ^' ^  ^  *  (  T*  )  , 


where  n  is  dynamic  viscosity  in  gm/cm-sec,  T  is  temperature 
in  *K,  T  =kT/ ^  is  reduced  temperature,  M  is  molecular 
weight,  (T  is  collision  diameter  in  angstroms  and  t/k  is 
the  potential  parameter  in  *K.  The  remaining  term  C^.^]., 
represents   an   artificial   quantity  for  a  hypothetical  pure 
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substance   with  molecules  of  molecular  weight  2h.l\^/  (.h  .+}\^) , 

that   interact   according   to  a  potential  curve  specified  by 

interaction    parameters    <r  and  ^  \2'      These   terms   are 

described   in   terms   of  the  Lennard-Jones  (6-12)  potential, 

where   (t   is  the  distance  of  closest  approach  (in  angstroms) 

of   two   molecules   which  collide  with  zero  initial  relative 

kinetic   energy,   and   t  is  the  maximum  energy  of  attraction 

of  the  two  molecules.  Ci^.^^i  ^^  given  by: 

C^^2^1^^®'^=266.93C2H^M2T/(H^+H2)3^/C(r^2^^'^^'^^*(Tj2*^ 

where   the   subscript   12   on   individual   terms  refers  to  a 

combined   value   for  species  1  and  2.  The  values  for  (r  _  and 

i.^     are   given  by  the  empirical  combining  laws  <r      =¥^{7   +<r    ) 

and    S:?"^  S  ^9^  '    respectively,   and   T^^    is   given   by 


Tj2  -kT/«j2 


is  given  by 
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The   thermal   conductivity  for  a  binary  mixture 


1+Z, 


^^ix^r 


in   cal/cm-sec-*K, 


h  *h 


where      X^= + + 


7 


L^^l^    CAj2^1    ^^2^1 


1          (  1  )             12    ( V )         2         (2) 
Y.= U         + L_±U^    + — = — U  , 
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'\h 


'\2\ 


'^2\ 


Z^=Xj^U^  ^^+2XjX2U^^^+X2^U^^^ 


1/12 


( 1)    '     ♦    ^r"     ♦    1  "i  ^"r"2^ 

"         =-^2    —    -^12    ^^M^ 

15    ^       12\   5    ^^       /  H2        2MjH2 
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15  ^^   12  \  5  ^^    /Hj    2HjH2 


15 


(M  +H  ) 2 

♦  /   12 


12 


''-'^'    *^i!B,,v« 


4H^H2  /  '^^i^i'^^2-'l  ^^^  ^ 


12 


32A.^*\  5  ^^   /   M.H, 
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( 2  2 )  ♦     *   ' 
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„        (T(H.+H^)/2H.M^) 


(T    2  ft 

12   12 


(T   '") 
^M2 


and 


Bj2*=C5.)<l-2>»-4l>"-3'»]/D<'-"V 


The  definitions  for  all  terms  are  similar  to  those  given 
for  the  dynamic  viscosity  relations. 

The  remaining  quantity  necessary  for  the  calcu- 
lation of  Prandtl  number  is  the  isobaric  specific  heat 
capacity   of   the   binary   mixture.   Huang   (1976)   gives   a 

relation   for  this  quantity:   c   ,  =Sm.c  ,/m=Sx,H,c  VSx.H,, 

^  -^  pmix  ipi  iipi  ii 

where   x    is   the   mole   fraction   of   species  i,  M   is  the 

molecular   weight   of   species   i   and   c     is  the  isobaric 

pi 

specific  heat  capacity  of  species  i.  The  required  units  of 
this  quantity,  for  compatibility  with  the  previous 
equations  are  cal/gm-*K. 

All  of  these  equations  were  formed  into  the 
Fortran   computer   program   presented  in  Appendix  A.  Outputs 
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have  been  created  (from  this  program)  for  temperatures  of 
2(2)0"K  and  300'-'K,  and  for  mixtures  of  helium-xenon  and 
helium-argon,  for  all  concentrations  from  100%  helium  to  0% 
helium.  A  composite  plot  showing  all  of  these  cases  is 
provided  as  Figure  II-3,  It  is  clear  from  these  curves  that 
the  Prandtl  number  variation  with  concentration  is  just  as 
predicted  by  the  qualitative  analysis,  given  in  the  first 
part  of  this  section. 

3,   Selection  of  a  Gas  Mixture 

Now  that  ue  have  shown  how  the  sound  speeds  and 
Prandtl  numbers  of  binary  gas  mixtures  vary  with  concentra- 
tion, it  is  necessary  to  consider  how  each  of  these 
parameters  affect  performance.  We  know  from  the 
thermoacoust ic  theory,  presented  in  section  A  of  this 
Chapter,  that  power  density  will  decrease  with  reduced 
sound  speed  (due  to  the  increased  mass),  and  that  COP  will 
increase  with  reduced  Prandtl  number.  However,  as  Prandtl 
number  decreases,  due  to  the  addition  of  a  heavy  gas,  the 
overall  mass  of  the  gas  increases,  and  sound  speed  is 
reduced.  It  is  therefore  necessary  that  the  mixture  have  as 
low  a  Prandtl  number  as  possible  while  the  mass  increase  is 
kept  to  a  minimum  (and  therefore  the  reduction  in  sound 
speed  is  minimized).  Looking  at  Figure  II-3  it  is  clear 
that  we  need  a  mixture  that  is  at  a  point  on  the  very 
steeply  sloped  section  of  the  curve,  as  the  concentration 
begins   to   deviate  from  pure  helium.  Our  selection  of  12,5% 
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Figure  II-3.  Prandtl  Number  vs.  Helium  Fraction 
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xenon  and  18.5%  argon  in  our  helium-xenon  and  helium-argon 
mixes,  respectively,  was  made  based  on  this  requirement. 
With  the  xenon  we  achieved  approximately  90%  of  the  maximum 
possible  Prandtl  number  reduction,  while  getting  only  40% 
of  the  possible  Mass  increase.  With  argon  we  once  again 
achieved  about  90%  of  the  maximum  possible  Prandtl  number 
reduction  but  got  approximately  60%  of  the  possible  mass 
increase.  The  results  of  these  selections  will  be  discussed 
in  Chapter  V. 

C.   STACK  DESIGN  THEORY 

The  purpose  of  this  section  is  to  explain  the  diff- 
erences between  the  two  stacks  built  for  this  project  and 
that  used  in  previous  experiments,  including  the  basis  for 
any  and  all  modifications.  To  provide  an  orderly  flow  of 
information,  the  differences  between  the  stacks  will  be 
fully  described  and  then  an  explanation  of  each  of  the  new 
features  will  be  given.  Details  of  stack  construction  are 
provided  in  Chapter  III. 

We  will  discuss  three  stacks  in  this  section:  the 
original  stack,  built  by  Dr.  Tom  Hofler,  will  be  referred 
to  as  the  "old"  stack,  the  first  stack  built  for  this 
experiment  will  be  referred  to  as  the  "new"  stack,  and  the 
second  stack  built  for  this  experiment  will  be  referred  to 
as  the  "dual"  stack.  The  actual  designs  for  all  of  the 
stacks   are   based   on  the  work  and  ideas  of  Dr.  Hofler.  All 
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of  the  stacks  have  the  same  outer  diameter  of  1,506  inches 
and  the  same  size  center  post  with  a  0.25  inch  diameter. 
Additionally  all  three  have  the  same  plate  thickness  of 
0.003  inches,  and  were  built  using  a  loom,  discussed  in 
Chapter  III,  with  a  line  spacing  of  0.2  inches.  The  primary 
differences  between  the  stacks  are  in  the  areas  of  average 
plate  spacing,  overall  stack  length,  and  final  arrangement 
of  the  fishing  line  used  to  provide  the  plate  spacing. 

The  old  stack  has  an  overall  length  of  3.089  inches  with 
an  average  plate  spacing  of  0.0140  inches,  with  no 
modifications  to  the  fishing  line  arrangement.  The  new 
stack  has  an  overall  length  of  3.066  inches  with  an  average 
plate  spacing  of  0.0133  inches.  Additionally,  it  has  a 
special  pattern  cut  into  the  fishing  line,  as  shown  in 
Figure  II-4.  The  dual  stack  has  an  overall  length  of  3.068 
inches  with  an  average  plate  spacing  of  0,016  inches  at  the 
hot  end  and  0,035  inches  at  the  cold  end.  There  is  no 
pattern  cut  into  the  fishing  line  of  the  dual  stack, 
however,  there  is  a  gap  of  approximately  0.75  inches 
between  the  rows  of  thin  lines  and  thick  lines,  as  shown  in 
Figure  II-5.  Additionally,  on  the  layer  of  plastic  with 
both  sizes  of  fishing  line,  every  other  piece  of  the  thick 
fishing  line  has  been  removed,  starting  at  a  point  18 
inches  from  the  end  of  the  plastic  sheet  that  is  at  the 
center  of  the  roll.  The  only  remaining  difference  between 
the   stacks  is  that  the  old  and  new  stacks  have  center  posts 
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of  a  cloth  phenolic  type  rod  and  the  dual  stack  has  a  nylon 
plastic  type  rod.  This  last  difference  is  strictly  due  to 
availability  of  material  at  the  time  of  construction  and  is 
not  considered  to  be  of  any  significance  in  relation  to 
refrigerator  or  stack  performance. 

Each  of  the  modifications,  relative  to  the  old  stack, 
Will  now  be  explained,  beginning  with  those  for  the  new 
stack.  The  changes  in  this  stack  are  primarily  intuitive  in 
nature,  based  on  experiments  by  Hofler  and  Wheatley,  and 
therefore  do  not  lend  themselves  to  quantitative  analysis. 
The  change  in  overall  length  was  selected  to  provide  a 
space  between  the  end  of  the  stack  and  primarily  the  cold 
heat  exchanger.  This  space  is  desired  to  preclude  any  flow 
restrictions  that  could  occur  with  the  stack  plates  resting 
against  the  plates  of  the  heat  exchanger  surface.  A  piece 
of  fishing  line  (approximately  0,020  inches  in  diameter) 
was  glued  to  the  end  of  the  stack,  to  fill  the  void  left  by 
the  reduced  length,  thereby  preventing  the  stack  from 
moving  around  in  the  resonator  relative  to  the  heat 
exchangers.  The  pattern  of  Figure  I  1-4  was  selected  for  two 
reasons.  Since  the  amount  of  heat  transport  is  related  to 
the  total  surface  area  of  the  plates,  the  removal  of  some 
of  the  fishing  line  would  increase  this  surface  area  and 
thereby  increase  the  heat  transport.  Also,  since  it  is 
virtually  impossible  to  have  precisely  uniform  spacing  of 
the   fishing   line   and  plates,  there  are  slight  differences 
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in  the  longitudinal  flow  channels  in  the  stack.  This  in 
turn  can  produce  flow  variations  and  pressure  differences 
from  one  channel  to  the  next.  Due  to  there  being  no 
continuous  pieces  of  fishing  line  along  the  full  length  of 
the  stack,  with  this  pattern,  the  channels  are  all 
interconnected  and  the  pressure  and  flow  differences,  if 
any,  will  be  minimized.  The  overall  result  of  this  pattern 
was  a  reduction  of  approximately  50%  in  the  amount  of 
fishing  line  present  in  the  stack.  This  reduction  should 
also  have  the  effect  of  reducing  the  overall  obstruction  to 
flow  in  the  resonator.  It  should  also  be  noted  that  this 
pattern  may  have  had  a  significant  effect  on  the  degree  of 
crush  experienced  during  the  rolling  process,  resulting  in 
the  smallest  average  plate  spacing  of  any  of  the  stacks. 
This  plate  spacing  is  very  important  in  the  thermoacoust ic 
performance  of  the  refrigerator. 

The  changes  in  the  dual  stack  are  somewhat  more  complex 
and  theoretical  in  nature.  Within  the  accuaracy  of 
construction  methods,  the  overall  length  of  the  dual  stack 
is  the  same  as  that  of  the  new  stack,  for  reasons  explained 
above.  The  removal  of  every  other  piece  of  the  thick 
fishing  line,  after  the  first  18  inches,  is  based  on 
reducing  the  overall  amount  of  fishing  line  in  the  stack, 
as  discussed  for  the  new  stack.  All  lines  were  left  intact 
for  the  first  18  inches  to  provide  a  more  uniform  plate 
spacing   toward   the   center  of  the  roll,  where  the  diameter 
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of   the   roll   is  not  very  large  with  respect  to  the  fishing 

line   spacing.   The   0.75   inch  gap  between  the  rows  of  thin 

and   thick   strings   was   used   for   two   reasons.  A  certain 

minimal   gap   is   necessary   to   prevent   the   inner  plastic 

layer,   on   the   hot   end,  from  butting  up  against  the  thick 

strings   when  the  ends  of  the  stack  are  pushed  together.  The 

extra   distance,   out   to   0,75   inches,   was   used   to  give 

approximately   the   same  total  length  of  string  as  we  had  in 

the  new  stack. 

Understanding   the   reasons   for   the  dual  nature  of  the 

stack   requires   some  knowledge  of  the  basic  characteristics 

of   large   temperature   span   refrigerators,   as  well  as  the 

basic   heat   transfer   process  that  occurs  in  the  stack.  The 

heat   flow   at   the   two   ends  of  the  stack  is  not  the  same, 

generally   being   significantly  higher  at  the  hot  end,  where 

heat   is  rejected.  This  difference  is  functionally  dependent 

on   efficiency   and   increases   with   reduced  efficiency.  To 

properly   explain  this  effect  some  terms  must  be  introduced. 

Although   the   following   relations   are  simplified  they  are 

adequate   for   the   explanation   of   this   effect.   The  heat 

transfer  capacity  q  is  given  by: 

q(x)=Q(x)/(r(x)-l)  =-Tr6^  p(x)u(x),       (II  -52  ) 

and  the  local  heat  flow  Q  is  given  by: 

Q(x)=-TT&^p(x)u(x)  (  r(x)-l),         (  11-53) 

where   r(x)=t7T/oT   ,^,   vT   ^^   is   the   maximum  temperature 

cr  it      cr  it 

gradient   in  the  x  direction,  IT  is  the  surface  area  per  unit 
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length  of  the  stack,  p(x)  is  the  acoustic  pressure 
distribution,  u(x)  is  the  acoustic  velocity  distribution, 
Sk  is  the  thermal  penetration  depth  (which  according  to 
Swift  (1989),  is  the  distance  that  heat  can  diffuse  through 
the  fluid  during  the  time  l/w),  and  x  is  the  longitudinal 
distance  along  the  stack.  In  a  standard  stack  u(x)  and  r(x)- 
1  are  the  only  parameters  that  vary  significantly  in  the 
longitudinal  direction.  The  quantity  r(x)-l  represents  the 
temperature  gradient  perpendicular  to  the  x  direction  (i.e. 
the  transverse  direction).  This  quantity  will  vary  greatly, 
depending  on  heat  load,  and  must  be  large  in  order  to 
transfer   large  amounts  of  heat.  In  a  large  temperature  span 


refrigerator    r(x)-l 


is   very   small   at  the  cold  end  and 


large  at  the  hot  end.  The  Q(x)  distribution  is  dominated  by 
r(x)-l  .  Conversely,  u(x)  is  generally  small  at  the  hot 
end  of  the  stack  and  large  at  the  cold  end.  It  is  u(x)  that 
dominates  the  q(x)  distribution.  These  effects  give  a  low 
heat  flow  at  the  cold  end,  where  the  heat  transfer  capacity 
is  highest,  and  a  high  heat  flow  at  the  hot  end,  where  the 
heat  transfer  capacity  is  lowest. 

The  object  of  the  dual  stack  design  is  to  reduce  this 
conflict  by  changing  the  distribution  of  q(x),  making  it 
more  uniform  and  similar  to  that  of  Q(x).  Since  r(x)-l  must 
be  geater  than  zero  to  transfer  heat,  and  is  very  dependent 
on  heat  load,  it  was  determined  that  q(x)  was  the  quantity 
that   should   be   modified.  To  accomplish  this,  the  acoustic 


6G 


velocity  distibution  had  to  be  countered  by  some  other 
effect.  The  only  parameter  affected  enough  by  stack  design 
to  accomplish  this  is  TT  By  varying  the  amount  of  surface 
area  along  the  stack  length  the  q(x)  distribution  can  be 
altered.  By  having  half  the  number  of  plates  at  the  cold 
end,  we  cut  the  surface  area  by  a  factor  of  two, 
essentially  countering  the  u(x)  distribution,  and  making 
the  capacity  for  heat  transfer  correspond  more  closely  to 
the  actual  heat  flow  conditions.  Ideally,  the  surface  area 
should  vary  continuously  in  the  longitudinal  direction  to 
precisely  counter  the  velocity  distribution  and  produce  an 
exactly  uniform  q(x)  distribution,  thereby  giving  optimum 
results.  This  situation  would  be  much  more  difficult  to 
construct  and  methods  that  have  been  discussed  to  date  have 
other  drawbacks  that  make  them  less  acceptable  than  the 
current  design.  Although  this  is  the  major  reason  for  an 
expected  increase  in  the  COP  with  the  dual  stack  design 
there  are  two  other  minor  effects  which  should  also  help. 
An  area  not  covered  by  the  simplified  equations  above  is 
that  of  viscous  friction  as  the  gas  moves  along  the  plates. 
The  energy  loss  due  to  this  friction  is  proportional  to 
both  the  square  of  the  acoustic  velocity  and  the  total 
surface  area  exposed  to  the  gas.  Since  the  surface  area  has 
been  cut  in  half  at  the  cold  end,  where  the  velocity  is 
highest,  this  source  of  energy  loss  should  be  significantly 
reduced,   thereby   increasing   the   COP.   The  second  area  of 
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interest  relates  to  the  uniformity  of  the  q(x) 
distribution.  Since  q(x)  depends  on  u(x),  anything  that 
makes  u(x)  more  uniform  will  also  make  q(x)  more  uniform, 
which  in  turn  tends  to  improve  the  COP.  The  velocity 
distribution  will  depend  somewhat  on  the  degree  of  blockage 
in  the  flow  path.  The  plates  have  an  "effective"  size  that 
is  larger  than  their  thickness,  due  to  viscous  effects.  By 
reducing  the  number  of  plates  the  degree  of  blockage  is 
reduced  and  the  velocity  distibution  is  more  uniform, 
thereby  tending  to  increase  the  COP. 
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III.   STACK  CONSTRUCTION 

The  thermoacoust ic  stack,  necessary  for  the  transport  of 
heat  as  discussed  in  chapter  II,  consists  of  a  continous 
roll  of  very  thin  plastic,  the  layers  of  which  are 
separated  by  spacers  attached  to  the  plastic  surface.  This 
geometry  was  selected  instead  of  a  parallel  plate 
arrangement  based  on  both  ease  of  construction  and 
structural  rigidity  of  the  finished  product.  When  the 
dimensions  of  a  stack  are  considered,  plate  spacing  of 
0,015  inches  (or  0.033  inches  depending  on  which  stack), 
plate  thickness  of  0.003  inches,  overall  diameter  of  1.506 
inches,  and  overall  length  of  about  3.066  inches,  it  is 
easily  seen  that  the  parallel  plate  geometry  would  be  much 
more  difficult  to  construct  and  extremely  fragile.  The 
plastic  used  in  this  experiment  was  a  polyester  Mylar  TH 
type.  There  was  no  particular  reason  for  this  selection 
other  than  ease  of  handling.  The  plate  spacers  consist  of 
strips  of  high  grade  monofilament  nylon  fishing  line  which 
are  affixed  to  the  plastic  surface  using  a  spray  adhesive. 

Prior  to  construction,  decisions  must  be  made  concerning 
overall  length  of  the  stack,  which  determines  the  width  of 
the  plastic  sheet,  required  diameter,  which  determines  the 
overall  length  of  the  plastic  sheet  and,  approximate  plate 
spacing,    which   determines   the   size   of   spacers   to   be 
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attached  to  the  plastic  surface.  For  this  project  two 
stacks  were  built.  The  first  had  an  overall  length  of  3.066 
inches,  with  a  diameter  of  1.506  inches  and  0.015  inch 
spacers.  The  second  was  a  dual  stack  with  the  same 
diameter,  an  overall  length  of  3.068  inches  and  spacers  of 
0.015  inches,  with  two  layers  of  plastic  along  half  of  its 
length,  and  spacers  of  0.033  inches  along  the  other  half. 
This  provided  a  variable  plate  spacing  for  reasons 
discussed  in  chapter  II.  The  overall  length  of  plastic 
sheet  needed  to  make  a  stack  can  be  determined  by 
mathematical  analysis  of  the  spiral  wound  roll  once  the 
dimensions  are  specified.  The  resulting  equation  for  length 
is  L=( T/d) (r . ^-r^M ,  where  d  is  the  sum  of  the  spacer  size 
and  plate  thickness,  r.  is  the  radius  of  the  stack  (and 
also  the  radius  of  the  resonator  bore),  and  r^  is  the 
radius  of  the  center  post,  on  which  the  plastic  is  rolled. 
The  center  posts  that  have  been  used  are  0.25  inch  diameter 
and  are  either  cloth  phenolic  or  a  nylon  type  plastic  rods. 
Usually  a  safety  margin  of  twelve  to  eighteen  inches  is 
added  to  the  length  calculated  by  the  above  equation.  This 
is  due  to  the  fact  that  a  precise  fit  is  desired  in  the 
resonator  bore  and  it  is  easier  to  trim  small  pieces  from 
an  oversized  roll  than  to  add  pieces  to  one  that  is 
undersized.  The  width  of  the  plastic  sheet  is  cut  very 
precisely  using  a  sharp  paper  cutter,  with  a  guide  that  has 
been   carefully  positioned  to  the  desired  width.  This  method 
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has  provided  results  with  an  accuracy  of  0,002  inches.  The 
overall  length  of  plastic  sheet  used  for  the  two  stacks  in 
this  experiment  was  120  inches  for  the  first  stack  and  two 
pieces  of  60  inches  each  for  the  dual  stack. 

Once  the  plastic  sheet  has  been  cut  to  size  and  fishing 
line  of  the  required  size  has  been  acquired,  construction 
can  begin.  The  spray  adhesive  used  in  this  experiment  was 
3H  brand  Super  77  spray  adhesive.  Since  attaching 
individual  pieces  of  fishing  line  would  be  very  time 
consuming  and  result  in  a  less  than  uniform  spacing  an 
alignment  device  is  used  to  assist  in  this  process.  An 
aluminum  loom  has  been  constructed  to  provide  uniform 
parallel  spacer  positioning  on  the  plastic  surface.  A  long 
piece  of  fishing  line  is  wound  around  the  hooks  of  the 
loom,  positioned  in  grooves  to  provide  paral lei  1  ism,  and 
fixed  at  both  ends  with  screws  to  hold  it  in  place.  The 
loom  has  an  open  front  and  back  to  provide  access  to  both 
sides  of  the  fishing  line.  Using  a  spray  shield  to  protect 
the  loom,  spray  adhesive  is  applied  to  the  front  side  of 
the  fishing  line.  The  loom  is  then  turned  over  and  aligned 
over  the  plastic  to  provide  parallel,  uniform  positioning 
of  the  strips  of  fishing  line.  Once  positioned  the  loom  is 
set  down  and  a  smooth  aluminum  block  is  used  to  press  the 
fishing  line  onto  the  surface.  The  working  surface  and 
press  used  in  this  experiment  were  smoothly  polished 
aluminum    blocks,    which    were   necessary   to   provide   a 
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relatively  uniform  compression  of  the  fishing  line  when 
pressing.  A  teflon  sheet  was  placed  under  the  plastic  and 
over  the  working  surface  to  prevent  the  strings,  which 
extend  beyond  the  plastic,  from  sticking  to  the  working 
surface.  A  sketch  of  the  loom,  with  fishing  line  attached 
is  shown  positioned  over  the  plastic  sheet  in  Figure  III-l. 
The  distance  between  strips  of  fishing  line  on  the  loom 
used  in  this  experiment  was  0.20  inches.  This  spacing  was 
selected  based  on  experiments  conducted  by  Dr.  Tom  Hofler 
at  the  Naval  Postgraduate  School.  The  loom  covers  a  length 
of  Six  inches  and  therefore  requires  twenty  repetitions  of 
this  procedure  to  complete  a  stack.  The  dual  stack  required 
more  repititions  due  to  there  being  two  different  string 
sizes  on  one  plastic  sheeet  and  a  second  layer  of  plastic 
of  half  the  width,  as  shown  in  Figure  III-2.  Once  the 
fishing  line  is  pressed  down  and  allowed  a  minimal  drying 
time  of  one  to  two  minutes  a  scalpel  or  razor  blade  is  used 
to  cut  the  fishing  line  off  at  the  edge  of  the  plastic.  It 
should  be  noted  that  if  the  line  is  wound  on  the  loom  with 
too  much  tension  the  plastic  will  tend  to  curl  when  the 
line  is  cut.  This  is  obviously  an  undesireable  situation  so 
care  must  be  taken  when  winding  the  line  on  the  loom  to 
provide  Just  enough  tension  to  maintain  the  strips  in  the 
grooves.  After  all  of  the  fishing  line  is  attached  to  the 
plastic  some  trimming  must  be  done  prior  to  rolling.  The 
fishing    line   is   carefully   cut   back   approximately   one 
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Figure  III-l.  Stack  Construction  Illustration 
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sixteenth  of  an  inch  from  each  edge,  to  prevent  potential 
flow  problems  at  the  heat  exchanger  interfaces. 
Additionally,  depending  on  the  stack  design,  a  pattern  is 
cut  in  the  strings  or  some  strings  are  removed,  for  reasons 
explained  in  the  stack  theory  of  chapter  II.  Great  care 
should  be  exercised  in  the  trimming  of  any  patterns  to 
avoid  making  any  large  cuts  in  the  plastic  surface  as  this 
could  produce  significant  flow  variations  during  operation 
of  the  unit. 

After  completing  this  process  the  stack  is  ready  to 
roll.  The  center  post  is  carefully  taped,  using  scotch 
tape,  to  the  end  of  the  plastic  sheet.  The  plastic  is 
carefully  rolled  up,  using  just  enough  tension  to  maintain 
the  shape  and  provide  uniform  layer  spacing,  without 
getting  too  much  crush,  which  would  result  in  a  layer 
separation  much  smaller  than  the  string  diameter.  This  is 
in  general  the  most  difficult  part  of  the  stack 
construction  process.  The  uniform  layer  spacing  is  checked 
visually  by  holding  the  stack  up  to  a  light  or  using  some 
sort  of  magnifying  device.  The  spacing  can  also  be  checked 
quantitatively  baaed  on  the  known  innner  post  diameter, 
the  measured  outer  diameter  and  a  count  of  the  number  of 
layers  along  several  radii.  This  will  provide  an  average 
layer  spacing  which  will  in  general  be  somewhat  smaller 
than  the  sum  of  the  fishing  line  diameter  and  the  plastic 
thickness.   This   is   due   to   a   certain  degree  of  crush  or 
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compression  during  the  pressing  and  rolling  processes. 
After  rolling,  even  with  a  great  deal  of  care,  there  will 
be  some  degree  of  spreading,  or  misalignment,  at  the  ends 
of  the  cylinder.  Using  two  smooth  surfaces  the  ends  of  the 
stack  can  be  pressed  and  the  layers  will  slide  into 
position  to  give  flat  ends  and  an  overall  length  as  precise 
as  the  width  of  the  plastic  sheet. 

Drawings  of  the  two  ends  of  the  dual  stack  are  shown  in 
Figures  III-3  and  III-4.  The  other  stack  has  the  same 
appearance  as  the  dense  or  hot  end  of  Figure  III-3. 
Polaroid  photographs  of  the  actual  dual  stack,  showing  a 
radius  from  the  center  post  out,  at  a  6.4x  magnification, 
are  provided  as  Figures  III-5  and  III-6.  The  finished  roll 
is  then  trimmed  down  in  diameter,  by  cutting  small  pieces 
off  the  length  of  the  plastic  strip,  to  fit  snugly  in  the 
resonator  bore.  The  final  product  is  taped  along  the  seam, 
with  scotch  tape,  to  hold  the  roll  in  the  cylindrical 
shape. 

Although  some  of  the  differences  between  the  two  stacks 
that  were  built  for  this  experiment  are  mentioned  in  here, 
the  details  of  and  reasons  for  the  differences  are 
presented  in  the  stack  theory  of  chapter  II. 
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Figure  III-5.  Photograph  of  Dual  Stack,  Cold 
End,  at  6.4x  Magnification 
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Figure  1 1 1-6.  Photograph  of  Dual  Stack,  Hot 
End,  at  6.4x  Magnification 
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IV.   EXPERIMENTAL  APPARATUS  AND  PROCEDURES 

In  this  Chapter  we  will  present  a  description  of  our 
experimental  apparatus,  including  the  electronics  used,  and 
a  discussion  of  our  experimental  procedure.  We  have  used 
the  prototype  thermoacoust ic  refrigerator,  and  associated 
electronics  package,  that  was  built  by  Hofler  (1986,1988). 
We  have  paraphrased  sections  of  his  discussion  of  the 
system  as  necessary  to  provide  a  basic  understanding  of  the 
apparatus  and  its  operation. 

A.   DESCRIPTION  OF  APPARATUS 

In  order  to  properly  describe  the  apparatus  it  is 
necessary  to  break  it  down  into  two  major  areas;  that  of 
the  high  intensity  driver,  with  its  instrumentation  for 
acoustic  power  measurements  and  that  of  the  refrigeration 
components  (resonator,  stack,  heat  exchanger,  etc.). 

1 .   The  Driver  and  its  Instrumentation 

The  driver  consists  of  a  modified  commercial 
loudspeaker,  shown  in  Figure  IV-1.  The  modifications 
consist  of  removal  of  the  original  plastic  housing  and 
fabric  dome  and  attachment  of  a  thin-wall  aluminum  cone- 
shaped  piston  and  a  surround.  Figure  IV-1  shows  the  unit 
with  a  voice  coil  and  former  1,  iron  pole  pieces  2,  magnet 
3,   aluminum   piston  4,  surround  5,  aluminum  pressure  vessel 
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Figure  IV-1.  Schematic  of  Driver  Apparatus  Showing 
Acoustic  Motor,  Transducers,  and  a  Portion  of  a 
Resonator.  Numbered  Parts  are  Identified  in  the  Text 
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6  and  7,  aluminum  plate  8,  cooling  water  tubes  9,  annular 
plate  (where  pressure  transducer  is  mounted)  10,  dynamic 
pressure  transducer  11  and  a  post  for  the  mounting  of  an 
acce lerometer  12.  The  unit  has  a  two  inch  diameter  voice 
coil  and  an  effective  piston  driving  area  of  0.93  square 
inches.  Additionally,  it  has  a  high  power  capability  that 
is  made  possible  by  suspending  ferrofluid  in  the  magnet 
gap,  to  conduct  heat  from  the  voice  coil,  A  capillary 
filled  hole  is  located  adjacent  to  the  surround  to  break 
the  gas  seal  provided  by  the  surround.  The  flow  impedance 
of  the  capillary  was  chosen  to  be  as  low  as  possible 
Without  affecting  the  acoustic  power  measurements.  The 
assembly  is  housed  in  the  aluminum  pressure  vessel  6  and  7, 
and  the  resonator  section  is  bolted  to  6.  The  heat  removed 
by  the  water  cooling  tubes  9,  is  conducted  through  the 
aluminum  housing  6  and  plate  8. 

The  dynamic  pressure  transducer  11,  is  physically 
attached  to  the  structure  of  the  driver  and  can  therefore 
be  calibrated  along  with  other  aspects  of  the  driver, 
independent  of  the  resonator.  The  transducer  itself  was 
designed  and  built  by  Dr.  Hofler.  It  consists  of  a  thin  y- 
cut  quartz  crystal  disk  with  evaporated  electrodes  and  a 
miniature  FET  impedance  converter  circuit  that  buffers  the 
output  signal.  The  dynamic  pressure  transducer  provides  a 
low  noise,  low  distortion,  phase  accurate  signal  which  is 
calibrated   to   determine   the   volume   velocity   in  a  small 
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calibration  cavity  that  can  "be  attached  to  the  driver. 
This,  in  conjunction  with  the  integrated  acce lerometer 
signal  (which  is  then  a  velocity  signal)  and  the  phase 
difference  between  the  two,  provides  for  an  accurate 
measurement  of  the  acoustic  power  provided  to  the 
resonator.  Hofler's  calibration  of  the  pressure  transducer 
was  accomplished  by  comparison  to  a  "standard"  transducer 
that  had  previously  been  calibrated  to  an  accuaracy  of  :1%, 
by  Swift,  et.  al.  (1982),  resulting  i.n  a  sensitivity  of 
1.292  (tl%)  volts/bar.  The  accuracy  of  the  acoustic  power 
measurement  is  t3%,  as  shown  by  Hofler  (1988). 

As  previously  noted  in  the  discussion  of  Figure  IV- 
1,  the  acce lerometer  is  mounted  on  the  post  12,  which  is 
attached  to  the  back  side  of  the  piston  4.  The  post  is 
necessary  to  thermally  isolate  the  accelerometer  from  the 
piston  (which  becomes  warm  during  operation)  and  to 
position  the  accelerometer  where  its  cable  has  clearance 
from  the  piston  wall.  The  accelerometer  is  a  shear 
piezoceramic  type.  The  piston's  velocity  signal  (which  is 
the  integrated  output  of  the  accelerometer)  is  calibrated 
with  the  dynamic  pressure  transducer  signal  using  a  small 
volume  of  known  acoustic  impedance,  to  get  an  accurate 
value  for  volume  velocity.  This  calibration  may  change  as 
the  properties  of  the  suround  change  with  use  and  must 
therefore  be  repeated  at  periodic  intervals,  to  maintain 
the   accuracy   of  results.  Repeated  calibrations,  over  a  two 
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year  period,  have  shown  that  the  sensitivity  of  the 
pressure  transducer  has  been  stable  within  about  0.1%, 
while  the  volume  velocity  calibration  value  has  shown  an 
overall  decrease  of  about  7%.  It  should  be  noted  that  the 
volume  velocity  value  dopped  just  over  4%  during  the  first 
four  months  of  driver  use,  and  dropped  only  3%  more  over 
the  next  18  months. 

2.   Refrigeration  Components 

The  main  components  of  the  refrigeration  section  are 
shown  in  Figure  IV-2.  The  resonator  (18,  20  and  21  in  the 
figure)  is  constructed  of  several  parts.  Where  temperature 
uniformity  is  desired  heavy-gauge  copper  is  used.  In 
section  18,  containing  the  stack  17,  there  is  a  large 
temperature  gradient  and  an  epoxy-f iberglass  substance  is 
used  to  provide  high  strength  with  very  low  thermal 
conductance.  This  epoxy-f iberglass  section,  by  itself  would 
be  very  permeable  to  gas  diffusion  if  left  alone.  To 
correct  this  situation,  Dr.  Hofler  evaporated  a  thin  film 
of  indium  over  the  outside  of  the  fiberglass  section,  and 
then  wrapped  more  fiberglass  over  the  metal  film. 

The  primary  refrigeration  component  is  the  stack  (17 
in  the  Figure  IV-2)  which  is  constructed  of  a  long  sheet  of 
plastic  that  is  spiral  wound  on  a  0.25  inch  diameter 
plastic  or  cloth  phenolic  rod  16.  The  layers  of  plastic  are 
separated  by  monofilament  fishing  line  spacers.  The  details 
of  stack  construction  were  presented  in  Chapter  III. 
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Figure  IV-2.  The  Primary  Refrigerator  Components, 
Shown  Attached  to  the  Driver  Assembly.  The  Numbered 
Components  are  Discussed  in  the  Text. 


85 


The  other  components  of  importance  in  the 
refrigeration  section  are  the  hot  and  cold  heat  exchangers 
15  and  22  in  Figuure  IV-2,  respectively.  The  hot  heat 
exchanger  provides  for  removal  of  heat  from  the  unit,  to 
prevent  overall  heating  of  the  stack  during  operation.  The 
cold  heat  exchanger  provides  good  thermal  contact  with  the 
applied  heat  load,  from  the  heater  coil  19.  Dr.  Hofler 
constructed  the  heat  exchangers  using  electroplating  and  a 
chemically  removable  aluminum  form,  shown  in  Figure  IV-3 
prior  to  removal  of  the  aluminum.  The  resulting  heat 
exchangers  had  parallel  copper  strips  with  the  planes  of 
the  strips  parallel  to  the  tube  axis.  The  details  of  this 
construction  are  provided  in  the  doctoral  dissertation  of 
Hofler  (1986). 

The  remainder  of  the  physical  apparatus  consists  of 
a  vacuum  system,  that  attaches  below  the  driver  housing, 
thereby  enclosing  the  entire  refrigeration  section,  and  a 
section  of  super-insulation  wrapped  around  the  resonator. 
All  of  this  is  to  provide  thermal  insulation,  thereby 
reducing  nuisance  heat  loads.  A  layer  of  super-insulation 
consists  of  a  layer  of  glass  micro-fiber  web  folloued  by  a 
layer  of  aluminized  polyester  film.  The  entire  resonator 
section   is   wrapped   with   10   to   12   of  these  layers.  The 

vacuum   system   uses   an   ion   gauge   to  monitor  the  vacuum, 

-4  -6 

which    varies    between    3x10     and   5x10     torr   during 

operation,  depending  on  operating  conditions  and  gas  types. 
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B.   DESCRIPTION  OF  ELECTRONICS 

The  electronics  portion  of  our  experimental  apparatus 
consists  of  two  primary  sections.  These  are  the  control  and 
protective  section  and  the  data  monitoring  section.  The 
control  and  protective  systems  will  be  discussed  here  and 
the  data  monitoring  aspects  will  be  covered  in  the  next 
section,  where  we  discuss  the  experimental  procedure. 

The  control  and  protective  electronics  system  provides 
feedback  control  of  the  drive  frequency  and  amplitude, 
protection  against  burn-up  of  the  driver  voice  coil,  and 
signal  conditioning  to  increase  the  precision  of  acoustic 
power  measurements.  The  protective  circuit  is  particularly 
important  in  our  situation  because  the  required  operating 
conditions  are  very  close  to  the  maximum  limits  of  the 
driver. 

Figure   IV-4   is  a  block  diagram  of  the  control  circuits 

for    drive    frequency    and    amplitude.   This   figure   is 

reproduced,   with  minor  equipment  updates,  from  the  doctoral 

dissertation   of   Tom   Hofler  (1986).  The  explanation  of  its 

operation,    for    the    most    part,   is   quoted   here   for 

completeness : 

"The  frequency  control  is  based  on  a  lock-in  amplifier 
that  is  used  as  a  phase  null  detector  with  an  analog 
output.  The  velocity  signal  is  fed  into  the  lock-in 
preamp,  and  the  pressure  signal  is  used  as  a  reference 
channel  input.  The  level  of  this  reference  signal  is 
adjusted  to  be  IV  rms,  in  order  to  minimize  lock-in 
phase  errors.  The  Y-channel  output  voltage  of  the  lock-in 
is  proportional  to  the  measured  phase  deviation  and  is 
fed   into   a  simple  integrating  controller.  The  controller 

88 


LJ 

t- 

1 

a; 

c   o 

O-J 

3   D 

.«  a 

(^   a. 

-J    c 

(J      QJ 

A 

v 

CJ 

-4^ 

c  c 

U- 

CJ 

eC 

C_       3 

—    L- 

D     QJ 
Li_  C-3 

C_3 

C-3 

rv 

QJ       Q. 

3     c 

o  <a 

P 

c-  > 

A 

O     1— t 

.       . — 1 

Q_ 

c  s- 

c     o 

a 

OQ 

O       L 

o     E 

-4-3 

CJ      -4-3 

E —    or 

c 

o 

ji:  c 

.-<    0) 
O    Ol 

.-1         QJ 

..^ 

o  — 

CD   c- 

Q-  a   /-\ 

-J    — 1 

C3 

XI   3 

-p  — 1 

[.      o 

■OQ 

Ol 

(D 

ID     O 

Ol     c_ 

C       1_ 

(U  in 

u-f=> 

n  t,  e 
on  t 

-4-3        —1 

D     ' — ' 

-4-3 

c 

.-4          QJ 

3 
o  — 

1—1     CJ 

I-       O 
Ol      (_ 
QJ     -J 

O       Ol 

a_     a 

w  CO 

> 

(UT3 

, 

A 

rs 

-4J 

-4-3         4—4 

T 

kJ 

1 

c     o 

t— 1     CJ 

U3          U 

in    :> 

(0   <U 

t_ 

3  C 

01 

0) 

^" 

A 

in  E 

c^ 

c 

a     a 

o  C 

C     -J 

Z>     C-3 

-  © 

' 

O       3 

C  -M 

Q. 

E 

-C       Ol 

A 

C- 

O  <D 

CJ     -^ 

CJ       OJ 

tC 

1         3 

CD     -4-3 

>-    c^ 

1_ 

O   £- 

c 

O       QJ 

-4-3              D 

(DO, 

c_ 

Q 

1 — 1 

0) 

u 

C  — p 

en 

E 

c 
CJ      o 

U   O 

a 

a:    CJ 

O 

c 

u 

0)    (U 

:> 

s:  (D 

Ol 

o 

C_        C    -4-3 

->                ^ 

1-T) 

QJ 

_J 

QJ    C     3 

3 

-J 

c.-   o   a 

r^ 

S) 

(D-C    C 

'^ 

o  -^ 

CD 

az  c_3 1— I 

1 — 1 

az 
to 

^^ 

c 

(_        .-4 

O       D 

e  Q_ 

a. 

E   -J 
a     3 

I 

\ 

03  e 
<-< 

QD 

0)       Q. 

■:>   iJ3 

CD  <D  ■ 
—   c- 

c     c 

Q_     H- 1 

A 

Q  3 

O          Q. 

u  ^   E 

in 

A 

Jd  m 

o  cs)    a 

A 

O   (D 

j:  cvj    D 

/\ 

—  CL 

-J    -H     c_ 

— 1                                       /      \ 

1— 1           Q_ 

/  \ 

m 

L 

A 

■  c 

l^ 

•>r  (D 

1 

>  >^c 

QJ 

1     [_ 

—  o  o 

C       T— 

(U     (U  ^-1 

D     —1 

(D   (D  -^ 

1 — 1  _o    a 

(n      D 

c-  3  ca 

(U     QJ    c 

u   E   cji 

U     Q   "H 

It)       CL 
Q)       Ol 

c    .-^ 

3   D~c- 
QOO    (D 

ac    i_  LH 

Q_     CO 

tutuO 

83 


output  is  connected  to  the  vol tage-contro 1 led-f requency 
(VCF)  input  of  the  signal  generator,  and  the  loop  is 
closed  by  the  driver-resonator  system.  The  phase, 
including  calibration  corrections,  is  set  on  the  lock-in 
so  that  the  desired  Y-channel  output  voltage  is  zero  when 
the  pressure  and  velocity  are  in  phase.  The  controller 
has  a  set  point  voltage  fixed  at  zero,  so  the  frequency 
is  controlled  to  maintain  the  resonance  condition  of  the 
pressure  and  velocity  being  in  phase.  The  pressure 
amplitude  control  group  begins  by  adjusting  the  amplitude 
of  the  pressure  signal  to  a  pre-defined  level  with  an 
adjustable  gain  stage.  This  ensures  that  the  loop  gain 
and  stability  margin  are  constant  for  different  drive 
levels.  After  the  pressure  signal  is  adjusted  in 
amplitude,  it  is  rectified  and  filtered,  and  this  DC 
signal  is  used  by  another  integrating  controller  to 
produce  a  control  voltage  for  an  adjustable  gain  control 
Circuit  ( AGC ) .  The  output  of  the  signal  generator  is  then 
passed  through  the  AGC  and  into  the  power  amplifier 
input.  The  loop  is  again  closed  by  the  acoustic  system, 
with  its  pressure  amplitude  thereby  being  regulated. 

Integrating  the  acce lerometer  signal  is  convenient 
since  it  is  the  velocity  that  is  the  important  quantity 
in  a  power  measurement. ..." 

To   finish   the   explanation   we   will   need   to  discuss  the 

velocity   Signal  measurement  since  the  method  used  by  Hofler 

was   not   used   in   our   experiment.   In   our   experiment  we 

measured   the  velocity  amplitude  using  the  lock-in  directly. 

It   was   determined  that  this  would  not  introduce  any  errors 

into   the   power  measurements,  as  previously  thought.  Unlike 

the   previous   apparatus,   our   lock-in   amplifier  uses  true 

Sine   wave   conversion,   resulting   in  complete  rejection  of 

the   harmonic   distortion  in  the  acce lerometer  signal.  Using 

the   measured   dynamic  pressure  and  velocity,  along  with  the 

sensitivities   of   the  transducer  and  acce lerometer ,  and  the 

volume   velocity   calibration   factor   we   can   calculate  an 

accurate   acoustic   power,   as  long  as  the  frequency  control 
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Circuit   is  maintaining  the  phase  difference  between  the  two 
signals  at  zero. 

The  power  protection  circuit  monitors  the  power  amp- 
lifier current  and  will  trigger  a  crowbar  circuit  to  blow  a 
fuse  if  either  of  two  setpoints  are  exceeded.  If  a  set  peak 
current  value  is  exceeded  the  fuse  is  blown,  to  protect 
against  large  fast  transients.  If  a  set  rms  current  value 
is  exceeded  the  fuse  will  also  be  blown,  to  protect  against 
currents  below  the  peak  limit  but  high  enough  to  burn  out 
the  vo  ice  coil, 

C.   EXPERIMENTAL  PROCEDURE 

In  this  section  we  will  discuss  three  basic  areas  of  our 
experimental  procedure.  First  is  the  setup  of  the  apparatus 
for  a  specific  data  run.  Then  we  will  cover  the  procedure 
for  starting  up  and  operating  the  refrigerator.  We  will 
conclude  with  a  discussion  of  the  actual  data  that  is 
recorded.  In  the  operating  section  we  will  cover  two 
different  methods  since  the  procedure  was  changed  after  the 
second  data  run,  for  reasons  that  will  be  explained  later. 

Before  discussing  the  actual  setup  it  should  be  noted 
that  all  temperatures  that  we  recorded  were  measured  using 
commercial  type-E  thermocouple  wire  with  an  accuracy  of 
il.7'-C  or  il%  of  the  reading,  whichever  is  greater.  All  of 
the  thermocouple  wires  are  fed  through  a  Hewlett  Packard 
HP3421A   data   acquisition   control   unit  to  a  computer,  for 
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continuous    display 


on 


the    computer    monitor   during 


operat  ion . 

1 .   Setup  of  Apparatus 

The  setup  of  the  system  for  a  particular  data  run 
consists  of  filling  the  system  with  the  desired  gas, 
pressurizing  the  system  up  to  the  desired  pressure  of  10 
bar  and  starting  up  the  external  vacuum  system. 

The  resonator  is  filled  through  a  connection  at  the 
top  of  the  driver  housing,  however  the  controlling  factor 
is  the  capillary  near  the  driver  surround  that  was 
previously  discussed.  The  pressur izat ion  rate  is  monitored 
on  a  differential  pressure  gauge  since  the  surround  is 
capable  of  sustaining  no  more  than  0.25  bar  of  differential 
pressure  without  permanent  damage.  The  system  is  filled 
using  a  process  of  evacuation  and  purging  that  is  repeated 
twice  before  the  system  is  brought  up  to  operating 
pressure.  This  procedure  is  designed  to  remove  any  gas 
other  than  the  desired  gas  (particularly  air  if  the  system 
has  been  opened)  prior  to  system  operation.  After  the 
second  evacuation  the  system  is  pressurized  to  10  bar.  The 
evacuation  and  purging  process  consumes  the  better  part  of 
a  day  due  to  the  low  flow  rates  that  are  possible  within 
the  allowed  differential  pressure  range. 

To  save  time,  the  external  vacuum  system  can  be 
started  in  parallel  with  the  above  procedure.  Since  the 
diffusion   pump  of  the  vacuum  system  creates  a  great  deal  of 


92 


heat,  an  external  cooling  supply  is  used  to  keep  it  at  a 
reasonable  temperature.  Therefore,  to  start  up  the  vacuum 
system  the  cooling  system  is  lit  off,  the  main  vacuum  pump 
suction  is  closed,  the  vacuum  pump  is  started  and  the  pump 
suction  is  slowly  opened,  to  prevent  rapid 
depressur izat ion  and  damage  to  the  previously  discussed 
super-insulation.   Then   the   diffusion   pump   is  turned  on. 

Once   the   system   is  up  to  pressure  and  the  vacuum  is  below 

_3 
1x10    torr  it  is  ready  to  start  up. 

2.   System  Startup  and  Operation 

To    start   the   refrigerator   the   two   controllers 

(frequency   and   pressure   amplitude)  are  turned  off  and  the 

signal    generator   is   set   up   for   a   sine   wave   at   the 

approximate   resonant   frequency   for   the  resonator  and  gas 

type   being  used.  The  drive  amplitude  is  set  to  zero  and  the 

power   is   turned   on.   The   amplitude   is  brought  up  slowly 

until   a   drive   current   of   about  0.5  amps  is  achieved.  At 

this   point   the  frequency  is  manually  adjusted  to  resonance 

(i.e.   zero   quadrature   output)   and   the  frequency  control 

circuit    is    turned   on.   Then   the   amplitude   is   slowly 

increased   to   the   desired   value   (if   possible)   and   the 

reference   channel  is  adjusted  to  precisely  1.00  volts.  Then 

the    pressure   amplitude   control   circuit   is   turned   on. 

Depending   on  the  desired  pressure  amplitude  and  gas  type  it 

is   not  always  possible  to  reach  the  desired  operating  point 

at   initial   startup.   It   is   therefore  necessary  to  select 
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some  lower  amplitude,  until  the  engine  cools  down  and  the 
desired  amplitude  can  be  reached.  The  limiting  factor  is 
generally  the  piston's  displacement  amplitude.  If  the 
maximum  displacement  amplitude  is  exceeded  significant 
damage  to  the  surround  occurs.  Once  the  desired  amplitude 
is  achieved  the  refrigerator  runs  unattended  until  steady 
state  is  reached.  In  this  case,  steady  state  is  defined  as 
less  than  a  one  to  two  percent  deviation  in  the  cold  end 
temperature  over  a  90  minute  period.  This  is  determined  by 
evaluating  the  continuous  computer  monitor  display  of 
temperatures,  which  shows  approximately  90  minutes  of  data 
history  when  set  up  for  display  at  four  minute  intervals. 
This  is  the  basic  procedure  for  startup  and  setup  for  a  no 
(applied)  heat  load  data  point.  In  general  equilibration 
takes  from  8  to  24  hours  depending  on  the  pressure 
amplitude  setting. 

We  have  used  two  different  methods  for  obtaining 
data  with  heat  loads  applied  to  the  cold  end.  The  first 
method  is  that  used  by  Hofler  (1986),  where  a  designated 
heat  load  is  applied  to  the  heater,  at  the  cold  end,  and 
equ i 1 i-brat ion  is  based  on  temperature  stabilization,  as 
discussed  above  for  the  no  load  case.  This  heat  load  is  set 
by  applying  a  fixed  DC  current  and  voltage  to  the  heater 
coil.  It  was  determined  after  two  data  runs  that  a  better 
method  might  be  available.  This  better  method,  designed  and 
set   up   by   Dr.   Hofler,   invloves   a   temperature  feedback 
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control  network  which  allows  setting  a  desired  steady  state 
cold  end  temperature.  The  controller  then  establishes  the 
current  flow  to  the  heater  that  is  necessary  to  achieve  and 
maintain  that  temperature.  The  equilibration  point  is  no 
longer  based  on  temperature  but  is  instead  based  on  the 
heater  current  deviation  being  less  than  approximately 
i0.5mA  over  a  90  minute  period.  This  corresponds  to  being 
within  better  than  ±1%  over  the  90  minute  period,  for  the 
currents  we  have  measured.  The  control  network  involves  a 
rather  complex  Fortran  computer  program,  written  by  Dr. 
Hofler,  that  reads  the  heater  thermocouple,  at  an  operator 
selected  time  interval,  and  then  performs  a  crude  time 
integral  of  the  difference  between  the  monitored 
temperature  and  the  selected  set  point  temperature.  The 
program's  integrator  output  variable  is  used  to  numerically 
control  the  heater  current  from  the  HP6633A  system  D.C. 
power  supply.  The  program  allows  the  operator  to  select  the 
control  time  interval  (usually  5  seconds),  the  periodicity 
of  display  (usually  4  minutes),  the  desired  temperature,  a 
proportional  constant  in  volts/deg  (usually  0.0),  and  an 
integration  constant  in  vo 1 ts/deg-sec  (usually  0.01).  The 
latter  two  of  these  control  the  stabilization  time  and  the 
degree  of  stabilization.  The  usual  values  that  were  given 
are  those  that  have  shown  the  best  results  during  our 
exper  iments . 
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This  temperature  feedback  control  method  was 
determined  to  be  better  for  the  STAR  project  for  two  basic 
reasons.  First,  it  provided  a  more  rapid  equilibration 
because  the  refrigerator  was  loaded  more  heavily  at  first, 
to  drive  the  temperature  of  the  cold  end  to  the  desired 
temperature,  and  then  the  load  was  let  off  only  as 
necessary  to  maintain  that  temperature.  This  reduced  the 
equilibration  time  to  6  to  8  hours  in  most  cases.  Secondly, 
this  method  conserves  power  due  to  the  reduced  equil- 
ibration time.  This  is  considered  to  be  vital  to  the  STAR 
project  due  to  its  reliance  on  batteries  for  power,  which 
gives  it  a  limited  operating  lifetime. 
3.   Data  Recording 

Once  steady  state  has  been  established  a  data  set  is 
recorded  in  a  notebook.  A  data  set  consists  of  nine 
readings,  as  well  as  date  and  time.  First,  the  temperatures 
of  the  hot  and  cold  heat  exchangers  are  recorded.  There  are 
two  thermocouples  for  each  heat  exchanger,  one  being  at  the 
center  and  one  at  the  wall.  Due  to  an  intermittent  failure 
of  the  thermocouple  at  the  center  of  the  cold  heat 
exchanger  we  have  only  recorded  the  other  three 
temperatures.  Then  we  record  the  drive  frequency,  the  mean 
pressure,  the  amplitude  of  the  dynamic  pressure  signal  and 
the  amplitude  of  the  lock-in  velocity  signal.  Then, 
depending  on  which  method  of  operation  was  used,  we  record 
the   DC   heater   voltage   and  current,  or  just  the  DC  heater 
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current.  For  the  latter  case  we  use  the  known  heater 
resistance  versus  temperature  characteristic  to  calculate 
the  heater  load  (instead  of  voltage  times  current).  After 
the  data  are  recorded  we  change  either  the  dynamic  pressure 
amplitude  or  the  heat  load,  and  wait  for  the  system  to 
equilibrate.  A  standard  data  run  for  a  given  gas  type  and 
stack  involves  8  to  10  data  points  and  takes  about  a  week, 
from  start  to  finish.  The  difference  in  the  number  of  data 
points  is  based  on  the  different  temperatures  achievable  at 
different  pressure  amplitudes,  while  staying  within  the 
piston's  displacement  amplitude  limit. 

In  addition  to  the  above  data,  we  also  record  the 
vacuum  in  the  system  and  the  temperature  and  resonant 
frequency  of  a  gas  analysis  tube,  similar  to  that  described 
by  Polturak,  et.  al.  (1986),  and  Garrett,  et.  al,  (1986), 
that  we  attached  to  the  system.  These  are  monitored  just 
for  completeness  and  are  not  directly  involved  in  any 
performance  evaluation.  With  the  exception  of  one  data  run, 
all  of  our  gases  had  analysis  certifications  for  purity, 
and  these  values  were  used  for  calculations.  In  one  case  we 
mixed  our  own  helium-argon  gas  and  used  the  gas  analyzer, 
built  by  Dr.  Hofler,  to  determine  the  mixture 
concentrat  ion. 

The  final  bit  of  data  that  is  recorded  involves  the 
rate  at  which  the  engine  warms  up,  after  being  fully  cooled 
down   and   turned   off.   This   data  is  analyzed  and  together 
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with  the  total  heat  capacity  of  the  cold  end  of  the  system, 
is  used  to  determine  a  value  for  the  heat  leak  out  of  the 
refrigerator.  The  heat  leak  determined  by  this  method  is 
then  reduced  by  2.4  mw/'-c,  which  is  the  portion  that  has 
been  determined  to  be  due  to  the  thermal  conductance  of  the 
stack  and  gas,  which  is  internal  and  intrinsic  to  the 
stack.  This  heat  leak  is  then  included  in  the  total  heat 
load. 
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V.   DATA  ANALYSIS  AND  RESULTS 

In  this  chapter  we  will  discuss  the  methods  and  equa- 
tions that  were  used  in  analyzing  our  data.  This  will  be 
followed  by  a  numerical  and  graphical  presentation  of 
various  aspects  of  our  experiment.  It  should  be  kept  in 
mind  that  the  purpose  of  these  experiments  was  to  improve 
the  coefficient  of  performance  (COP)  of  the  prototype 
refrigerator,  built  by  Hofler  (1986), 

A.   DATA  ANALYSIS 

All  of  our  data  analysis  was  done  through  the  use  of  a 
Microsoft  Hultiplan  spreadsheet  program  and  graphical 
display  of  the  numbers  that  were  thereby  generated.  This 
method  was  used  because  we  had  voluminous  amounts  of  data 
where  many  of  the  required  calculations  were  repetitive  in 
nature.  Taking  into  account  all  of  the  quantities,  both 
recorded  and  calculated,  we  ended  up  with  a  spreadsheet  of 
25  columns  by  122  rows.  It  should  be  noted  that  many  of 
these  quantities  were  calculated  for  informational  purposes 
and  are  not  directly  related  to  our  final  evaluation  of  the 
refrigerator's  coefficient  of  performance. 

Our  analysis  consists  of  comparing  a  normalized 
(normalized   to   the   Carnot   COP   for   the  same  temperature 
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difference)  coefficient  of  performance,  for  the  refrig- 
erator with  modifications  (COPR),  to  the  same  quantity  for 
the  original  Hofler  refrigerator.  Since  we  have  modified 
both  the  thermoacoust ic  stack  and  the  gas  type,  we  have 
three  basic  areas  for  comparison.  We  will  compare  the 
results  with  different  stacks  and  the  same  gas  type,  with 
different  gases  and  the  same  stack,  and  then  with  different 
stacks  and  gases  combined.  In  order  to  make  these 
comparisons  we  have  created  plots  of  COPK  versus  total  heat 
load  (Qtotal).  In  addition  to  these  we  have  created  plots 
of  temperature  ratio  (Tc/Th)  versus  Qtotal,  to  validate  our 
data.  Based  on  the  experiments  of  Hofler  (1986),  these 
latter  plots  are  expected  to  be  linear. 

In  order  to  create  these  plots  we  need  to  know  the 
values  for  COPR,  Qtotal,  Tc  and  Th .  The  values  for  Tc  and 
Th  are  simply  raw  data  taken  during  the  experiments.  We 
have  used  the  temperature  at  the  center  of  the  hot  heat 
exchanger  as  Th ,  and  the  temperature  at  the  wall  of  the 
cold  heat  exchanger  as  Tc  (due  to  the  intermittent  failure 
of  the  cold  heat  exchanger's  center  thermocouple). 

The  value  for  Qtotal  has  two  components.  The  first  is 
the  column  labelled  QHEAT  in  Appendix  C.  This  is  the  power 
applied  to  the  heater,  located  below  the  cold  heat 
exchanger,  in  watts.  It  is  calculated  by  either  multiplying 
the  heater  DC  voltage  times  the  heater  DC  current,  or  by 
multiplying   the   heater   resistance   by   the   square  of  the 
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heater  DC  current  (P=VI  or  ?=l^'R  for  DC).  The  voltage  and 
current  values  are  raw  data.  The  resistance  values  were 
taken  from  a  plot  of  Tc  vs.  R  created  from  the  two  sets  of 
data  where  we  recorded  both  heater  voltage  and  current 
(R=V/I).  This  plot  is  shown  as  Figure  V- 1 .  Both  of  these 
methods  were  used  as  discussed  in  Chapter  IV.  The  second 
component  of  Qtotal  is  the  heat  leak,  which  is  labelled 
HEATLK  in  Appendix  C.  As  discussed  in  Chapter  IV,  this 
value  is  determined  by  analysis  of  the  rate  at  which  the 
refrigerator  warms  up  once  the  power  is  turned  off.  Based 
on  this  warmup  data  and  using  the  heat  capacity  of  the  cold 
portion  of  the  resonator  a  value  for  the  total  heat  leak  to 
the  cold  end  of  the  refrigerator  is  determined  (in  mW/*c). 
This  number  is  then  reduced  by  2.4  mW/*c,  which  is  the 
thermal  conduction  of  the  stack  and  the  gas.  The  resulting 
quantity  is  the  conduction  external  to  the  stack  and  is 
accuarate  to  within  8%.  The  value  for  Qtotal  is  then  the 
sum  of  QHEAT  +  HEATLK.  The  error  in  QHEAT  is  negligible, 
but  the  error  in  HEATLK  is  not.  The  error  in  Qtotal  is 
about  8%  when  QHEAT  is  zero  and  is  about  1.5%  to  2%  for 
values  of  Qtotal  near  heat  loads  of  optimum  efficiency. 

The  remaining  quantity  necessary  for  the  calculation  of 
COP  is  the  acoustic  power  (or  work).  We  use  a  relation 
Similar  to  that  derived  by  Hofler  (1986).  The  difference 
being  that  we  no  longer  use  the  velocity  signal  from  the 
bandpass   filter;  using  the  velocity  signal  from  the  lock-in 
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instead.   Also,   we   have   added  two  correction  terms  to  the 
equation.  We  are  using  an  equation  of  the  form: 

w= 


where   Hp    is   the  dynamic  pressure  transducer  sensitivity, 

M.,   is   the   volume  velocity  sensitivity,  Vp  is  the  measured 

pressure   signal,   V, .   is  the  measured  velocity  signal,  (  1- 

t.)       is   an  accelerometer  non-linearity  correction  term,  and 
A 

€yy   is   a   volume   velocity   calibration   correction   term. 

-5 
t.  =  (7. 4x10     secVcm)a,   where   a   is   the   acceleration  in 
A 

_3 

cm/sec^   (xl0   ),   as   shown   in   column   13   of  Appendix  C. 

^   =0. 05( f-550) ,   where   f   is   the   frequency   in  hertz,  as 

shown   in   column   4   of   Appendix  C.  The  ratio  H,j/Hp  is  the 

result   of   the   small   cavity   volume  velocity  calibration, 

discussed   previously,   and   is  therefore  just  a  number  that 

we   put   into   the   equation  prior  to  calculation.  Since  the 

refrigerator   was   disassembled   on  several  occasions  during 

our   experiments   this   calibration   was   performed   several 

times.   The   calibration   value   changes   over   time   as  the 

driver   is  used  and  therefore  this  number  M../Mp  is  corrected 

in   our  equations  for  all  data  taken  after  each  calibration. 

The   quantity   Mp   is  also  a  calculated  number  that  is  input 

into   the   equations,   since  it  is  just  the  dynamic  pressure 

transducer   sensitivity.   The   actual  equations  that  we  have 

used  are  of  the  following  form: 

W=59810(  l-7.4xl0~^a)Vp(0.003Vj^j  )/(  1839.  8-0.  05  (  f -550  )  )  , 
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where   59810   is   (1/Mp)^   including   unit   conversions,  (  1- 

7,4x10   a)   is   a   correction   term   for   accelerometer  non- 

liearities,   as  discussed  above,   (0.003V,,)  is  the  measured 

velocity   signal  (with  unit  conversions  in  the  coefficient), 

1839.8    is   (H.j/Hp)   the   result   of   our   volume   velocity 

calibration,   and   (0. 05( f-550) )   is   a   correction   to   the 

volume   velocity  term  due  to  the  calibration  being  frequency 

dependent,   as  discussed  above.  The  accelerometer  correction 

(1-^.)   is   about   2%   at  most,  and  the  frequency  correction 

(t   )   is   about   1%   at   most.  The  frequency  correction  was 

obtained   for   the   case   of  pure  helium.  However,  given  the 

much   larger  frequency  range  of  the  gas  mixtures,  the  amount 

of   correction   in   some   cases   is   too   small,   giving  a 

calibration   error   of   about  1%,  The  total  worst  case  error 

for  W  is  about  5%  with  3.5%  being  a  more  typical  value. 

Using   the   Qtotal   and  W  values  calculated  by  the  above 

methods   we   get  COP=Qtotal/W,  To  get  COPR  we  need  to  divide 

this   by  the  Carnot  COP  for  the  same  temperature  difference. 

The   Carnot   COP   is  known  to  be  COP      ^=Tc/(Th-Tc ) ,  which 

carnot 

means  COPR  can  be  derived  as  COPR=COP( (Th/Tc ) -1 ) .  This  is 
the  equation  we  have  used  to  determine  COPR,  On  our 
spreadsheet,  in  Appendix  C  this  is  the  C0PR2  column.  The 
worst  case  error  in  the  absolute  values  for  COPR  may  be  as 
high  as  15%  however,  typical  values  are  about  5%.  Also,  the 
relative  error  involved  in  comparing  data  taken  under 
different  conditions  is  probably  considerably  less  than  5%, 
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We  now  have  all  of  the  values  necessary  to  evaluate  the 
effects  of  our  modifications  on  the  refrigerator's 
coefficient  of  performance. 

B.   RESULTS  OF  EXPEKIMENT 

As  stated  earlier,  we  will  break  this  discussion  into 
three  major  sections.  First  we  will  discuss  the  results  of 
our  stack  modifications.  This  will  be  followed  by  a 
discussion  of  performance  with  different  gases  and  the  same 
stack.  We  will  conclude  with  the  overall  improvements  using 
our  best  stack/gas  combination.  Also,  a  brief  discussion 
concerning  the  results  with  different  pressure  amplitudes 
and  different  resonator  lengths  will  be  provided. 

Since  we  have  data  from  two  new  stacks  and  two  gas 
mixtures,  as  well  as  two  different  pressure  amplitudes  and 
two  resonator  lengths,  available  for  comparisons,  we  must 
be  selective  in  our  presentation,  so  as  not  to  have  an 
excessive  number  of  plots  that  show  the  same  result.  To 
fulfill  this  goal  we  will  present  one  or  two  sample  sets  of 
curves,  with  analysis,  for  each  major  category  of  interest, 
and  provide  numerical  results  for  the  other  cases  within 
that  category,  where  apprpriate.  In  general  a  set  of  curves 
includes  a  plot  of  COPE  versus  Qtotal  and  a  plot  of  Tc/Th 
versus  Qtotal. 

Prior  to  presenting  any  of  our  figures  we  shall  intro- 
duce  the   notation   that  is  used  in  the  titles  and  legends. 
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All  of  our  plots  are  designed  to  provide  data  comparisons 
and  as  such  each  contains  more  than  one  curve.  The  legend 
lists  the  specific  curves  using  the  following  notation: 

GAS  TYPE— STACK—RESONATOR  LENGTH— PRESSURE  AMPLITUDE, 
where  abbreviations  have  been  used  in  all  cases.  For  gas 
types  we  have  helium,  helium-argon,  and  helium-xenon  which 
are  represented  by  He,  HeAr,  and  HeXe,  respectively.  For 
the  stacks  we  have  old,  new,  and  dual  which  are  represented 
by  0,  N,  and  D,  respectively.  The  resonator  length  is 
designated  by  s  or  1  for  short  or  long  and  the  pressure 
amplitude  is  represented  by  the  percent,  relative  to  the 
mean  pressure  (Po/Pm),  by  the  number  2  or  3, 
1 .   Stack  modification  results 

To  illustrate  the  effects  of  our  stack  modifications 
we  have  elected  to  compare  the  data  for  pure  helium,  at  2% 
pressure  amplitude,  with  the  short  resonator  installed. 
Figure  V-2  shows  the  temperature  ratio  (Tc/Th)  versus  total 
heat  load  (Qtotal)  curves  for  all  three  stacks  under  these 
conditions.  It  is  fairly  clear  from  the  curves  that  the 
data  are  essentially  linear  for  all  three  cases,  as 
expected.  The  curves  also  show  that  for  a  given  temperature 
ratio  the  old  stack  can  support  more  heat  load  than  the 
dual  stack  and  the  new  stack  can  support  more  heat  load 
than  either  of  the  others.  This  statement  can  be  reversed, 
such  that  for  a  given  heat  load  the  temperature  spanned  is 
greatest   for   the   new   stack  and  lowest  of  the  dual  stack. 
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Figure  V-2.  Temperature  Ratio  vs.  Qtotal,  For  Helium,  Short  Resonator 
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These  curves  also  show  that  our  data  are  following  the 
expected  trends  for  the  refrigerator,  in  general.  Figure  V- 
3  shows  COPR  versus  Qtotal  for  all  three  stacks.  It  is 
clear  from  this  Figure  that  the  COP  is  significantly 
improved  with  the  new  stack  and  vastly  degraded  with  the 
dual  stack.  Evaluating  the  difference  at  the  peaks  of  the 
curves,  we  see  that  the  dual  stack  has  a  peak  COPR  of  0.034 
while  the  old  stack  has  a  peak  at  0.067  and  the  new  stack 
peaks  at  0.078.  This  shows  a  16%  improvement  for  the  new 
stack  and  a  51%  reduction  for  the  dual  stack.  Conducting 
this  same  analysis  for  the  remainder  of  our  data  sets 
yields  an  average  improvement  of  14%  with  the  new  stack  and 
an  average  reduction  of  57%  with  the  dual  stack. 

It  is  clear  from  these  results  that  the  dual  stack 
has  not  performed  as  predicted  while  the  new  stack  has 
provided  significant  improvements.  It  is  believed  that  the 
basic  theory  used  to  design  the  dual  stack  is  valid,  in-so- 
far  as  the  concept  is  concerned.  However,  it  is  obvious 
that  we  have  underdesigned  the  stack  by  either  omitting  an 
important  factor  or  over  simplifying.  At  the  time  of  design 
and  construction  of  the  stack,  it  was  our  belief  that  plate 
spacing,  so  long  as  it  was  not  too  small,  would  not  have 
any  detrimental  effects.  Based  on  this  belief,  and  an 
initial  evaluation  of  the  new  stack's  performance,  we 
decided  to  build  the  dual  stack  with  the  hot  end  having  the 
same   fishing   line   size   as   that   used   in  the  new  stack. 
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Figure  V-3.  COPR  vs.  Qtotal  for  Helium,  Short  Resonator 
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However,  with  the  dual  layer  at  the  hot  end  and  the  thick 
fishing  line  at  the  cold  end,  we  ended  up  with  a  20% 
increase  in  plate  spacing  at  the  hot  end  and  a  163% 
increase  in  plate  spacing  at  the  cold  end  (relative  to  the 
spacing  in  the  new  stack).  Based  on  the  extremely  poor 
results  with  the  dual  stack  it  appears  that  the  magnitude 
of  the  plate  spacing  is  more  critical  than  we  had 
originally  thought.  The  theory,  presented  in  Chapter  II, 
shows  that  a  plate  spacing  of  2  to  4  thermal  penetration 
depths,  with  the  entire  resonator  cross-section  filled  with 
plates  is  desired  for  the  best  heat  pumping 
characteristics.  However,  our  analysis  of  the  theory  did 
not  predict  this  massive  drop  in  COP  with  a  larger  plate 
spacing.  To  some  extent  this  plate  spacing  effect  is 
illustrated  by  the  new  stack's  performance,  as  well. 
Although  our  intuitive  reasoning  for  the  design  of  the  new 
stack  predicted  improved  performance,  the  results  are 
better  than  expected.  Combining  this  with  the  extremely 
poor  performance  of  the  dual  stack,  and  the  fact  that  the 
new  stack  has  the  smallest  plate  spacing  of  the  three 
stacks,  we  can  see  that  plate  spacing  is  much  more 
important  than  previously  thought.  Rough  calculations  for 
helium  show  that  the  average  plate  spacing,  in  number  of 
thermal  penetration  depths  is  3.3  for  the  new  stack,  3.5 
for  the  old  stack,  4.0  for  the  hot  end  of  the  dual  stack, 
and  8.8  for  the  cold  end  of  the  dual  stack. 
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2.   Results  for  Gas  Mixtures 

For  this  analysis  we  have  elected  to  present  two 
sets  of  data,  both  of  which  are  for  the  new  stack.  This  is 
primarily  due  to  the  fact  that  the  new  stack  is  the  only 
one  for  which  we  have  any  helium-argon  data.  Also,  based  on 
the  stack  analysis,  the  new  stack  is  the  best  one 
available.  For  this  analysis  we  will  discuss  the  results 
for  2%  pressure  amplitude  and  both  the  short  and  long 
resonator.  The  long  resonator  consists  of  a  one  inch  spacer 
that  lengthens  the  resonator  and  thereby  changes  the 
position  of  the  stack  realtive  to  the  acoustic  standing 
wave. 

Figures  V-4  and  V-5  show  the  plots  of  Tc/Th  versus 
Qtotal  for  the  short  and  long  resonators,  respectively.  It 
is  clear  from  these  curves  that  all  of  the  data  are  linear, 
as  expected.  It  is  interesting  to  note  that  in  Figure  V-2 
all  of  the  curves  for  pure  helium  were  essentially 
parallel,  even  though  the  data  was  for  different  stacks. 
Now,  in  Figures  V-4  and  V-5  we  see  that  the  curves  are  no 
longer  parallel  for  the  different  gases,  even  though  the 
data  are  all  for  the  same  stack.  This  is  believed  to  be 
related  to  the  fact  that  the  heat  pumping  capacity  is 
related  to  the  PU  product  (pressure  times  volume  velocity). 
For  helium  gas  at  constant  pressure  amplitude  the  PU 
product  will  vary  in  a  consistent  manner,  with  respect  to 
temperature,   regardless   of   stack   design.   For   different 
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Figure  V-4.  Temperature  Ratio  vs.  Qtotal,  New  Stacl<,  Short  Resonator,  All  Gases 
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Figure  V-5.  Temprtature  Ratio  vs.  Qtotai,  New  Stack,  Long  Resonator,  All  Gases 
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gases  however,  the  velocities  will  vary  because  of 
different  sound  speeds,  and  therefore  even  at  constant 
pressure  amplitude  the  PU  product  will  vary  differently, 
producing  curves  with  different  slopes. 

The  plots  of  COPR  versus  Qtotal  for  the  short  and 
long  resonator  are  shown  in  Figures  V-6  and  V-7, 
respectively.  It  is  clear  from  all  of  these  curves  that 
helium-argon  produces  a  certain  degree  of  improvement  in 
COP,  over  pure  helium,  and  that  helium-xenon  produces  the 
highest  COPR.  For  the  short  resonator  curves  of  Figure  V-6 
we  see  that  the  peak  COPR  values  are  approximately  0.078 
for  helium,  0.102  for  helium-argon,  and  0.129  for  helium- 
xenon.  From  these  values  we  see  that  helium-argon  gives  a 
31%  increase  over  helium,  while  helium-xenon  provides  a  65% 
improvement  over  helium  and  a  27%  improvement  over  helium- 
argon.  It  should  be  noted  that  the  peak  COPR  values  are 
shifted  to  lower  heat  load  values  for  the  gas  mixtures. 
This  is  believed  to  be  due  to  a  reduction  in  this  PU 
product  that  was  discussed  above.  Similar  calculations  for 
other  stacks  and  other  pressure  amplitudes  show  an  average 
improvement  in  the  peak  COPR  values  of  81%,  when  comparing 
helium-xenon  to  helium.  There  is  no  other  short  resonator 
data  available  for  helium-argon.  The  long  resonator  curves 
of  Figure  V-7  have  been  included  primarily  because  they 
illustrate  the  peak  COPR  value  that  has  been  measured  to 
date   of   0.197   for  helium-xenon.  Although  the  helium-argon 
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Figure  V-6.  COPR  vs.  Qtotal,New  Stack,  Short  Resonator,  All  Gases 
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Figure  V-7.  COPR  vs.  Qtotal,  New  Stack,  Long  Resonator,  All  Gases 


curve  appears  somewhat  incomplete  we  feel  safe  in  assuming, 
based  on  all  of  our  other  curves  and  data,  that  the  curve 
is  within  a  few  percentage  points  of  its  peak  value  at  the 
last  plotted  data  point.  We  will  therefore  use  that  datum 
for  our  peak  value.  These  curves  show  the  same  general 
trends  as  those  of  Figure  V-6,  however  the  magnitudes  of 
the  improvements  are  somewhat  reduced.  For  peak  values  of 
0,132  for  helium,  0.1595  (from  data  of  Appendix  C)  for 
helium-argon  and  0.197  for  helium-xenon,  we  see  that  helium- 
argon  only  produces  a  21%  improvement  over  helium  and 
helium-xenon  produces  a  49%  improvement  over  helium  with  a 
24%  improvement  over  helium-argon.  It  is  believed  that  the 
improvement  percentages  with  the  long  rresonator  are 
reduced  because  the  long  resonator  itself  produces  about  a 
62%  improvement  over  the  short  resonator,  thereby  placing 
the  refrigerator  closer  to  a  maximum  possible  or  limiting 
COPR  value.  Thus  the  absolute  improvements  we  achieved  have 
a  smaller  effect,  percentage  wise. 
3,   Overall  Improvement 

In  this  section  we  present  the  overall  result  of  the 
thesis.  That  is,  the  overall  improvement  in  the  COP  of  the 
prototype  refrigerator,  taking  into  account  both  the  stack 
modifications  and  the  binary  gas  mixtures.  To  illustrate 
this  case  we  have  elected  to  compare  the  results  with  the 
old  stack  and  pure  helium  at  2%  pressure  amplitude  to  those 
achieved   with  the  new  stack  and  helium-xenon  at  2%  pressure 
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amplitude.  These  sets  of  data  were  selected  because  the 
original  refrigerator  used  only  the  old  stack  and  pure 
helium,  while  our  best  results  were  achieved  with  the  new 
stack  and  helium-xenon.  We  have  not  used  the  long  resonator 
data,  even  though  that  gives  our  highest  COPR,  because  we 
did  not  take  any  long  resonator  data  with  the  old  stack. 
However,  for  a  number  for  the  old  stack,  long  resonator, 
pure  helium  setup  we  can  use  the  0.126  peak  COPR  given  by 
Hofler(1986)  in  his  doctoral  dissertation.  Comparing  this 
to  our  best  case  peak  COPR  of  0. 197  we  get  a  56% 
improvement. 

For  the  comparison  of  our  experimental  data  we  will 
present  the  COPR  versus  Qtotal  plot  only,  since  the  Tc/Th 
versus  Qtotal  plots  have  already  been  discussed  for  this 
data.  Figure  V-8  shows  our  plot  for  overall  improvement. 
With  peak  COPR  values  of  0.129  for  the  helium-xenon  and 
0.067  for  the  helium  we  see  that  the  overall  improvement  in 
the  refrigerator's  coefficient  of  performance,  due  to  our 
modifications,  is  93%. 

4.   Other  Effects  of  Interest 

There  are  two  other  areas  of  interest  that  have  been 
mentioned  briefly  during  our  discussions  of  results.  These 
are  the  effects  produced  by  changing  resonator  length  and 
those  produced  by  changing  the  pressure  amplitude. 

In  Figure  V-9  we  present  the  data  for  all  three 
gases   in  the  new  stack  with  both  short  and  long  resonators. 
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Figure  V-8.  COPR  vs.  Qtotal,  Showing  Overall  Performance  Improvement 
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Figure  V-9.  COPR  vs.  Qtotal,  Showing  Resonator  Effects 


120 


From  this  we  can  see  that  just  the  shift  from  short  to  long 
resonator  produces  increases  in  COPE  as  follows:  70%  for 
helium,  59%  for  helium-argon  and  52%  for  helium-xenon.  This 
rather  large  increase  in  COP  is  due  primarily  to  the  shift 
in  the  position  of  the  stack  relative  to  the  standing  wave. 
With  the  long  resonator  installed,  the  hot  end  of  the  stack 
is  moved  closer  to  the  velocity  antinode  and  the  heat 
transfer  power  is  increased  significantly,  thereby 
producing  the  large  increase  in  the  COP  and  the  observed 
shift  of  the  COPR  peak  towards  higher  heat  loads. 

In  Figures  V-10  and  V- 1 1  we  present  a  comparison 
between  both  gas  types  and  pressure  amplitudes.  In  Figure  V- 
10  we  have  helium-xenon  with  the  short  resonator,  both 
pressure  amplitudes  and  all  three  stacks.  In  Figure  V-11  we 
have  all  of  the  same  curves  for  pure  helium.  In  all  of 
these  curves  it  is  obvious  that  at  higher  pressure 
amplitudes  the  peak  COPR  value  has  shifted  toward  higher 
heat  loads.  Based  on  this  and  the  discussion  of  Figures  V-4 
and  V-5  it  is  clearly  shown  that  heat  transfer  power  (or 
heat  pumping  capacity)  increases  with  pressure  amplitude 
(due  to  the  PU  product  increasing).  It  is  interesting  to 
note  that  the  peak  COPR  value  is  in  fact  lower,  at  higher 
pressure  amplitudes,  for  all  of  the  cases  for  which  we  have 
data,  with  the  exception  of  the  pure  helium,  short 
resonator,  old  stack  case.  In  his  doctoral  dissertation, 
Hofler   reported   that   the   peak   COPR  value  increased  with 
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Figure  V-10.  COPR  vs.  Qtotal  for  HeXe,  Short  Resonator,  All  Stacks  and  Amplitudes 
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Figure  V-11.  COPR  vs.  Qtotal  for  He,  Short  Resonator,  All  Stacks  and  Amplitudes 
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pressure  amplitude  for  all  of  his  data.  In  fact  our  data 
with  his  stack  and  helium  also  follow  this  trend. 
Additionally,  Hofler  (1986)  reported:  "there  is  also  some 
indication  that  the  magnitude  of  these  maxima  approach  some 
constant  value  at  higher  amplitudes."  Our  results  for  gas 
mixtures  and  other  stack  designs  tend  to  contradict  this, 
and  lead  toward  a  rather  alarming  conclusion  that  there  is 
some  relatively  low  pressure  amplitude  at  which  the  COPR 
maxima  peak  and  above  which  the  COPR  maxima  drop  off  again. 
This  seems  to  indicate  that  there  may  be  some  non- 
linearity,  such  as  turbulence,  that  limits  the  best  COP's 
to  relatively  low  amplitudes. 
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VI.   CONCLUSIONS  AND  RECOHHENDATIONS 

The  primary  purpose  of  this  thesis  is  to  provide 
modifications  for  the  Space  Thermoacoust ic  Refrigerator 
(STAR)  that  will  improve  its  thermodynamic  efficiency, 
specifically  its  coefficient  of  performance  (COP).  Through 
the  use  of  a  helium-xenon  binary  gas  mixture  and  a  new 
stack  design,  with  less  string  and  reduced  plate  spacing, 
we  have  produced  a  93%  increase  in  the  peak  value  of  COP. 

Along  with  this  improvement  we  have  raised  several 
questions  that  are  as  yet  unanswered.  Based  on  our  results 
with  a  long  resonator,  relative  to  those  with  a  short 
resonator,  it  would  appear  that  there  is  some  optimum 
resonator  length  and  design  that  must  be  found.  Further 
testing  with  different  resonator  spacer  thicknesses,  must 
be  performed  to  resolve  this  issue.  The  unexplained  drop 
in  peak  COP  relative  to  Carnot  (COPR),  with  increasing 
pressure  amplitude  opens  the  door  for  more  testing  with 
different  setups  and  many  pressure  amplitudes.  It  is  clear 
that  the  binary  gas  mixture  of  helium-xenon  is  superior 
under  all  of  the  conditions  that  we  have  observed.  However, 
the  failure  of  the  dual  stack  and  the  better  than  expected 
results  with  the  new  stack  leave  the  final  stack  design  as 
an  open  item.  It  is  clear  that  another  dual  stack  and  a 
different   version   of   the   new  stack  should  both  be  built, 
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with  smaller  plate  spacings,  to  determine  the  optimum  stack 
design.  It  is  the  answers  to  these  question  that  will  yield 
the  design  capable  of  producing  the  maximum  thermodynamic 
efficiency  for  the  Space  Thermoacoust ic  Refrigerator 
(STAK). 

Although  this  project  deals  specifically  with  the  space 
application  of  our  refrigerator,  there  is  no  reason  that  it 
should  be  limited  to  that  application.  Due  to  the  lack  of 
moving  parts  that  are  present  in  most  other  cooling  engines 
(i.e.  sliding  seals)  this  unit  should  be  extremely 
reliable.  It  also  has  reasonable  efficiency  relative  to 
most  commercial  refrigeration  devices.  The  fact  that  it 
does  not  use  chlorof lourocarbons  is  a  great  advantage, 
especially  with  the  close  scrutiny  that  this  ozone 
destroying  substance  is  currently  receiving.  Therefore, 
this  type  of  refrigeration  device  could,  and  should  be 
developed  for  commercial  use  and  should  be  considered  as  a 
possible  candidate  for  virtually  any  type  of  cooling  need. 
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APPENDIX  A.   PRANDTL  NUMBER  PROGRAM 

This  appendix  provides  a  program  and  documentation  for 
the  calculation  of  Prandtl  numbers  for  binary  gas  mixtures. 
The  program,  written  by  the  author,  is  designed  to 
calculate  Prandtl  numbers  for  all  concentrations  of  a 
binary  mixture  of  Helium  and  one  other  gas.  The  program 
is  written  in  Fortran  for  use  on  the  Naval  Postgraduate 
school  mainframe  computer,  with  a  Watfor-77  compiler.  Two 
output  files  are  created  when  the  program  is  run.  The 
first,  called  "PRANTL  OUTPUT  Al",  is  a  listing  of  Prandtl 
numbers  versus  helium  fraction  for  further  use  in  a 
plotting  program,  to  create  Prandtl  number  versus  helium 
fraction  curves  for  the  gas  mixture  selected.  The  second, 
called  "PRANTL  RESULTS  Al"  is  a  tabular  listing  of  Prandtl 
number,  thermal  conductivity  (in  cal/cm-sec-K) ,  and  dynamic 
viscosity  (in  gm/cm-sec)  versus  helium  fraction  for  the 
selected  gas  mixture.  The  listings  go  from  pure  helium  to 
zero  percent  helium  in  one  percent  intervals  to  provide 
easy  access  to  any  specific  mixture.  All  of  the  theory  and 
equations  necessary  to  create  the  program,  other  than 
knowledge  of  a  programming  language,  can  be  found  in 
"MOLECULAR  THEORY  of  GASES  and  LIQUIDS"  by  Hirschf e Ider , 
Curtiss,  and  Bird,  John  Wiley  &  Sons,  Inc,  New  York,  1954. 
The   primary  source  was  chapter  eight  of  the  text,  hereafter 
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referred  to  as  Hirschf elder ,  although  section  and  equation 
numbers  will  be  specified  where  necessary  for 
clar i  f icat ion. 

As  written,  this  program  was  designed  for  personal  use, 
exclusively  to  derive  the  data  necessary  to  create  the 
previously  mentioned  curves  of  Prandtl  number  versus  helium 
fraction.  It  has  been  modified  to  include  the  second 
tabular  output  for  completeness  and  further  use  by  others 
in  this  field.  The  program  is  currently  written  to  require 
hard  input  and  therefore  requires  editing  of  the  main 
program  for  each  change  of  gas  type  or  temperature.  In 
addition  to  satisfying  the  specified  needs,  this  method  was 
used  due  to  there  being  thirteen  input  parameters  that  vary 
With  gas  type  and  temperature,  five  of  which  require  table 
lookups  in  Hirschf elder .  This  program  could  be  modified  for 
incorporation  into  a  larger  program  as  a  subroutine, 
provided  the  specific  gas  types  and  temperatures  of 
interest  are  completely  specified  and  all  of  the  input 
parameters  are  looked  up  and  incorporated  into  a  data  array 
for  reading.  This  has  not  been  done  to  date  as  the  current 
useage  did  not  require  these  features. 

As  previously  stated  the  primary  source  for  this  program 
is  chapter  eight  of  Hirschf elder ,  which  provides  the 
equations  necessary  to  calculate  various  transport 
coefficients,  such  as  the  coefficient  of  viscosity  and  the 
coefficient    of    thermal    conductivity.    The    equations 
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presented  in  chapter  eight  rely  heavily  on  the  elementary 
theory  of  transport  phenomena  presented  in  section  1.2  of 
the  text  as  well  as  the  rigorous  kinetic  theory  for 
monatomic  gases  presented  in  chapter  seven.  The  intricate 
details  of  this  theory  will  not  be  presented  here,  however 
a  discussion  of  the  methods  used,  including  the  main 
assumptions  and  restrictions,  will  be  provided  as  well  as  a 
definition  of  each  of  the  terms  or  symbols  in  the  program. 

The  development  of  the  kinetic  theory  for  gases  is  based 
on  knowledge  of  the  distribution  function  that  represents 
the  number  of  molecules  of  a  specific  species,  in  a  unit 
volume  element  about  a  point  in  space,  with  velocities  in  a 
unit  range  about  a  specified  velocity,  at  a  given  instant 
in  time.  Since  the  usual  definitions  for  transport 
coefficients  only  apply  under  conditions  of  equilibrium  (or 
only  slightly  different  from  equilibrium)  this  restriction 
applies  to  the  results.  Under  this  limit  the  distribution 
function  is  nearly  Haxwellian  and  is  solvable  by  a 
perturbation  method  developed  by  Chapman  and  Enskog.  The 
solutions  are  then  used  to  obtain  expressions  for  transport 
coefficients  in  terms  of  a  set  of  integrals,  Omega(l,s), 
which  involve  explicitly  the  dynamics  of  molecular 
encounters.  The  Chapman-Enskog  theory  relies  on  several 
assumptions  and  therefore  has  limited  applicability.  Each 
assumption  will  be  briefly  discussed  and  related  to  the 
specific  situation  of  our  experiment. 
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Only  binary  collisions  are  considered,  which  limits  us 
to  densities  low  enough  that  three  body  collisions  are 
unimportant.  The  theory  predicts  that  viscosity  is 
independent  of  pressure  at  constant  temperature. 
Experiments  on  Nitrogen  have  shown  that  only  a  4%  deviation 
in  viscosity  occurs  as  pressure  is  varied  form  1  to  60 
atmospheres.  Our  operating  condition  of  10  atmospheres  is 
well  within  this  and  should  therefore  closely  follow  the 
theory  in  this  respect. 

The  theory  further  limits  itself  to  Classical  Mechanics 
and  thereby  precludes  low  temperature  situations  where 
quantum  effects  become  significant.  These  quantum  effects 
are  less  than  1%  for  helium  above  200  Kelvin  (and  even 
smaller  for  heavier  isotopes).  The  coldest  temperature  we 
have  achieved  to  date  has  been  193  Kelvin  under  a  no  load 
condition  so  our  normal  operations  are  all  above  200  Kelvin 
and  we  are  not  concerned  with  quantum  effects. 

The  Chapman-Enskog  theory  uses  a  series  approximation, 
to  the  above  mentioned  distribution  function,  in  solving 
the  Boltzmann  equation  and  then  uses  only  the  first 
approximation,  which  is  only  valid  when  the  gradients  of 
the  physical  quantities  (i.e.  density,  velocity, 
temperature)  are  small.  The  changes  of  these  properties 
over  one  mean  free  path  must  be  small  with  respect  to 
unity.  Based  on  the  mean  free  paths  for  gases  at  pressures 
above   one   atmosphere   this   condition   is  satisfied  except 
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under  conditions  of  extreme  gradients  such  as  those  present 
in  shock  waves.  Since  we  do  not  have  any  shock  waves  we  are 
not  restricted  by  this. 

The  theory  further  assumes  that  the  dimensions  of  the 
containing  vessel  and  any  obstacles  within  it  are  large 
compared  to  the  mean  free  path,  which  provides  no 
limitations  for  our  gases  and  our  pressures.  Although  the 
Chapman-Enskog  theory  applies  strictly  to  monatomic  gases 
the  results  have  been  shown  to  be  good,  even  for  polyatomic 
gases,  provided  that  the  molecules  are  not  to  non- 
spherical.  Finally,  the  development  of  the  theory  results 
in  the  set  of  integrals  Omega(l,s),  which  are  present  in 
all  of  the  equations  for  transport  coefficients.  These 
integrals,  in  turn  are  dependent  upon  the  potential 
function  that  represents  the  molecular  interactions. 
Throughout  the  derivation  in  Hirschf elder ,  and  for  all  of 
the  equations  used  here,  the  Lennard-Jones  (6-12)  potential 
is  assumed. 

Since  some  of  the  other  terms  and  symbols  depend  upon 
the  Omega  terms  in  the  program,  we  will  start  with  a 
discussion  of  where  the  Omega  terms  come  from.  These 
Omega(l,s)  terms  represent  a  set  of  integrals,  linear 
combinations  of  which  are  used  to  evaluate  the  more  complex 
bracket  integrals  arrived  at  in  solving  for  the  expansion 
coefficients  of  the  Sonine  polynomials.  The  transport 
coefficients   are   expressed    in    terms   of   the   Sonine 


polynomial  expansion  coefficients  and  can  then  be  expressed 
in  terms  of  these  OmegaCl.s)  integrals.  In  Hirschfelder 
there  is  also  an  Omega-star ( 1 , s )  symbol  which  represents 
the  deviation  of  a  particualr  molecular  model  (in  our  case 
the  Lennard-Jones  (6-12)  potential)  from  the  idealized 
rigid  sphere  model  of  molecular  interactions.  All  of  the 
Omegas  in  the  program  represent  an  Omega-star  value.  The 
Omega  integrals  are  functions  of  the  reduced  temperature 
(T*=kT/n  and  are  tabulated  in  table  I-H  of  Hirschfelder. 

Where  more  than  one  variable  of  a  given  type  exists, 
such  as  sigmal,  sigma2,  and  sigml2,  the  1  and  2  represent 
the  value  of  the  quantity  for  pure  gas  1  or  pure  gas  2 
(which  is  always  helium  for  our  purposes),  respectively. 
The  12  represents  the  value  for  gases  1  and  2  together. 
This  is  true  for  all  cases  except  ETA12  and  LAM12-  These 
two  variables  represent  artificial  quantities  for  a 
hypothetical  pure  substance  with  molecules  of  molecular 
weight  2(H1) (H2)/(H1+H2)  that  interact  according  to  a 
potential  curve  specified  by  interaction  parameters  sigmal2 
and  n2.  These  parameters  are  defined  in  terms  of  the 
Lennard-Jones  (6-12)  potential  in  chapter  one  of 
Hirschfelder,  where  sigma  is  the  distance  of  closest 
approach  (in  angstroms)  of  two  molecules  which  collide  with 
zero  initial  relative  kinetic  energy  and  t  is  the  maximum 
energy  of  attraction  of  the  two  molecules.  Table  I-A  of 
Hirschfelder   provides   values  for  sigma  and  i/k    for  various 
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substances,  where  k  is  the  Boltzmann  constant  and  i/k  has 
units  of  Kelvin.  The  parameters  ^12  and  siginl2  are  arrived 
at  through  empirical  combining  laws  which  relate  force 
constants  between  unlike  molecules  to  those  between  like 
molecules.  These  laws  are  sigm  12=(  1/2) ( sigmal+slgma2)  and 
U2=(U*^2)  .  The  quantities  A12  and  B12  are  merely 
frequently  encounterd  ratios  of  different  Omegas  that  have 
been  tabulated  separately  as  functions  of  T*  in  table  I-N 
of  Hirschf elder. 

The  remaining  terms  or  symbols  in  the  program  are  fairly 
straightforward,  based  on  the  above  explanations  and  each 
type  will  now  be  defined: 

PR=prandtl  number 

HE=helium  fraction 

LAM=thermal  conductivity  coefficient  (in  cal/cm-sec-K) 

ETA=coef f icient  of  dynamic  viscosity  (in  gm/cm-sec) 

H=molecular  weight 

CP=specific  heat  capacity  at  constant  pressure 
( in  cal/gm-K) 

T=absolute  temperature  (in  Kelvin) 

Xl=fraction  of  non-helium  gas 

X2=helium  fraction 

I=counting  variable  for  the  Do  loop 

The  output  parameters  are  Prandtl  number,  helium 
fraction,  thermal  conductivity,  and  dynamic  viscosity.  They 
are   represented   by   PR( I ) . HE( I ) , LAHIX( I ) ,   and   ETAHIX(I) 
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respectively.  The  remaining  parameters  are  all  merely 
mathematically  defined  intermediate  steps  in  the 
calculations  of  the  output  parameters  and  have  no  physical 
meaning.  The  terms  used  for  viscosity  are  XETA,  YETA,  and 
ZETA,  and  the  terms  used  for  thermal  conductivity  are  Ul, 
U2,  UY,  UZ,  XLAH,  YLAH,  and  ZLAH. 

This  section  will  now  provide  the  most  rececnt  version 
of  the  program.  Since  the  explanation  of  the  input  is 
fairly  lengthy  it  will  be  covered  first,  rather  than 
disrupt  the  flow  of  the  program.  This  will  include  a  list 
of  all  required  input  parameters,  their  sources  where 
appropriate,  and  a  sample  set  of  inputs  that  were  used. 
There  are  thirteen  input  quantities,  A12,  B12,  0MEG12,  Ml, 
H2,  CPl,  CP2,  SIGMAl,  SIGMA2,  SIGH12,  OMEGA  1 ,  0HEGA2,  AND 
T.  After  selecting  Ml,  M2,  and  T  you  have  specified  your 
gas  combination  and  temperature  of  interest.  Values  for  CP 
may  be  obtained  from  various  sources,  such  as  the  CRC 
tables,  thermodynamics  textbooks,  or  standard  tables  of 
thermodynamic  properties.  The  values  for  SIGMA  as  well  as 
i/k  are  obtained  from  table  I-A  of  Hirschf e Ider .  Using  i/k 
for  each  gas  and  tl2=(tl*^2)  you  calculate  Tl*,  T2*,  and 
T12*  (where  T*=kT/0.  Now  using  the  appropriate  T*  you  go 
to  table  I-M  of  Hirschfelder  for  the  OMEGA  values  and  table 
I-N  for  A12  and  B12.  This  provides  a  complete  set  of  input 
quantities.  The  following  table  is  a  sample  set  of  inputs 
for  two  gas  combinations  at  two  different  temperatures: 
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GAS  COMBINATION  AND  TEHPEKATURE 


QUANTITY 

OF  INTEREST 

Ml 

M2 

SIGMAl 

SIGMA2 

(  i/K)  1 

(«/K)2 

OMEGA  1 

0HEGA2 

CPl 

CP2 

( €/K)12 

T12* 

A12 

B12 

0MEG12 

He-Xe 


200 


He-Xe 


300 


He-Ar 


200 


He-Ar 


300 


131.3 

131.3 

39.948 

39.948 

4.0026 

4.0026 

4.0026 

4.0026 

4.055 

4.055 

3.418 

3,418 

2.556 

2.556 

2.556 

2.556 

229 

229 

124 

124 

10.22 

10.22 

10.22 

10.22 

1.700 

1.394 

1.275 

1.  104 

0.7467 

0.7485 

0.7467 

0.7485 

0.0382 

0.0382 

0. 12477 

0. 12477 

1.2512 

1.2512 

1.2512 

1.2512 

48.38 

48.38 

35.60 

35.60 

4.  134 

6.201 

5.618 

8.427 

1.098 

1. 1034 

1.  102 

1.  108 

1.093 

1 . 0904 

1.091 

1.0905 

0.9632 

0.8916 

0.9080 

0,8470 

PROGRAM  PRANTL 


C   First   the  array  variables  are  dimensioned  to  allow 
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C   calculating  and  plotting  the  desired  quantities.  Then 
C   the  real  quantities  are  defined  since  some  begin  with 
C   letters  that  would  make  them  integers  by  default, 
C 
C 

DIMENSION  PR( 150), HE ( 150),LAHIX( 150,ETAHIX( 150) 

REAL  H 1 , H2 , OMEGA  1 . 0MEGA2 , OMEG 12 , LAM  1 , LAM2 , LAM  12 , PR , HE 

REAL  LAMIX,ETAMIX 
C 
C 

C   Now  comes  the  input  section  where  the  counting  variable 
C   is  initialized  to  zero  and  all  required  data  is  input. 
C 
C 

A  12= 1,1080 

B 12= 1.0905 

OMEG 12=0. 8470 

1=0 

M 1=39. 948 

M2=4.0028 

CP1=0.  12477 

CP2=1.2512 

SIGMA1=3.418 

SIGMA2=2.556 

SIGM12=0.5*(SIGMA1+SIGHA2) 

OMEGA  1=1.  104 
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OMEGA2=0.7485 


T=300 


C   Now  the  output  files  are  created  and  opened  and  headings 

C   are  placed  where  appropriate. 

C 


0PEN(67,FILE='PRANTL   OUTPUT   AIM 
0PEN(68,FILE='PRANTL    RESULTS    AIM 

WRITE(68,*) *THIS  OUTPUT  FILE  IS  FOR  A  MIXTURE  OF 
*GASES  WITH' 
WRITE(88, 1)M1,H2 

1  FORMAT  ( 'ATOMIC  WEIGHTS  OF  MF10.7,'  AND',  F10.  7,  \  AND 

*  A  '  ) 
WRITE(68,2)T 

2  FORMAT ( 'TEMPERATURE  OF  ',F11.7,'  DEGREES  KELVIN.') 
WRITE(68,*) 

WRITE(68,*)'    PRANTL         HE         THERMAL 

*  DYNAMIC 


WRITE(68,*)'         NUMBER 

*  VISCOSITY' 
WRITE(68,*) 
WRITE(68,*) ' 

*  GM/CM-SEC 
WRITE(68,*) 


FRACTION    CONDUCTIVITY 


CAL/CM-SEC-K 
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C   Now  the  Do  loop  is  opened  to  perforin   the  required 

C   calculations.  The  equations  are  from  chapter  eight  of 

C   Hirschfelder  and  the  equation  numbers  are  as  follows: 

C   ETAl  and  ETA2  equation  8,2-18 

C   ETA12  equation  8.2-21 

C   LAMl  and  LAM2  equation  8.2-31 

C   LAM12  equation  8,2-35 

C   U1,U2,UY,UZ,XLAM,YLAH,  and  ZLAM  equation  8.2-36 

C   XETA,YETA,  and  ZETA  EQUATION  8,2-22 

C   The  equation  for  Prandtl  number  is  ETAHIX*CPMIX/LAmX 

C 

C 

DO  10  X1=0.0, 1.005,0.01 
1  =  1  +  1 

ETA  1=2. 6693E-5*SQRT ( H 1*T ) / (SIGMA 1*SIGHA 1*0HEGA 1 ) 
ETA2=2. 8693E-5*SQKT ( H2*T ) / ( SIGHA2*SIGHA2*0MEGA2 ) 
ETA  12=2, 8893E-5*SQRT(2*M1*H2*T/(H1+M2) )/(SIGH12* 
*SIGH12*0MEG12) 
LAM  1=1. 989 1E-4*SQRT ( T/M 1 ) / ( S I GMA 1*S I GHA 1*0HEGA 1 ) 
LAH2= 1 . 989 1E-4*SQRT ( T/M2 ) / ( SIGHA2*SIGMA2*0MEGA2 ) 
LAM  12= 1 . 989 1E-4*SQRT ( T* ( H 1+M2 ) / ( 2*H 1 *M2 ) ) / ( ( S I GM 12 
*:«c*2)*0MEG12) 

Ul=(4*A12/15)-(  (  (B12/5)  +  (  1/12)  ))lc(Hl/H2)  )  +  (  (  (H1-M2) 
***2)/(2*Hl*M2) ) 
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U2=(4*A12/15)-( ( (B12/5)+( 1/12) )*(M2/H1) )+( ( (H2-M1) 
***2)/(2*Hl*H2) ) 

UY=( (4*A12/15)*( ( (H1+M2)**2)/(4*M1*H2) )*(LAH12*LAH 
*12/(LAM1*LAH2) ) )-( (B12/5)+( 1/12) )-( ( ( 12*B12/32*A12 
*)-(25/(32*A12) ) )*( ( (H1-H2)**2)/(H1*H2) ) ) 

UZ=( (4*A12/15)*( ( (H1+M2)**@)/(4*H1*M2) )*( (LAH12/ 
*LAm)  +  (LAM12/LAH2)  )-l)  )-(B12/5)-(  1/12) 

X2=1-X1 

XETA= ( X  1*X 1/ETA 1 )  +  ( 2*X 1*X2/ETA 12 )  +  ( X2*X2/ETA2 ) 

YETA=  ( 3*A  12/5 )  *  (  (  X  1*X  1*H  1/  (  H2*ETA  1 )  )  +  ( X2*X2>tcH2/  (Ml* 
*ETA2)  )  +  (2*Xl5tcX2*(  (  M  1+H2  )**2  )*(  ETA  12*ETA  12  )  /  (  4*H  1*H2* 
*ETA1*ETA2*ETA12) ) ) 

ZETA= ( 3*A  12/5 ) * ( ( X 1*X 1*H 1/H2 )  +  ( X2*X2*H2/H 1 )  +  ( ( 2*X  1* 
*X2 )*((((( M 1+H2 ) **2 ) / ( 4*H 1*M2 ) ) * ( ( ETA 12/ETA 1 )  +  ( ETA  12 
*/ETA2) ) )-l) ) ) 

XLAH= ( X 1*X 1/LAM 1 ) + ( 2*X 1*X2/LAH 12 ) + ( X2*X2/LAM2 ) 

YLAH=  (  X  1*X  1*U  1/LAM  1 )  +  (  2*X  1*X2*UY/LAM  12 )  +  ( X2>»cX2*U2/ 
*LAH2) 

ZLAH= ( X  1*X 1*U 1 )  +  ( 2*X 1*X2*UZ )  +  ( X2*X2*U2 ) 

PR( !)=( ( 1+ZETA)*(XLAH+YLAH)*( (X 1*H 1*CP 1 ) + ( X2*H2*CP2) 
*))/(( XETA+YETA ) * ( 1+ZLAH ) * ( ( X 1*H 1 ) + ( X2*H2 ) ) ) 

HE(I)=X2 

LAMIXC I )=( 1+ZLAH )/(XLAH+YLAM) 

ETAHIX(I)=( 1+ZETA) /(XETA+YETA) 
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C   Now  we  will  write  to  the  output  files  the  data  that  has 

C   been  calculated. 

C 

C 

WRITE(67,5)PR(I ) , HE ( I) 

WRITE(68,6)PR( I),HE( I),LAMIX( I),ETAHIX( I) 

5  FORMAT(F10.7,F10,2) 

6  FORMAT  (  F 10.  7,  F  12.  2,  E 15.  7,  E 15.  7) 
10  CONTINUE 

STOP 
END 
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APPENDIX  B,   PLOTTING  PROGRAM 


This  appendix  provides  a  program,  written  by  the 
author,  that  will  draw  curves  of  Prandtl  number  versus 
helium  fraction  from  the  output  file,  "PRANTL  OUTPUT  Al", 
created  by  running  the  program  "PRANDTL  FORTRAN  Al"  that 
was  presented  in  Appendix  A.  This  program  must  be  compiled 
in  FORTVS  and  is  specifically  designed  to  use  DISSPLA  on 
the  Naval  Postgraduate  School  mainframe  computer.  Various 
DISSPLA  subroutines  are  called  by  the  program  and  each 
will  be  explained  with  comment  statements.  DISSPLA 
requires  the  data  that  is  to  be  plotted  to  be  in  an  array 
format.  Each  of  the  arrays  used  here  are  set  to  a  maximum 
of  150  points.  If  modifications  are  made  to  "PRANDTL 
FORTRAN  Al"  such  that  more  than  150  points  are  calculated 
the  array  sizes  must  be  redefined.  Two  character  strings 
are  used  in  this  program.  They  are  NEW1$  and  NEW2$,  both 
of  which  are  used  to  rename  files.  The  first  changes  the 
input  filename  into  the  format  required  for  DISSPLA  and 
the  second  changes  the  file  back  to  its  original  name  for 
storage.  The  only  parameters  used  are  PR(I)  and  HE(I) 
which  represent  Prandtl  number  and  helium  fraction, 
respectively.  Both  of  these  quantities  are  real  numbers 
and  are  defined  as  such  in  the  program.  This  completes  the 
introduction   and   explains   the   first   three  lines  of  the 
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program.   The   remainder  of  this  appendix  will  be  a  listing 
of   the  program  with  explanatory  comments.  Latest  update  of 
\   program,  March  1988. 

Program  GKFPRN 

character*80  NEW1$,NEW2$ 

dimension  PR( 150) , HE( 150) 

real  PR, HE 
c 

c  The  next  three  lines  rename  the  input  file  to  the 
c  format  necessary  for  DISSPLA  and  open  it  for  use. 
c 

NEW  1$=' rename  V/*  PRANTL  OUTPUT  AlV/'  file 
*XYZ1  Al' 

call  excms(NEWl$) 

open(87,f ile='XYZl» ) 
c 

c   In  the  next  seven  lines  a  Do  loop  is  used  to  read 
c   all  of  the  data  from  the  input  file  into  the 
c   specified  arrays.  Then  the  input  file  is  closed  and 
c   renamed  it  to  its  original  name.  The  format  line 
c   ensures  the  data  is  read  from  the  file  in  the  same 
c   format  as  it  was  written  into  the  file. 
c 

1  format(f 10.7, f 10.2) 

do  10  1=1, 100 
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read(67, 1)PK(I),HE(I) 
10     continue 
20  close(67) 

New2$=' rename  V/'  file  XYZl  AlV/*  PRANTL 
*OUTPUT  Al' 

call  excms(NEW2$) 
c 

c   The  remainder  of  the  program  calls  DISSPLA  subroutines 
c   to  create  and  plot  the  desired  curves.  Call  tek618 
c   sets  up  the  required  machine  for  the  plotting  and 
c   links  you  to  the  plotter.  You  must  be  on  a  tekSlS 
c   terminal  to  use  this  program.  Call  blowup  is  used  to 
c   magnify  the  plot  to  near  full  page  size,  for  clarity. 
c   Call  page  sets  up  the  page  size,  in  inches.  Call  area2d 
c   sets  up  a  two  dimensional  plot  area  of  the  specified 
c   size,  in  inches.  Call  xname  and  yname  allow  you  to 
c   label  the  axes.  Call  headin  allows  naming  the  plot, 
c   Call  graf  sets  up  the  physical  origin  and  the  maxima 
c   for  the  axes.  The  'scale*  portion  of  this  call  allows 
c   DISSPLA  to  select  appropriate  increments  for  the  axes 
c   based  on  the  minima  and  maxima  specified.  Call  thkfrm 
c   and  call  frame  place  frames  around  the  inner  plot  area 
c   and  the  entire  plot,  of  a  specified  thickness.  Call 
c   curve  tells  the  plotter  to  plot  the  specified  arrays  of 
c   data  points.  The  arrays  must  be  called  in  the  proper 
c   order  (y-axis,  x-axis)  to  get  the  correct  curve.  Call 
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c  endpl  ends  the  specific  plot  and  sets  up  a  new  page  for 
c  another  plot.  Call  donepl  tells  the  device  you  are  done 
c   plotting,  . 


cal 
cal 
cal 
cal 
cal 
cal 
cal 
*100 
cal 
cal 
cal 
cal 
cal 
cal 
stop 
end 


tek618 
blowup( 1.5) 
page( 11,8.5) 
area2d(8. ,8. ) 

xname ( 'hel ium  f raction$* , 100) 
yname ( 'Prandt 1  number^' , 100) 

headin ( *Prandt 1  number  vs.  helium  fraction^', 
1.  ,  1) 

graf (0.0, 'scale' , 1.0,0.4, 'scale' ,0,7) 
thkfrm( .02) 
frame 

curve(HE,PR, 100,0) 
endpl (0) 
donepl 
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APPENDIX  C.   RAW  DATA 

This  Appendix  provides  a  listing  of  our  spreadsheet, 
discussed  in  Chapter  V,  including  all  recorded  data  and 
calculated  values  that  we  have  used.  Many  of  the  values  are 
not  related  to  our  coefficient  of  performance  evaluation, 
but  were  calculated  for  informational  purposes. 

The  Appendix  is  broken  down  into  three  sets  of  three 
pages  each,  in  order  to  fit  the  24  column  by  122  row  spread- 
sheet within  the  desired  margins.  Each  set  of  three  pages 
covers  the  entire  24  columns  for  the  specified  data  sets, 
with  nine  columns  per  page.  We  will  discuss  each  of  the 
columns,  omitting  the  first  which  specifies  the  data  set. 

Columns  1,2,3,4,6,8,9,10  and  11  are  all  raw  data  and  are 

self   explanatory.  Column  5  is  the  measured  mean  pressure  of 

column   4,  after  conversion  from  milli-volts  to  psia  (at  1.5 

psia   per   mv).   Column   7   converts   the   measured   dynamic 

pressure   signal   of   column  6  to  percent  pressure  amplitude 

(p^/p  =1586xC6/C5) .   Column   12  is  an  acceleration  term  that 
H)      m 

is  calculated  for  use  in  the  work  equation  as  a  correction 
for  accelerometer  non-linearities  ( acc=0. 2848xC3xC8/ 170. 5) . 
Column  13  is  the  absolute  temperature  ratio 
(Tp/T„=(C2+273.  15)/(Cl+273.  15) ) .  Column  14  is  the  applied 
heater  load  in  watts  (QHEAT=C9xC 10x0. 00 1  or 
QHEAT=C10xC10xCll/10^,   as   appropriate).   Column   15  is  the 
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work  in  watts  ( W=59810x( 1-7. 4x  10  ^xC 12)xC6x0. 003xC8/C 1839. 8- 
0.05(03-550)],  as  discussed  in  Chapter  V.  Column  16  is  the 
coefficient  of  performance,  neglecting  the  heat  leak  from 
the  cold  end  of  the  refrigerator  (COP=C 14/C 15) .  Column  17 
is  the  coefficient  of  performance  relative  to  Carnot, 
neglecting  the  heat  leak  (COPR=C 16( ( 1/C 13) -1 ) ) .  Column  18 
is  the  heat  leak  in  watts,  as  discussed  in  Chapters  IV  and 
V  (HEATLK=0.0064(C1-C2) ) .  Column  19  is  the  coefficient  of 
performance,  taking  into  account  the  heat  leak 
(C0P2=(C14+C18)/C15) .  Column  20  is  the  coefficient  of 
performance  relative  to  Carnot,  taking  into  account  the 
heat  leak  (C0PR2=C 19( ( 1/C 13)-1) ) .  Column  21  is  the  total 
heat  load  in  watts  (Qtot=C  14+C 18) .  Column  22  is  the  overall 
temperature  difference  (DELTA  T=C1-C2).  Column  23  is  an 
approximate  viscous  penetration  depth  for  the  given 
temperature,  pressure,  frequency  and  gas  type  and  an 
approximate  value  for  dynamic  viscosity  (VISC 
PEN=SQRT( (Cl+273. 15) (0.005014692852)/(C3xC5) ),  where  the 
coefficient  takes  into  account  all  constants  and  the  value 
used  for  dynamic  viscosity).  Column  24  is  an  approximate 
thermal  penetration  depth  calculated  from  &  =i  o-  (THERM 
PEN=1.219820636xC23) .  The  equations  for  viscous  and  thermal 
penetration  depths  use  a  single  value  for  and  <r  within 
each  data  set,  respectively,  instead  of  a  value  calculated 
for   each   data  point.  This  is  done  for  ease  of  calculation. 
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The   resulting   values   are  therefore  not  precise  but  merely 
reasonable  approximations. 

The  spreadsheet  of  raw  data  and  calculated  values  is 
provided  on  the  next  nine  pages,  and  does  not  follow  on 
this  page  for  formatting  reasons. 
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TH 

TC 

FREQ 

PMEAN 

PMEAN 

VPRESS 

Po/Pb 

VLOCKIN 

VHEAT 

deg  C 

deg  C 

Hz 

nVrns 

psia 

Vruis 

7. 

(iiVrms 

V 

25.1 

-58.6 

522.03 

100. 17 

150.26 

0.28391 

2.997 

367.3 

0.000 

22.4 

-57.2 

523.20 

99.98 

149.97 

0. 18870 

1.996 

241.6 

0.000 

old 

22.6 

-44.7 

537.10 

100.00 

150.00 

0. 18890 

1.997 

262.1 

6.870 

stack 

22.9 

-34.0 

548.70 

99.97 

149.96 

0.18860 

1.995 

277.0 

9.714 

pure  He 

23.2 

-24.7 

558.70 

99.98 

149.97 

0.18940 

2.003 

289.1 

11.897 

short 

20.3 

-13.5 

570.40 

101.06 

151.59 

0. 18870 

1.974 

307.0 

13.780 

resonator 

23.1 

-13.4 

570.40 

100.40 

150.60 

0.18900 

1.990 

309.0 

13.780 

26.7 

-43.0 

539.60 

99.93 

149.90 

0.28280 

2.992 

391.5 

13.780 

22.8 

-62.2 

518.10 

99.96 

149.94 

0. 18960 

2.006 

228.2 

0.000 

24.2 

-71.3 

508.40 

100.05 

150.08 

0.28390 

3.000 

326.3 

0.000 

22.6 

-67.5 

507.94 

99.54 

149.31 

0.18902 

2.008 

233.1 

25.9 

-68.6 

506.90 

100.20 

150.30 

0.28332 

2.990 

349.8 

-44.7 

533.70 

104.17 

156.26 

a. 18873 

1.916 

256.1 

-34.0 

544.90 

102.50 

153.75 

0.18903 

1.950 

230.2 

new 

23.6 

-24.9 

554.49 

103.20 

154.80 

0. 18947 

1.941 

293.1 

stack 

24.1 

-13.5 

562.52 

102.00 

153.00 

0. 18939 

1.963 

319.8 

pure  He 

26.9 

-43.1 

535.55 

100.10 

150. 15 

0.28340 

2.993 

416.0 

short 

25.3 

-73.5 

501.57 

100.04 

150.06 

0.28367 

2.998 

337.2 

resonator 

-68.6 

507.10 

99.97 

149.96 

0.18900 

1.999 

230.1 

20.7 

-34.0 

544.50 

100.09 

150.14 

0.09465 

1.000 

147.7 

23.7 

-24.5 

554.62 

100.05 

150.08 

0.18886 

1.996 

309.0 

23.7 

-13.4 

565.44 

100.11 

150.17 

0. 18922 

1.998 

331.7 

23.0 

-44.6 

533.81 

100.07 

150.11 

0.18899 

1.997 

272.5 

23.7 

-68.9 

235.50 

98.67 

148.01 

0.28390 

3.042 

169.3 

0.000 

24.1 

-60.5 

239.86 

98.13 

147.20 

0.28310 

3.050 

181.4 

6.871 

old 

24.2 

-51.9 

244.21 

97.84 

146.76 

0.28310 

3.059 

194.1 

9.706 

stack 

24.3 

-42.7 

248.82 

97.54 

146.31 

0.28380 

3.076 

208.6 

11.904 

He-Xe 

22.4 

-38.4 

250.88 

97.94 

146.91 

0. 18870 

2.037 

142.2 

8.416 

short 

-70.1 

234.76 

98.84 

148.26 

0.28330 

3.031 

167.9 

0.000 

resonator 

21.0 

-62.5 

238.85 

99.14 

148.71 

0.19000 

2.026 

118.6 

0.000 

22.2 

-17.5 

260.29 

98.60 

147.90 

0.1B950 

2.032 

159.3 

11.893 

22.9 

-77.6 

229.24 

100.99 

151.49 

0.28380 

2.971 

160.3 

24.1 

-60.5 

238.23 

101.12 

151.68 

0.2B340 

2.963 

181.7 

neM 

23.9 

-51.9 

242.48 

101.03 

151.55 

0.28291 

2.961 

193.5 

stack 

24.5 

-42.8 

246.82 

101.73 

152.60 

0.28291 

2.940 

202.8 

He-Xe 

22.1 

-TO  T 

--■Ua  >J 

248.93 

101.63 

152.45 

0.18910 

1.967 

141.5 

short 

on  -r 

-50.0 

243.33 

101.71 

152.57 

0.13870 

1.962 

130.0 

resonator 

21.9 

-67.5 

234.57 

102.03 

153.05 

0.18910 

1.960 

113.2 

23.5 

-77.8 

229,07 

100. 16 

150.24 

0.28340 

2.992 

161.6 

23.3 

-80.3 

227.13 

89.60 

134.40 

0.28421 

3.354 

183.8 
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IHEAT 

RHEAT 

ACC 

TC/TH 

QHEAT 

WORK 

COP 

COPR 

HEATLK 

IDA 

ohms 

cg5E-3 

watt 

watt 

Q=heater 

cl7/c3rn 

watt 

0.00 

320.28 

0.7194 

0.0000 

9.9215 

0.0000 

0.0000 

0.5692 

0.00 

211.14 

0.7307 

0.0000 

4.3736 

0.0000 

0.0000 

3.5413 

old 

72.81 

235. 15 

0.7725 

0.5002 

4.7429 

0.1055 

0.0311 

0.4576 

stack 

102.82 

253.88 

0.8078 

0.9988 

4.9991 

0.1998 

0.0475 

0.3369 

pure  He 

125.82 

269.30 

0.3384 

1.4969 

5.2343 

0.2860 

0.0551 

0.3257 

short 

145.60 

292.51 

0.8848 

2.0064 

5.5306 

0.3628 

0.0472 

0.2298 

resonator 

145.60 

294.41 

0.8768 

2.0064 

5.5747 

0.3599 

0.0506 

0.2432 

146.00 

352.87 

0.7676 

2.0119 

10.5129 

0.1914 

0.05S0 

0.4740 

0.00 

197.49 

0.7123 

0.0000 

4. 1544 

0.0000 

0.0000 

0.5780 

0.00 

277.10 

0.6783 

0.0000 

8.8393 

0.0000 

0.0000 

0.6494 

0.00 

94.00 

197.77 

0.6954 

0.0000 

4.2636 

0.0000 

0.9000 

0.5676 

0.00 

94.00 

296. 18 

0.6840 

0.0000 

9.5190 

0.0000 

0.0000 

0.5954 

39.50 

94.40 

223.31 

0.7709 

0.7562 

4.6697 

0.1619 

0.0481 

0.4278 

116.90 

94.48 

255.04 

0.3062 

1.2911 

5.1085 

0.2527 

0.0603 

new 

133.20 

94.55 

271.47 

0.8366 

1.6775 

5.3509 

0.3135 

0.0612 

0,3056 

stack 

158.30 

94.64 

300.49 

0.S735 

2.3716 

5.3244 

0.4072 

0.0590 

0.2369 

pure  He 

158.00 

94.40 

372.36 

0.7667 

2.3566 

11.2885 

0.2088 

0.0635 

0.4410 

short 

0.00 

94.00 

282.51 

0.6690 

0.0000 

9.1956 

0.0000 

0.0000 

0.d224 

resonator 

0.00 

94.00 

194.91 

0.6916 

0.0000 

4.2091 

0.0000 

0.0000 

0.5746 

0.00 

94.48 

134.34 

0.8139 

0.0000 

1.3606 

0.0000 

0.0000 

0.3446 

138.00 

94.55 

236.27 

0.8376 

1.8006 

5.6163 

0.3206 

0.0621 

0.3037 

156.80 

94.64 

313.29 

0.8750 

2.3268 

6.0303 

0.3359 

0.0551 

0.2337 

99.50 

94.40 

242.98 

0.7717 

0.9346 

4.9701 

0. 1330 

0.0556 

0.4259 

0.00 

72.97 

old 

102.94 

stack 

126.09 

He-Xe 

89.11 

short 

0.00 

resonator 

0.00 

125.77 

66.60 

0.6881 

0.0000 

4.6249 

0. 0000 

0. 0000 

0.5926 

72.68 

0.7154 

0.5014 

4.9399 

0.1015 

0.0404 

0.5414 

79,18 

0.7441 

0.9992 

J«  jiDC'u 

0.1891 

0.0650 

0.4870 

86.70 

0.7748 

1.5010 

5.6900 

0.2638 

0.0767 

0.4283 

59.59 

0.7943 

0.7499 

2.5844 

0.2902 

0.0752 

0.3891 

65.34 

0.6852 

0.0000 

4.5//2 

0.0000 

0.0000 

0.5971 

47.32 

0.7161 

0.0000 

2.1716 

0.0000 

0.0000 

0.5344 

69.48 

0.8656 

1.4964 

2.9152 

0.5133 

0.0797 

0.2541 

0.00 

94.00 

61.38 

0.6605 

0.0000 

4.4137 

0.0000 

0.0000 

0.ul31 

95.70 

94.16 

72.32 

0.7154 

0.3624 

4.9930 

0.1727 

0.0687 

a. 5161 

new 

119.30 

94.29 

/  O  •  0  ; 

0.7448 

1.3420 

5.3063 

0.2529 

3.0866 

3.4624 

stack 

141.10 

94.41 

33.61 

0.7/39 

1  010  L 
1 . u/  /u 

5.5593 

0.3381 

0.0983 

0.4105 

He-Xe 

96.80 

94.45 

5S.84 

0.7954 

0.3850 

2.5979 

0.3407 

0.0876 

0.3684 

short 

73.50 

94.32 

52.34 

0.7553 

0.5095 

2.3324 

0.2139 

0.0693 

0.4410 

resonator 

B.S0 

94.00 

44.35 

0.6970 

0. 0000 

2.0797 

0. 0000 

0.0000 

0.5453 

0.00 

94.00 

61.83 

0.6585 

0.0000 

4.4430 

0.0000 

0.0000 

0.6179 

3.00 

94.00 

69.73 

0. 6505 

0.0000 

J. 0646 

0. 0000 

0.0000 

0.6320 
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C0F2 

C0PR2 

Qtot 

DELTA  T 

Vise  PEN  THERM  PEN 

w  cl9 

c20/carn 

watt 

deg 

cm 

cm 

0.0574 

0.0224 

0.5692 

83.70 

0.003307 

0.010133 

0.1233 

0.0456 

0.5413 

79.60 

0.008268 

0.010086 

old 

0.2020 

0.0595 

0.9578 

67.30 

0.008162 

0.009957 

stack 

0.2772 

0.0659 

1.3357 

56.90 

0.008081 

3.309857 

pure  He 

0.3482 

0.0671 

1.3226 

47.90 

0.008012 

0.009773 

short 

0.4043 

0.0526 

2.2362 

33.80 

0.007848 

0.009573 

resonator 

0. 4044 

0.0568 

2.2546 

36.50 

3.007911 

0.309651 

0.2365 

0.0716 

2.4858 

69.70 

0.008203 

0.310006 

0.1391 

0.0561 

0.5780 

85.00 

0.008315 

0.010143 

0.0735 

0.0348 

0.6494 

95.50 

0.036410 

3.010259 

0.1331 

0.0583 

0.5676 

90.10 

0.008413 

3.010262 

0.3625 

0.0289 

0.5954 

94.50 

0.008440 

3.010296 

0.2535 

0.0754 

1.1339 

67.90 

0.008031 

0.009796 

0.3237 

0.0778 

1.6534 

57.50 

3.008017 

3.009779 

new 

0.3706 

0.0724 

1.9831 

48.50 

3.007921 

9.309663 

stack 

3.4473 

0.3649 

2.6085 

37.60 

3.007917 

0.009653 

pure  Ke 

0.247S 

0.0754 

2.7976 

73.00 

0.303229 

0.313038 

short 

0.0677 

0.0335 

0.6224 

93.80 

0.008483 

0.010343 

resonator 

0. 1365 

3.0609 

3.5746 

91.20 

0.008402 

0.010249 

0.2533 

0.0579 

0.3446 

54.70 

3.003077 

3.009853 

0.3746 

0.0726 

2.1343 

48.20 

0.008045 

3.009314 

0.4246 

0.3606 

2.5606 

37.10 

0.007966 

0.009717 

0.2737 

0.0810 

1.3605 

67.60 

0.008190 

0.009991 

0.1231 

0.0531 

0.5926 

92.60 

0.006535 

0.012596 

0.2111 

0.0840 

1.0423 

84.60 

0.006498 

0.012524 

old 

0.2813 

0.0967  .. 

1.4362 

76.10 

0.006450 

0.012432 

stack 

0.3392 

3.0986 

1.9293 

67.00 

0.006401 

3.012338 

He-Xe 

0.4407 

3.1141 

1.1391 

60.30 

3.306341 

0.012223 

short 

3.1305 

0.0599 

3.5971 

93.33 

0.006534 

0.012595 

resonator 

0.2461 

0.3975 

3.5344 

33.50 

0.006444 

0.012421 

0.6005 

3.0932 

1.7505 

39.70 

0.006203 

0.011955 

3. 1339 

3.3714 

0.G131 

103.50 

3.006533 

0.312603 

0.2761 

e.ie98 

i.3784 

84.60 

3.306422 

3.012378 

new 

3.3400 

0.1165 

J. 3044 

75.80 

3.006-367 

0.012272 

stack 

3.4119 

e.1203 

2.2902 

67.30 

0.006295 

0.012134 

He-Xe 

3. 4825 

0.1241 

1     -iC-Tir 
1 .  .<.JOuJ 

60.40 

0.006246 

0.012040 

short 

0.3990 

0.1293 

0.9506 

72.33 

0.006317 

0.012177 

resonator 

0.2622 

0.1140 

0.5453 

89.40 

0.006420 

3.012374 

0. 1391 

3.3721 

0.6179 

101.33 

0.336575 

0.012672 

0. 1248 

0.3670 

0.6-320 

103.60 

0.006978 

0.013451 
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TH 

TC 

FREQ 

PMEAN 

PMEAN 

VPRESS 

Po/PiD 

VLOCKIN   VHEAT 

deg  C 

deg  C 

Hz 

raVrms 

psia 

Vrms 

7. 

aiVrms     V 

23.7 

-72.5 

291.33 

100.34 

151.26 

0.28311 

2.968 

190.1 

new 

24.4 

-60.4 

299.42 

100.49 

150.74 

0.28301 

2.973 

208.1 

stack 

24.2 

-51.8 

305.07 

100.83 

151.32 

0.28321 

2.968 

220.3 

He-Ar 

25.0 

-42.9 

310.85 

101.84 

152.76 

0.28303 

2.939 

229.8 

(227.) 

22.3 

-38.2 

313.60 

100.82 

151.23 

0. 18950 

1.987 

160.6 

short 

22.5 

-50.1 

306.06 

100.25 

150.38 

0. 18378 

1.991 

149.3 

resonator 

22.2 

-65.0 

296.28 

100.25 

150.38 

0.18921 

1.996 

135.1 

24.0 

-72.2 

291.54 

100.25 

150.33 

0.28333 

2.988 

191.7 

22.9 

-60.3 

470.43 

100.40 

150.60 

0.18916 

1.992 

212.8 

23.9 

-44.3 

485.56 

99.93 

149.97 

0. 18902 

1.999 

249.2 

new 

23.9 

-34.0 

495.98 

99.95 

149.93 

0. 18920 

2.001 

276.7 

stack 

24.6 

-25.1 

504. 18 

99.30 

149.70 

0. 18902 

2.003 

298.7 

pure  He 

24.9 

-13.7 

514.63 

100.38 

150.57 

0.18934 

1.994 

■ji.b .  0' 

long 

22.9 

-61.2 

469.90 

100.10 

150.15 

0. 18932 

2.000 

212.2 

resonator 

27.5 

-43.1 

488. 19 

99.92 

149.88 

0.28343 

2.999 

336.0 

23.3 

-32.3 

498.54 

100.00 

150.00 

0.23365 

2.999 

425.6 

28.8 

-32.3 

498.52 

99.93 

149.90 

0.28359 

3.001 

426.2 

26.1 

-62.0 

468.88 

99.98 

149.97 

0.23335 

2.997 

316.7 

24.4 

-64.0 

284.28 

99.94 

149.91 

0.28315 

2.996 

183.5 

new 

24.7 

-60.3 

286.60 

99.32 

149.73 

0.28300 

2.998 

191.7 

stack  ■ 

25.2 

-51.8 

291.71 

99.70 

149.55 

0.28310 

3.002 

208.1 

He-Ar 

25.3 

-43.0 

296.91 

99.64 

149.46 

0.28268 

3.000 

226.7 

(19/:) 

22.9 

-38.4 

299.55 

100.09 

150.14 

0. 18942 

2.001 

159.6 

long 

22.9 

-50.1 

292.58 

99.99 

149.99 

0.18923 

2.001 

142.0 

resonator 

22.1 

-60.9 

286.15 

100.02 

150.03 

0.13982 

2.007 

126.1 

24.4 

-64.9 

203.80 

99.98 

149.97 

0.28274 

2.990 

182.1 

21.8 

-63.1 

207.92 

99.78 

149.67 

0.18940 

2.007 

94.1 

new 

22.6 

-50.1 

213.39 

99.94 

149.91 

0. 18902 

2.000 

108.7 

stack 

22.4 

-38.2 

218.27 

99.91 

149.87 

0. 18948 

2.005 

124.2 

He-Xe 

23.2 

-26.2 

222.90 

99.83 

149.75 

0.18916 

2.003 

133.6 

long 

23.0 

-14.0 

227.38 

99.91 

149.87 

0.18957 

2.006 

154.2 

resonator 

23.8 

-2.3 

231.63 

100.01 

150.02 

0. 13923 

2.001 

167.3 

21.7 

-63.0 

207.98 

99.85 

149.78 

0.18941 

2.006 

94.3 

23.0 

-41.0 

490.81 

100.01 

150.02 

0. 18892 

1.997 

193.3 

23.7 

-34.1 

497.57 

99.34 

149.76 

0.13838 

1.995 

208.3 

dual 

23.5 

-24.9 

506.39 

99.33 

149.75 

0.18973 

2.010 

226.1 

stack 

24.1 

-13.1 

517.38 

100.56 

150.84 

0. 13866 

1.984 

241.1 

pure  He 

24.2 

-2.7 

526.79 

100.35 

150.53 

0. 18900 

1.991 

260.2 

long 

27.9 

-13.3 

517.66 

100.27 

150.41 

0.28320 

2.936 

372.1 

resonator 

27.2 

-24.1 

507.60 

100.16 

150.24 

0.28346 

2.992 

OHDa  O 

27.1 

498.35 

100.13 

150.27 

0.28369 

2.994 

322.8 

26.0 

-Oq.6 

493.36 

100.20 

150.30 

0.23325 

2.989 

309.7 

24.8 

-JJ.  J 

472. 10 

59.68 

89.52 

0.17804 

3. 154 

J  jO  ■  0 
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I  HEAT 

RHEAT 

ACC 

TC/TH 

QHEAT 

WORK 

COP 

COPR 

HEATLK 

mA 

ohms 

cgsE-3 

watt 

watt 

Q=heater 

cl7/carn 

watt 

0.00 

94.00 

92.51 

0.6759 

0.0000 

5.2182 

0.0000 

0,0000 

0.5530 

new 

80.10 

94.16 

104.03 

0.7150 

0.6041 

5.7066 

0.1059 

0.0422 

0.4918 

stack 

107.70 

94.29 

112.26 

0.7444 

1.0937 

6.0427 

0.1310 

0.0621 

0.4408 

He-Ar 

128.40 

94.41 

119.32 

0.7723 

1.5565 

6.2969 

0.2472 

0.0729 

0.3938 

(22::) 

91.20 

94.45 

34.13 

0.7952 

0.7856 

2.9544 

0.2659 

0.0685 

0.3509 

short 

67.10 

94.32 

76.53 

0.7544 

0.4247 

2.7462 

0,1546 

0.0503 

0.4211 

resonator 

0.00 

94.00 

66.36 

0.7048 

0.0000 

2.4835 

0,0000 

0.0000 

0.5058 

0.00 

94.00 

93. 35 

0. 6763 

0.0000 

5.2659 

0. 0000 

0.0000 

0.5530 

0.00 

94.16 

167.22 

0.7190 

0.0000 

3.9001 

0.0000 

0.0000 

9.5242 

123.00 

<?4.40 

202. 12 

0. 7637 

1.3594 

4.5538 

0.2985 

0.0893 

0,4323 

new 

159.30 

94.48 

229.24 

0.8051 

2.3385 

5.0523 

0.4723 

0.1145 

0.3643 

stack 

133.50 

94.55 

251.56 

0.0331 

3.1837 

5.4409 

0.5851 

0.1172 

0,3131 

pure  He 

212.60 

94.64 

230.50 

0.3705 

4.2776 

5.9424 

0.7193 

0.1071 

0,2432 

long 

0.00 

94.14 

166.56 

0.7159 

0.0000 

3.8926 

0.0000 

0.0000 

0,5293 

resonator 

190.50 

94.40 

314.77 

0.7652 

3.4258 

10.4831 

0,3266 

0. 1002 

0,4443 

237.60 

94.49 

354.42 

0.7^77 

Ja  JC'4j 

11.5415 

0,4622 

0.1172 

0,3349 

233. 10 

94.49 

354.90 

0.7977 

5.3568 

11.5549 

0. 4636 

0,1176 

0,3349 

0.00 

94.12 

248.04 

0.7056 

0.0000 

8.6416 

0,0000 

0.0000 

0,5550 

0.00 

94.08 

87.14 

0.7029 

0.0000 

5.0383 

0,0000 

0.0000 

0,5746 

new 

70.80 

94.16 

91.77 

0.7146 

0.4720 

5.2597 

0,3897 

0.0358 

0.5525 

stack 

130.60 

94.29 

101.40 

0.7419 

1.6082 

5.7084 

0.2817 

0.0980 

0.5005 

He-Ar 

168.63 

94.40 

112.43 

0.7699 

2.6834 

6.2051 

0.4325 

0.1293 

0.4472 

(197.) 

121.50 

94.44 

79.36 

0.7929 

1.3941 

2.9346 

0.4751 

0.1241 

0.3985 

long 

81.70 

94.32 

69.40 

0.7534 

0.6296 

2.6099 

0.2412 

0.0789 

0,4745 

resonator 

0.00 

94.15 

60.27 

0.7189 

0.0000 

2.3261 

0.0000 

0.0000 

0.5395 

3.00 

94.06 

86.33 

0.6999 

0.0000 

4.9933 

0.0000 

0. 0000 

0.5305 

0.00 

94.15 

32.68 

0.7122 

0.0000 

1  TTIp 

0.0000 

0.0000 

0.5349 

new 

33.40 

94.29 

38.75 

0.7542 

0.6558 

1.9959 

0.3286 

0.1071 

0.4580 

stack 

120.20 

94.42 

45.23 

0.7950 

1.3642 

2.2352 

0.5970 

0.1540 

0.3818 

He-);e 

146.40 

94.54 

51.60 

0.8333 

2.0263 

2.5450 

0.7962 

0.1593 

0.3112 

long 

171.70 

94.67 

58.57 

0.8751 

2.7910 

2.8365 

0.9840 

0. 1405 

resonator 

190.60 

94.80 

64.92 

0.9121 

3.4439 

3,0300 

1.1182 

0.1077 

0. 1644 

0.03 

94.15 

32.76 

0.7127 

0.0000 

1.7356 

0.0000 

0.0000 

0.5336 

0.80 

94.39 

162.98 

0.7339 

0.0000 

3.6356 

0. 3000 

0.0000 

0.3776 

72.70 

94.46 

173,54 

3.3053 

0.4992 

3.3053 

0.1312 

0.0317 

0.3410 

dual 

114.53 

94.56 

191.25 

0.3369 

1.2397 

4.1467 

3.2990 

0.0533 

0.2856 

stack 

148.30 

94.67 

203.36 

0.3749 

2.0321 

4.3914 

0.4741 

9.0673 

9.2195 

pure  He 

175.70 

94.80 

223.96 

0. 9095 

2.9265 

4.7417 

9.6172 

0.0614 

0.1537 

long 

207.00 

94.67 

321.75 

0.3632 

4.0565 

10.0371 

9.4021 

0.063S 

0.2431 

resonator 

154.30 

94.57 

292.77 

0.3292 

2.2662 

9. 3872 

3,2414 

0.0497 

0. 3027 

85.70 

94. 4d 

263.71 

0.7985 

0.6933 

8.7963 

0,0789 

0.0199 

0.3570 

0.00 

94.41 

255.22 

0.7341 

0,0000 

8.4337 

0,0000 

0.0000 

0.3611 

0.00 

94.23 

273.84 

0.7305 

0.0000 

6.0381 

9.0000 

0.0000 

0.4733 
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C0P2 

C0PR2 

Qtot 

DELTA  T 

Vise  PEN  THERM  PEN 

w  cl9 

c23/carn 

■^tt 

deg 

cm 

cm 

3.1069 

0.0513 

3.5580 

96.20 

0.007043 

0.010759 

new 

S.1921 

3.0765 

1.09o0 

34.30 

0.036967 

3.310644 

stack 

0.2539 

0.3372 

1.5345 

76.00 

0.006837 

0.310521 

He-Ar 

S.3097 

0.0913 

1.9503 

67.90 

0.006799 

0.010337 

i22'.) 

3.3847 

0.0990 

1.1365 

60.50 

0.006773 

0.010347 

short 

3.3080 

0.1032 

0.3457 

72,60 

0.806377 

0.010507 

resonator 

S.2B3d 

0.0853 

0.5053 

87.20 

0.006936 

3.010673 

0.1060 

3.0507 

3.5530 

96.20 

3.007364 

0.013792 

0.0525 

2.5242 

33.20 

0.003739 

0.010623 

0.3936 

0.1134 

1.7922 

63.70 

0.008604 

0.313496 

new 

0.5450 

3.1319 

^«  i'j-'j-j 

57.90 

3.003515 

0.010387 

stack 

0.6427 

0.1238 

3.4963' 

49.73 

0.008462 

0.310322 

pure  He 

0.7608 

0.1132 

4.5208 

38.60 

0.008355 

0.313192 

iung 

0.1361 

0.3540 

3.5293 

84.10 

0.008727 

0.310645 

resonator 

3.3693 

3.1133 

3.  S706 

70.60 

0.338636 

0.010534 

0.4955 

3.1257 

5.7192 

61.10 

3.038561 

3.313442 

0.4969 

0.12ul 

5.7417 

61.10 

0.008564 

3.010446 

0.0642 

3.2263 

0.5550 

38.10 

0.008739 

0.010720 

0.1140 

0.0482 

0.5746 

38.40 

0.007475 

3.311314 

new 

0.1948 

3.0773 

1.0245 

85.00 

0.007453 

0.011231 

stacK 

3.3694 

2.1087 

77.00 

0.007393 

0.011193 

He-Ar 

0.5045 

0. 1508 

3. 1306 

63.30 

0.007343 

0.011114 

(19;U 

0.6109 

3.1595 

1.7926 

61.30 

0.007253 

0.010986 

long 

0.4230 

3.1334 

1.1041 

73.00 

0.007343 

0.311122 

resonator 

0.2319 

0.0907 

0. 5395 

33.00 

0.007419 

0. 01 1229 

0.1162 

3.3498 

3.5805 

39.30 

0.307480 

3.011322 

0.3089 

0.1248 

0.5349 

84.93 

0.006394 

0.013288 

new 

0.5531 

0.1319 

1.113S 

72.70 

3.006809 

0.813124 

stack 

0.7640 

0.1971 

1.7460 

60.60 

3.006731 

0.012974 

He-,^e 

0.9185 

0.1337 

49.40 

0.006673 

0.012861 

long 

l.B6fcl 

0. 1522 

3.3241 

37.00 

0.006602 

0.012724 

resonator 

1.1715 

3.1129 

3. 6034 

26.10 

3.006546 

0.012618 

8.3075 

0.1239 

0.5336 

34.70 

0.036893 

0.313279 

0. 1039 

0.0286 

3.3776 

64.00 

0.008544 

0.313422 

y.223S 

0.0534 

3.8403 

57.80 

3.003503 

3.010372 

dual 

0.3678 

0.3717 

1.5253 

43.40 

0.003426 

3.310279 

stack 

0.5241 

0.0753 

2.3015 

37.23 

3.003314 

3.013142 

pLirs  he 

3.a507 

0.0u47 

3.3852 

26.90 

0.008250 

0.318063 

lOng 

0.4262 

0.0676 

4.2996 

41.20 

0. 003377 

0.310219 

resonator 

0.2737 

0.0564 

2.5689 

51.30 

0.303455 

0.018313 

0.1194 

0.0301 

1.0507 

60.50 

3.333530 

3.310486 

0.0452 

0.0124 

3.3311 

64.60 

0.308557 

0.010433 

0.37B3 

3.0289 

0.4738 

30.30 

3.311312 

3.313798 
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TH 

TC 

FREQ 

PMEAN 

PMEAN 

ypRESS 

Po/Pm 

VLOCKIN   VHEAT 

deg  C 

deg  C 

Hz 

rnVrms 

psia 

Vrms 

7. 

(iiVrras     V 

23.7 

-41.0 

217.92 

100.12 

150.18 

0.28365 

2.996 

144.7 

dual 

24.6 

-31.2 

222.07 

100.73 

151.10 

0.28308 

2.971 

154.9 

stack 

24.8 

-21.1 

226.16 

100.63 

150.95 

0.28335 

2.977 

168.1 

He-Xe 

25.6 

-11.4 

230.13 

101.49 

152.24 

0.28289 

2.947 

177.9 

long 

22.2 

-40.6 

217.91 

97.26 

145.89 

0. 18935 

2.058 

99.0 

resonator 

23.0 

-28.0 

223.30 

100.87 

151.31 

0.18890 

1,980 

104.0 

22.5 

-14.0 

228.81 

99.64 

149.46 

0. 18928 

2.009 

118.6 

24.4 

-41.0 

217.90 

99.70 

149.55 

0.28332 

3.005 

144.7 

23.8 

-36.6 

^O'O .  j6 

99.70 

149.55 

0.28328 

3.004 

170.3 

dual 

24.0 

-27.9 

242.61 

100.30 

150.45 

0.23351 

2.989 

177.1 

stack 

24.2 

-16.9 

247.39 

100.56 

150.84 

0.23333 

2.979 

186.6 

He-Xe 

24.4 

-7.0 

251.73 

101.53 

152.30 

0.23349 

2.952 

193.8 

short 

22.0 

-34.6 

239.37 

99.43 

149.15 

0. 18893 

2.009 

113.6 

resonator 

-23.8 

244.41 

100.44 

150.66 

0.18864 

1.986 

118.5 

22.1 

-11.9 

249.70 

101.03 

151.55 

0. 18907 

1.979 

125.9 

24.1 

"■■jO»  J 

238.77 

99.40 

149.10 

0.28244 

3.004 

170.6 

25.9 

-31.0 

545.83 

100.49 

150.74 

0.28347 

2.983 

350.1 

26.6 

554.24 

98.71 

148.07 

0.28264 

3.028 

374.1 

dual 

26.6 

-15.0 

562.40 

99.23 

148.85 

0.28305 

3.016 

386.1 

stack 

27.2 

-7.1 

570.22 

99.49 

149.24 

0.23270 

3.004 

400.1 

pure  He 

23.1 

-34.6 

541.73 

100.09 

150.14 

0.18935 

2.000 

224.7 

short 

23.4 

-22.9 

554.27 

101.09 

151.64 

0.13919 

1.979 

235.3 

resonator 

23.4 

-12.9 

564.27 

99.59 

149.39 

0. 18939 

2.011 

254.4 

23.9 

-3.0 

574.07 

100.11 

150.17 

0.18897 

1.996 

264.1 

23.6 

5.0 

581.68 

99.83 

149.75 

0.18951 

2.007 

276.7 

23.4 

-32.2 

543.70 

99.83 

149.75 

0.13909 

2.003 

227.8 
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IHEAT 

RHEAT 

ACC 

TC/TH 

QHEAT 

WORK 

COP 

COPR 

HEATLK 

inA 

ohms 

cg5E-3 

watt 

watt 

Q=heat8r 

cl7/carn 

watt 

0.00 

94.33 

0.7321 

0.0000 

3.975S 

0.0000 

0.0000 

0.3921 

dual 

100.00 

94. 4S 

57.46 

0.3126 

0.9448 

4.2465 

0.2225 

0.0513 

0.3331 

stack 

147.60 

94.60 

63.50 

0.8460 

2.0609 

4.6112 

0.4469 

0.0314 

3.27S2 

He-Xe 

179.50 

94.70 

68.39 

0.8762 

3.0513 

4,3708 

0.u264 

0.0385 

0.22-2 

long 

0.00 

94.39 

36.04 

0.7374 

0.0000 

1.3131 

0.0000 

0.0000 

0.3£0b 

resonator 

77.10 

94.52 

33.79 

0.3278 

0.5619 

1.9052 

0.2949 

0.0613 

0.3091 

121.10 

94.66 

45.33 

0.3765 

1.3882 

^a 1 /6j 

0.6379 

0.0893 

0.2212 

0.00 

94.33 

52.67 

0.7802 

0.0000 

3.9711 

0.0000 

0.0000 

0.3963 

0.00 

94.43 

67,36 

0.7966 

0.0000 

4.6704 

0.0000 

0.0000 

0.  :QtJ0 

dual 

69.00 

94.52 

71.77 

3,3253 

0.4500 

4.3600 

0.0926 

0.Qi9b 

0.3093 

stack 

106,40 

94.64 

77,11 

0,8613 

1.0714 

5.1161 

0.2094 

0.0336 

0.2450 

He-Xe 

130.30 

94.75 

31.49 

0,3945 

1.6037 

5.3154 

0.3026 

0.0357 

0.1871 

short 

0.00 

94,46 

45,42 

0.8082 

0.0000 

2.0313 

0.0000 

0.0000 

0.3373 

asonator 

57.90 

94,57 

48.38 

0.3440 

0,3170 

2.1676 

0.14.63. 

0.0270 

0.2743 

33.70 

94.70 

C-i  Ci 

0.3343 

0.6634 

2.3073 

0.2375 

0.0374 

0. 1026 

0.00 

94,43 

63.04 

0.7963 

0.0000 

4.6648 

0.0000 

0.0000 

0.3600 

0.00 

94,49 

319.20 

0.3097 

0.0000 

9.5091 

0.0000 

0.0000 

0.3653 

36.50 

94,57 

346.34 

0.8349 

0.7076 

10.1127 

0.0700 

0.0133 

0.3173 

Guai 

119.70 

94. 6o 

362.71 

0.8612 

1.3563 

10.4416 

0.1299 

0.0209 

0.2671 

stack 

143.60 

94.74 

381.09 

0.8853 

1.9536 

10.7941 

0.1810 

0.0233 

0.2202 

pure  He 

0.00 

94.45 

203.33 

0.3052 

0.0000 

4.1120 

0.0000 

0.0000 

3.3704 

short 

69,40 

94.57 

217.85 

0.3439 

0.4555 

4.2992 

0.1059 

0.0196 

0.2972 

resonator 

96,60 

94.67 

239.73 

0.3776 

0.8834 

4.6467 

0. 1901 

0.0265 

0.2330 

115.00 

94.83 

0.9094 

1.2541 

4.3096 

0.2608 

0.0260 

0.1727 

129.40 

94.90 

263.85 

0.9373 

1.5890 

5.0436 

0.3147 

0.0210 

0.1194 

0,00 

94.47 

206.88 

0.3125 

0.0000 

4.1621 

0.0000 

0.0000 

0.3570 
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C0P2 

C0PR2 

Qtot 

DELTA  T 

Vise  PEN  THERM  PEN 

w  cl9 

c20/carn 

watt 

deg 

cm 

cm 

0.0986 

0.0275 

0.3921 

64.70 

0.006744 

0.012999 

dual 

0.3021 

0.3697 

1.2S29 

55.30 

3.006671 

0.012358 

stack 

3.5073 

0.3924 

2.3391 

45.90 

0.006616 

0.312752 

He-he 

0.6725 

0.0951 

3.2755 

37.00 

0.006539 

0.012&04 

long 

3.2093 

0.0565 

0.3806 

62.30 

3.306826 

0.013156 

rsBonator 

0.4571 

0.0951 

9. 8709 

51.00 

0.006630 

0.012779 

0.7395 

0.1042 

1.6394 

36.50 

0.006584 

0.012691 

0.0993 

a.02si 

3.3963 

65.40 

0.S067&7 

0.013043 

0.0771 

3.219/ 

0.3600 

60.40 

3.0064q1 

3.312453 

dual 

3.1562 

0.2331 

0.7593 

51.90 

0. 006339 

0.012315 

stack 

8.^573 

0.0413 

1.3164 

41.10 

0.306321 

0.312134 

He-xe 

9.3379 

9.3399 

1.7958 

31.40 

0.006239 

0.312025 

ihori 

3.1621 

3.0335 

0.3373 

56.60 

0.006439 

0.012411 

resonator 

0.2733 

0.0505 

0.5913 

46.10 

0.906343 

9.012226 

3.3753 

9.0433 

3.3661 

34.00 

3.006255 

0.312357 

0.0772 

0.0197 

0.3600 

60.40 

0.006471 

3.012472 

0.0334 

0.0390 

0.3653 

56.90 

9.008122 

3.009907 

0.1014 

0.0201 

1.0254 

49.50 

0.338142 

0.009932 

dual 

0. 1555 

0.3251 

1.6234 

41.60 

3.003062 

3.309334 

stack 

0.2014 

3.0260 

2.1733 

34.39 

3.308004 

0.009763 

pure  He 

0.0901 

3.0213 

0.3704 

57.70 

3.308131 

0.009918 

short 

3.1751 

0.0324 

0.7527 

46.30 

9.005032 

0.009761 

resonator 

0.2403 

3.3335 

1.1165 

36.30 

0.007991 

0.009747 

0.2967 

0.0295 

1.4268 

26.90 

0.307903 

0.009647 

0.3384 

0.3226 

1.7035 

13.60 

0.007363 

3.009592 

0.0858 

3.0198 

0.3570 

55.  &0 

0.S0S131 

3.009913 
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